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PREFACE 


This  report  describes  the  results  of  the  development,  experimental  calibration  and 
predictions  of  a  physiologically  based  pharmacodynamic  model  simulating  biological  effects  of 
oxidative  stress  induced  by  chemicals  in  vitro  and  in  vivo.  This  is  one  of  a  series  of  technical 
reports  and  publications  describing  results  of  a  collaborative  effort  conducted  by  ManTech 
Environmental  Technology,  Inc.,  Toxic  Hazards  Research  Unit,  located  at  Wright-Patterson  Air 
Force  Base,  and  by  Occupational  and  Environmental  Health  Directorate,  Toxicology  Division, 
and  aimed  at  pharmacodynamic  description  of  biological  effects. 

The  animals  used  in  this  study  were  handled  in  accordance  with  the  principles  stated  in  the 
Guide  for  the  Care  and  Use  of  Laboratory  Animals  prepared  by  the  Committee  on  Care  and  Use 
of  Laboratory  Animals  of  the  Institute  of  Laboratory  Animal  Resources,  National  Research 
Council,  National  Academy  Press  1996,  and  the  Animal  Welfare  Act  of  1966,  as  amended. 

Research  performed  by  ManTech  Environmental  Technology  was  conducted  under 
Department  of  the  Air  Force  Contract  No.  F41624-96-C-9010.  Lt  Col  Terry  A.  Childress, 
Director  of  the  Toxicology  Division,  served  as  Contract  Technical  Monitor. 
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INTRODUCTION 


CHEMICALLY  INDUCED  OXIDATIVE  STRESS 

Pro-oxidant  chemicals  are  those  compounds  that  may  bring  about  a  state  of  excess  one- 
electron  oxidations,  either  directly  or  indirectly  via  metabolic  breakdown.  Oxidative  stress  is  a 
pathophysiological  process  in  which  the  balance  between  pro-oxidants  and  antioxidants  in  tissue 
is  shifted  towards  pro-oxidants  (Figure  1;  Byczkowski  and  Channel,  1996).  Pro-oxidant 
chemicals  may  be  provided  by  environmental,  occupational,  or  therapeutic  exposure  to 
xenobiotics  (Kehrer,  1993),  may  arise  from  dietary  polyunsaturated  fat  (Gower,  1988;  Finley 
and  Otterbum,  1993;  Haegele  et  al.,  1994),  or  may  be  produced  endogenously  during 
physiological  function  of  aerobic  cells  (Byczkowski  and  Gessner,  1988).  From  the  primary  pro¬ 
oxidant  chemical,  further  metabolic  reactions  generate  free  radicals  (defined  as  molecules  or 
groups  of  atoms  with  one  unpaired  electron),  and  then,  an  avalanche-type  process  (e.g.,  lipid 
peroxidation)  may  release  secondary  and  tertiary  free  radicals  (Figure  2;  for  review  see 
Roberfroid  and  Calderon,  1994).  Additional  factors,  such  as  aging  (Stadtman  et  al.,  1993)  or 
dietary  deficiencies,  may  augment  the  oxidative  stress  status.  Depletion  of  cellular  antioxidants, 
as  well  as  defective  enzymatic  scavenging  systems,  further  increase  oxidative  stress  and  may 
enhance  damage  to  cellular  components.  Oxidative  stress  can  be  reversed  by  natural  and 
synthetic  antioxidants  (Williams,  1993;  Papas,  1993;  Pratt,  1993).  Possible  definitions  and 
implications  of  oxidative  stress  and  the  basic  literature  on  its  measurement  were  reviewed  in  a 
recent  publication  by  Byczkowski  and  Channel  (1996). 

Several  chemicals  can  cause  an  oxidative  stress  directly  (e.g.,  by  generating  free  radicals 
during  metabolism  by  cytochrome  P450  (CYP)  or  during  redox  cycling;  Kappus,  1986)  or 
indirectly  (e.g.,  by  stimulating  a  respiratory  burst  in  inflammatory  cells;  Kulkami  and 
Byczkowski,  1994a).  In  addition,  several  natural  and  synthetic  peroxides  (e.g.,  tumor-promoting 
organic  hydroperoxides;  Taffe  et  al.,  1987;  Timmins  and  Davies,  1993)  can  be  cleaved  by  trace 
amounts  of  the  transition  metals,  directly  producing  highly  reactive  free  radicals  without 
involvement  of  enzymatic  metabolic  pathways  (Kulkami  and  Byczkowski,  1994b).  Chemically 
induced  oxidative  stress  causes  derangement  of  antioxidant  mechanisms  in  tissues  (Videla  et  al., 
1990),  may  lead  to  lipid  peroxidation  (Comporti,  1985),  inhibition  of  cellular  enzymatic 
activities  (e.g.,  CYP  activity;  Willis,  1980)  and  may  result  in  cell  injury  (de  Groot  and  Littauer, 
1989).  It  has  been  demonstrated  that  lipid  peroxidation  may  cause  necrotic  tissue  damage  rather 
than  be  an  effect  of  necrosis  (Biasi  et  al.,  1995). 
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Antioxidant  Defense  System 


Nonenzymatic  Enzymatic 


Figure  1.  Cellular  balance  between  pro-oxidants  capable  of  inducing  oxidative  stress  (’Ag02  - 

singlet  oxygen;  OH*  -  hydroxyl  radical;  R*  -  carbon-,  nitrogen-  or  oxygen-centered  free 
radical;  02*“-  superoxide  anion  radical;  -0-0*  -  peroxyl  radical;  H2O2-  hydrogen 
peroxide),  and  intracellular  antioxidants  (SOD  -  superoxide  dismutase;  Catalase;  GSH- 
Px  -  glutathione  peroxidase;  GSH  -  glutathione;  Tocopherol  -  vitamin  E;  Ascorbate  - 
vitamin  C;  |3-Carotene;  modified  from  Byczkowski  and  Channel,  1996;  according  to 
Bray  and  Betteger,  1990). 

LIPID  PEROXIDATION 

Lipid  peroxidation  is  a  pathological  process  leading  to  a  unique  form  of  hepatocellular  injury 
implicated  in  the  genesis  of  liver  necrosis  evoked  by  several  pro-oxidant  chemical  hepatotoxicants 
(e.g.,  carbon  tetrachloride  -  CCI4,  yellow  phosphorous,  ethanol,  etc.;  Kulkami  and  Byczkowski, 
1994a),  and  it  may  be  linked  with  carcinogenicity  (Byczkowski  and  Channel,  1996).  Lipid 
peroxidation  is  characterized  by  the  formation  of  conjugated  dienes,  formation  of  thiobarbitunc 
acid  reactive  substance  (TBARS;  mainly  malondialdehyde),  and  the  exhalation  of  alkanes  (e.g., 
ethane).  TBARS  and  ethane  are  generated,  among  other  stable  products,  during  the  propagation  and 
termination  of  lipid  peroxidation  process  (Figure  3).  The  alkanes  are  formed  in  biological  systems 
through  peroxidation  of  the  omega-3  (ethane)  or  omega-6  (pentane)  fatty  acids  and  the  subsequent 
beta-scission  decomposition  of  the  intermediate  hydroperoxides  (Gardner,  1989). 
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Pro-Oxidant  Chemicals  and  Free  Radicals  Involved  in  Oxidative 

Stress 


Pro-Oxidant  Chemicals 


Chemical  and 

Metabolic 

Generation 


Direct 


Indirect 


V7 

Primary  Free  Radicals 

Peroxyl  (ROO- ) 

Alkoxy!  (RO-) 

Phenoxyl  (0O  ■ )  1 

Alkyl  (RC-) 
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Activation  of 
Natural  Cellular 
Sources 


etc. 


Antioxidants 


Interactions  with  Biomolecules 

i  i  i 

Gene  and/or  Enzyme  DNA  Lipid 


Secondary  Free  Radicals 

And  Reactive  Oxygen  Species 
n  Superoxide  (Oj ) 

i/  Hydrogen  Peroxide  (HjCVj) 
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Tertiary  Free  Radicals 

n  Lipid  Peroxyl  (LOO-) 
. — p1.  Lipid  Alkoxy  I  (LO-) 

Activation/inactivation  Damage/Repair  Peroxidation  Lipienyi  (L-) 

I  I  | 

Biomarkers  of  Oxidative  Stress 


Figure  2.  Free  radicals  and  reactive  oxygen  species  involved  in  the  chemically  induced  oxidative 
stress  (modified  from  Byczkowski  and  Channel,  1996,  according  to  Trush  and  Kensler, 
1991). 

BIOMARKERS  OF  LIPID  PEROXIDATION 

Based  on  the  findings  of  Riely  et  al.  (1974)  that  mice  treated  with  carbon  tetrachloride  have  an 
increased  amount  of  the  exhaled  ethane  in  vivo,  numerous  studies  have  been  conducted  in  which 
ethane  and/or  pentane  were  measured  as  indices  of  lipid  peroxidation  or  surrogate  biomarkers  of 
tissue  damage  by  oxidative  stress.  Increased  ethane  exhalation  was  also  found  by  Cojocel  et  al. 
(1989)  as  a  consequence  of  lipid  peroxidation  in  mice  treated  with  trichloroethylene  (TCE). 
Consequently,  volatile  alkane  (ethane  and/or  pentane)  detection  in  expired  air  has  been  used  for 
some  time  as  a  non-invasive  technique  to  measure  lipid  peroxidation  in  whole  animals  or  human 
subjects  (Refat  et  al.,  1991;  Kazui  et  al.,  1992;  Arterbery  et  al.,  1994;  Guilbaud  et  al.,  1994).  Ethane 
exhalation  is  more  reliable  as  an  index  of  lipid  peroxidation  than  pentane,  because  CYP-mediated 
metabolism  of  ethane  is  substantially  slower  than  that  of  pentane  (Smith,  1991). 
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A  Scheme  of  Lipid  Peroxidation  Process 

\/=\M/  RC0°H 

v  v  v  Linolenic  Acid  (C18:3,w3) 


.  ^H: 


’  A  A  /RC00H 

\/ v/\/v  Conjugated  Diene  Radical 


°2^J 
•oo  * 


GSH 


R-COOH 

Peroxyl  Radical 


Endoperoxides 


RCOOH  | 

Hydroperoxide  XBARS 

Ogpx 

Me2+ ^  Hydroxy  Fatty  Acid 

9  ▼  0ROOOH 

AA^V  Alkoxyl  Radical 

k  O 

II  A  , 

/VVV 

'  Fatty  Acid  Aldehyde 
^3^2*  Ethy‘  Rac,ical 


Me  1^ — ►  CH_CH  Ethane 
▼  _  3  3  1 

Ethylene 

Figure  3.  A  simplified  scheme  of  lipid  peroxidation  reactions  which  lead  to  production  of  thio- 
barbituric  acid  reactive  substances  (TBARS)  and  to  generation  of  ethane  (according 
to  Sagai  and  Ichinose,  1980).  Me2+  -  transition  metal  cation. 

Typically,  in  the  assay,  exhaled  air  is  probed  for  ethane  by  gas  chromatography  using  flame 
ionization,  photoionization,  or  ion  trap  detectors  (Rneepkens  et  al.,  1994).  Because  of  its  direct 
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relation  to  lipid  peroxidation  (Jeejeebhoy,  1991;  Figure  3),  we  chose  ethane  exhalation  assay  as  an 
end  point  for  development  of  the  computer-aided  PBPD  model  for  simulation  of  the  biological 
effects  caused  by  CCI4,  TCE,  and  other  pro-oxidant  chemicals  (tert-butyl  hydroperoxide  and 
bromotrichloromethane,  BrCCl3). 

The  measurement  of  rates  of  TBARS  generation  in  vitro  was  used  previously  for  calibration  of 
the  biologically  based  pharmacodynamic  model  (BBPD)  of  lipid  peroxidation  induced  by  tert-butyl 
hydroperoxide  and  BrCCl3  in  mouse  liver  slices,  described  elsewhere  by  Byczkowski  et  al.  (1995; 
1996).  In  the  present  report,  a  physiologically  based  pharmacodynamic  (PBPD)  model  of  lipid 
peroxidation  is  described  which  employs  ethane  exhalation  as  a  measurable  end  point  in  vivo  and 
incorporates  the  previously  developed  in  vitro  model  as  a  mechanistic  module.  The  developed 
PBPD  model  has  been  linked  with  a  PBPK  sub-model  that  describes  the  local  concentrations  of 
CCI4  and  TCE  in  the  liver.  The  resultant  hybrid  PBPD  model  may  be  used  for  a  pharmacodynamic 
description  of  oxidative  stress,  dose-response  characterization,  and  risk  characterization  of  pro¬ 
oxidant  chemicals. 

RISK  CHARACTERIZATION  OF  PRO-OXIDANT  CHEMICALS  BASED  ON  THEIR 
MODE  OF  ACTION 

Health  risk  from  chemicals  depends  on  both  the  extent  of  exposure  and  a  dose-response 
relationship,  which  reflects,  in  turn,  the  mode  of  action  of  chemicals  (Figure  4). 

Risk  Characterization  Combines  Dose-Response  with  Exposure 


Figure  4.  A  paradigm  for  risk  characterization. 

Because  biological  responses  of  tissues  and  organs  are  mechanistically  linked  to  local 
concentration  of  the  active  form  of  xenobiotic,  the  internal  dose  of  chemical  that  reaches  a 
particular  physiological  compartment  must  be  used  for  any  meaningful  risk  characterization. 
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Without  understanding  "what  the  particular  dose  of  a  chemical  can  do  to  the  organism"  the  whole 
process  of  risk  characterization  is  useless  to  a  risk  manager  (Figure  4). 

In  addition  to  the  question  about  the  direct  biological  effect  of  the  delivered  concentration  of 
an  activated  (free  radical)  form  of  the  pro-oxidant  chemical,  other  important  questions  are:  i.  ”How 
is  the  target  organ  protected  against  free  radical  insult?”;  ii.  “How  fast  are  normal  (physiological) 
processes  of  autooxidation?”;  and  iii.  “How  fast  do  the  activating  enzymes  (CYP)  degrade?” 

For  quantitative  characterization  and  modeling  of  the  dose-response  for  pro-oxidant  chemicals 
(e.g.,  CCI4,  BrCCl3,  or  TCE)  it  was  necessary  not  only  to  answer  these  questions,  but  also  to 
include  the  mode  of  action  and  to  determine  the  exact  chain  of  events  in  the  chemical  interaction 
with  the  biological  system  (Figure  5). 


Conceptual  Framework  of  Oxidative  Stress  Modeling 


Figure  5.  A  conceptual  framework  for  quantitative  modeling  and  dose-response  characterization  of 
the  chemically  induced  oxidative  stress. 

At  first,  the  internal  delivered  dose  of  pro-oxidant  chemical  has  been  described  as  a  function  of 
exposure  with  an  appropriate  pharmacokinetic  (PK  and/or  PBPK)  model.  Then,  the  local 
concentration  of  free  radicals,  generated  by  the  pro-oxidant  chemical,  has  been  described  as  a 
function  of  the  local  dose  of  the  chemical  with  an  appropriate  BBPD  model.  Next,  using 
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mechanistic  information,  oxidative  stress  and  lipid  peroxidation  have  been  described  as  functions  of 
free  radical  concentration.  Finally,  the  health  effects  may  be  quantitatively  described  by 
biologically  based  dose-response  (BBDR)  sub-model  as  continuous  or  discrete  phenomena, 
depending  on  the  magnitude  of  oxidative  stress.  The  dependence  may  be  either  deterministic  or 
stochastic  in  nature.  These  considerations,  based  on  the  available  literature  and  experiments 
conducted  in  our  laboratory,  led  to  the  development,  calibration,  and  partial  verification  of  the 
hybrid  PBPD  model  for  chemically  induced  oxidative  stress.  The  resultant  computer-assisted 
simulation  model  described  in  this  report  may  be  useful  for  risk  characterization  of  pro-oxidant 
chemicals. 
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MATERIALS  AND  METHODS 


CHEMICALS 

All  chemicals  used  in  this  study  were  of  analytical  grade.  N-tert-butyl-a-nitrone  (PBN)  and 
2,2,5,5,-tetramethyl-l-pyrrolidinyloxyl-3-carboxyamide  (3-CAR)  were  purchased  from  Aldrich 
Chemical  Co.,  Inc.,  Milwaukee,  WI.  Dimethyl  sulfoxide  (DMSO),  BHT,  CCI4,  BrCCl3,  TCE, 
and  TBOOH  were  from  Sigma  Chemical  Co.,  St.  Louis,  MO. 

Dosage 

The  pro-oxidant  chemicals  tested:  carbon  tetrachloride  (CCI4),  bromotrichloromethane 
(BrCCl3),  trichloroethylene  (TCE),  and  tert-butyl  hydroperoxide  (TBOOH)  were  used  in  vitro  or 
in  vivo  in  the  following  concentrations  or  doses: 

Liver  slices  in  vitro  (recalculated  as  final  concentrations  in  the  medium): 

CCI4  -  0.1,  0.5, 1.0,  and  1.5  mM 
BrCCl3  -  0.1,  0.5,  1.0,  and  1.5  mM 
TCE  -  0.1,  0.5, 1.0,  and  1.5  mM 
TBOOH  -  0.1,  0.5,  1.0,  1.5  mM 

Mice  in  vivo  (single  i.p.  doses  expressed  per  body  weight): 

CCI4  -  0.075,  0.15,  0.3,  and  1.5  g/kg 
BrCCl3  -  0.025,  0.05,  0.1,  and  1.0  g/kg 
TCE-  0.26,  1.0,  and  2.6  g/kg 
TBOOH  -  0.025,  0.05,  0.1,  and  0.225  g/kg 

ANIMALS 

The  B6C3F1  male  mice  (Charles  River  Breeding  Laboratories,  Kingston,  NY),  25  -  38  g  body 
weight  were  used  throughout  the  study.  Routinely,  the  mice  were  provided  with  Purina  Formulab 
5008  standard  diet  and  Pseudomonas-free  softened  water  ad  libitum.  One  week  before  experiment, 
the  diet  was  switched  to  the  vitamin  A-  and  E-  deficient,  purified  diet  Purina  5827C-1 . 

MEASUREMENT  OF  LIPID  PEROXIDATION  IN  VITRO 

Precision-cut  slices  were  prepared  from  livers  of  B6C3F1  male  mice  and  maintained  using  the 
dynamic  roller  culture  method  (Sipes,  et  al.,  1987;  Brendel  et  al.,  1993).  The  mice  were  euthanized 
with  CO2,  their  livers  were  removed  and  placed  in  ice-cold  Sacks  buffer  (containing:  KH2PO4  0.75 
g/L,  K2HPO4  9.5  g/L,  NaHC03  1.2  g/L,  KHCO3  0.6  g/L,  mannitol  37.5  g/L,  and  MgCl2;  pH  7.4). 
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Liver  cores,  8  mm  diameter,  were  prepared  and  sliced  in  ice-cold  Sacks  buffer  using  a  Krumdieck 
tissue  slicer  (Alabama  Research  and  Development,  Munford,  AL;  Brendel  et  al.,  1987;  Krumdieck 
et  al.,  1980).  The  slices  were  loaded  on  rollers  (two  slices  per  roller)  in  ice-cold  Sacks  buffer.  The 
rollers  were  then  placed  in  scintillation  vials  containing  1.7  mL  of  Waymouths  MB  752/1  media  at 
37  °C  (Formula  78-5107EC,  without  phenol  red,  pH  7.4,  Gibco  BRL,  Grand  Island,  NY), 
supplemented  with  NaHCC>3  1.3  g/L,  HEPES  2.38  g/L,  NaCl  0.292  g/L,  1-glutamine  0.35  g/L, 
gentamycin  sulfate  50  mg/L,  and  10%  fetal  bovine  serum  (Hyclone,  Logan,  UT),  and  capped  with  a 
scintillation  vial  cap  with  1/4"  hole  for  gas  exchange.  The  vials  were  placed  in  a  Dynamic  Roller 
Culture  Incubator  (Vitron,  Tucson,  AZ)  and  gassed  with  95%  02/5%  CO2  for  a  2-h  preincubation 
period. 

After  a  2-h  preincubation  period,  the  rollers  were  removed  from  the  vials,  placed  into 
prewarmed,  sealed  vials  containing  fresh  media  (pH  7.4),  and  dosed  through  the  septa  with  either 
vapors  (for  volatile  chlorinated  hydrocarbons:  CCI4,  B1CCI3,  TCE)  or  an  appropriate  dilution  (for 
water-soluble  compound,  TBOOH)  of  pro-oxidant  chemical  at  the  desired  final  concentration.  The 
dosed  vials  containing  rollers  were  then  returned  to  the  roller  culture  incubator.  Final 
concentrations  of  volatile  pro-oxidant  chemicals  in  the  media  were  calculated  using  partition 
coefficients,  determined  for  the  equilibrated  medium/air  system  at  37  °C  (for  CCI4  0.666  \+  0.05 
S.D.,  n=15],  for  BrCCl3  1.97  [±  0.23  S.D.,  n=18],  for  TCE  1.94  [±  0.17  S.D,  n=18]).  Zero  time 
controls  were  processed  immediately.  Then,  the  vials  were  removed  at  intervals  over  a  2-h 
incubation  and  slices  were  weighed  and  sonicated  in  their  own  media.  Finally,  aliquots  of  each 
sonicated  sample  were  removed  for  TBARS  assay  and  protein  content  measurements.  Samples  for 
TBARS  assay  were  added  to  ice-cold  D-PBS/GSH/EDTA  buffer  (pH  7.4)  containing  20  mg 
reduced  GSH  and  48  mg  EDTA  in  100  mL  D-PBS  (Dulbeco's  buffer;  Gibco  BRL,  Grand  Island, 
NY).  Incubations  were  repeated  several  times  with  different  liver  preparation  (typically,  n=4). 

Lipid  peroxidation  was  measured  by  the  formation  of  TBARS,  employing  the  fluorescence 
spectrophotometry  of  solvent  tissue  extracts  (Janero,  1990).  Essentially,  in  this  assay,  the  aldehyde 
products  generated  by  splitting  the  endoperoxide  alkoxyl  radicals  (formed  during  the  peroxidation 
of  unsaturated  fatty  acids;  mostly  malondialdehyde,  MDA)  reacted  with  thiobarbituric  acid  (TBA) 
to  yield  a  1:2  MDA:TBA  red,  fluorescent,  complex  (Janero,  1990).  Incubation  of  liver  slices 
without  chemical  inducer  did  not  significantly  increase  the  fluorescence  for  up  to  2  hours.  Since  the 
control  values  for  liver  slices  at  time  zero  were  subtracted  from  the  results,  under  conditions  of  the 
assay,  the  determined  amount  of  MDA:TBA  complex  reflected  the  extra  amount  of  lipid 
hydroperoxides  produced  in  addition  to  the  normal  physiological  background. 
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At  1  hour  and  2  hours  of  incubation,  the  samples  of  liver  slices  (control  and  treated)  were 
removed  for  viability  analysis.  The  viability  was  assessed  from  lactate  dehydrogenase  (LDH), 
aspartate  aminotransferase  (AST),  and  alanine  aminotransferase  (ALT)  leakage,  and  intracellular 
potassium  content.  The  enzyme  leakage  was  determined  using  a  Kodak  Ektachem  Analyzer  (model 
700XR)  for  aminotransferase  activities  and  DuPont  acaV  for  dehydrogenase  activity.  An  acceptable 
enzyme  leakage  level  for  precision-cut  liver  slices  was  assumed  to  be  less  than  20%  of  the  total 
content  of  enzymatic  activity.  Potassium  content  in  sonicated  tissue  samples  was  determined  using 
an  AVL  982-S  Electrolyte  Analyzer  (Roswell,  GA).  The  acceptable  level  of  intracellular  potassium 
content  in  precision-cut  liver  slices  was  assumed  to  be  greater  than  35  mM  K+/g  wet  weight.  If  the 
average  viability  tests  of  either  control  or  treated  liver  samples  did  not  meet  the  above  acceptable 
levels,  the  experimental  results  were  discarded. 

Free  Radical  Measurement 

Known  amounts  of  a  spin  label  (N-tert-butyl-a-nitrone,  PBN)  were  added  to  liver  slices  in 
Waymouth’s  media  and  incubated  with  or  without  addition  of  appropriate  pro-oxidant  chemical 
(1  mM  CCI4,  BrCCl3,  TCE,  TBOOH,  for  60  min.).  The  radicals  generated  by  these  chemicals 
formed  adducts  with  PBN  which  were  detected  by  the  electron  paramagnetic  resonance  (EPR) 
spectroscopy  (Buettner,  1987).  The  total  radicals  in  the  lyophilized  samples  of  liver  slices  were 
measured  using  a  Bruker  EMS  104  EPR  analyzer.  The  machine  parameters  for  the  EPR  analyzer 
were:  microwave  power,  25  mW;  sweep  width,  100  G;  modulation  amplitude,  4.02  G;  sweep 
time,  10.49  s;  filter  time  constant,  20.48  ms;  receiver  gain,  60.  The  spectra  were  measured  by 
peak  height  directly  from  the  EMS  104  EPR  analyzer  and  by  double  integration  with 
normalization  for  receiver  gain  using  the  EPR  program  (Bruker,  Billerica,  MS). 

All  results  were  recalculated  per  liver  dry  weight  and  analyzed  by  one-way  and  two- 
factorial  analysis  of  variance  using  the  statistics  package  Design  Ease®.  The  factors  were 
concentration  and  time.  Standard  deviations  and  regression  correlation  were  performed  using 
Sigma  Plot®. 

ANIMAL  TREATMENT  AND  MEASUREMENT  OF  LIPID  PEROXIDATION  IN  VIVO 
Animal  Treatment 

Male  B6C3F1  mice  (body  weight  29  -  32  g),  fed  for  one  week  the  vitamin  A-  and  E-  deficient 
diet,  were  used  to  conduct  the  exhalation  experiments  and  partition  coefficient  determinations. 
Animals  were  treated  intraperitoneally  (i.p.)  with  pro-oxidant  chemicals  at  appropriate  doses 
(calculated  in  g/kg  body  weight),  dissolved  in  0.2  mL  of  mineral  oil.  Immediately  after  treatment, 
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ethane  production  was  monitored  using  a  closed  gas  uptake  system  (Gargas  et  al.,  1986).  Four  or 
five  mice  were  placed  simultaneously  in  a  0.75  L  chamber  containing  50  grams  of  soda  lime  that 
absorbs  CO2  and  H2O.  Chamber  oxygen  concentration  was  monitored  (MDA  oxygen  analyzer, 
MDA  Scientific,  Lincolnshire,  IL)  and  kept  at  a  range  of  20-2 1 .5%  throughout  the  exposure.  Ethane 
exhalation  measurements  were  repeated  several  times  with  different  groups  of  mice  before  and  after 
the  treatment  (typically,  n=4). 

Air  Sample  Analysis 

Samples  were  collected  using  an  automatic  sampling  valve  (1  mL  sample  loop)  connected  to  a 
Hewlett  Packard  5890  GC.  Ethane  was  separated  from  other  respiratory  gases  by  a  6’  x  1/8" 
stainless  steel  column  packed  with  Chromosorb  102,  80-100  mesh  (Alltech,  Deerfield,  IL).  The 
column  temperature  was  50  °C,  and  the  injector  and  flame  ionization  detector  temperatures  were 
125  °C  and  200  °C,  respectively.  A  carrier  gas  (nitrogen)  flow  rate  was  set  at  20  mL/min  and  the 
air  plus  hydrogen  flow  rate  was  405  mL/min.  The  retention  time  was  2.1  min  after  the  valve 
opened.  A  Hewlett-Packard  3396  Integrator  was  used  to  measure  peak  heights;  and  ethane 
concentrations  were  calculated  using  a  calibration  curve  prepared  with  ethane  standards. 

A  background  noise  at  low  levels  of  ethane  as  well  as  the  ethane  detection  threshold  and  peak 
integration  by  GC  and  data-processing  software  were  responsible  for  a  quantification  threshold, 
below  which  any  measurement  of  ethane  concentration  was  uncertain.  Due  to  this  uncertainty,  a 
reliable  ethane  quantification  level  by  the  method  used  was  above  approximately  0.025  ppm. 
(minimum  ethane  concentration  integrated  as  a  peak  by  GC  +  2  SD). 

Partition  Coefficients 

Partition  coefficients  for  ethane  were  determined  in  our  laboratoiy  (Seckel  and  Byczkowski, 
1996)  using  a  modified  vial-equilibration  method  in  vitro,  described  by  Gargas  et  al.  (1989). 
Tissues  from  five  B6C3F1  mice  were  pooled,  homogenized,  and  aliquoted  into  12.4  mL  headspace 
analyzer  vials.  Tissues  analyzed  included:  blood,  liver,  fat,  kidney,  and  muscle.  Each  tissue  sample 
weighed  1  gram.  To  inhibit  spontaneous  lipid  peroxidation,  butylated  hydroxytoluene  (BHT;  Sigma 
Chemical,  St.  Louis,  MO)  was  added  in  the  amount  of  5  mg/g  tissue.  A  15  ppm  concentration  of 
ethane  gas  was  added  to  the  headspace  of  each  vial.  The  vials  with  constituents  were  vortexed  for  3 
h  at  37  °C.  A  sample  was  removed  from  the  headspace  of  each  vial  via  a  Hewlett  Packard  19395A 
autosampler.  A  Hewlett  Packard  5890  GC  was  used  to  analyze  samples  with  data  handled  by  a  P.  E. 
Nelson  Data  Aquisition  System  equipped  with  Turbochrom  (version  4.0)  software.  A  Poraplot  Q, 
25  m  x  0.53  mm  (Chrompak,  The  Netherlands),  was  used  for  chemical  separation  along  with  the 
following  GC  conditions:  oven  temperature  100  °C,  injector  temperature  100  °C,  flame  ionization 


11 


detector  250  °C,  and  N2  carrier  flow  through  column  and  headspace  sampler  5.3  cc/min.  Ethane 
retention  time  was  determined  at  1 .58  min. 

MATHEMATICAL  MODELING  AND  COMPUTER-ASSISTED  SIMULATIONS 

Ethane  metabolism  parameters  (Vmax  and  Km)  were  estimated,  fitted,  and  optimized  with  the 
PBPD  model  which  was  based  on  our  in  vivo  experimental  measurements  of  ethane  uptake  in  a 
closed  chamber,  using  the  method  for  volatile  chemicals  essentially  as  described  by  Gargas  et  al. 
(1986).  Partition  coefficients  were  determined  using  the  method  of  Gargas  et  al.  (1989),  as 
described  above.  The  PBPD  model  was  written  in  Advanced  Continuous  Simulation  Language 
(ACSL;  Mitchell  and  Gauthier  Associates,  Inc.  1993)  with  a  sub-routine  written  in  FORTRAN. 
The  simulations  were  performed  using  SIMUSOLV®  software  with  optimization  capabilities 
(Steiner  et  al.,  1990)  on  a  VAX/VMS  minicomputer.  Parameters  were  optimized  by 
SIMUSOLV®,  which  uses  the  log  likelihood  function  as  the  criterion.  Either  the  generalized 
reduced  gradient  method  for  single  parameter  optimization  or  the  Nelder-Mead  search  method 
for  multiple  parameters  optimization  was  used  to  adjust  the  values  (Steiner  et  al.,  1990). 
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RESULTS 

MODEL  STRUCTURE 
Modules 

The  PBPD  model  for  chemically  induced  oxidative  stress  was  composed  of  several  modules, 
interlinked  to  form  integral  sub-models  (Figure  6).  Each  sub-model  was  calibrated  and  verified 
individually  with  experimental  data  from  our  laboratory  and/or  from  the  available  literature. 

Modular  Structure  of  PBPD  Sub-Model  for  Chemically  Induced 

Lipid  Peroxidation _ 


Ethane  Exhalation 


Figure  6.  A  simplified  scheme  of  physiologically  based  pharmacodynamic  (PBPD)  sub-model 
for  chemically  induced  lipid  peroxidation. 

The  model,  describing  biological  effects  within  the  target  organ,  was  constructed  in  a  way 
compatible  with  physiologically  based  pharmacokinetic  (PBPK)  module  and/or  classic 
compartmental  pharmacokinetic  (PK)  module  which  estimated  the  internal  dose  of  a  chemical.  The 
biologically  based  pharmacodynamic  (BBPD)  module  described  bioactivation  of  a  chemical 
(decreased  by  cytochrome  P450  (CYP)  degradation  and  suicidal  inhibition),  production  of  free 
radicals  (quenched  by  antioxidants),  and  lipid  peroxidation  (enhanced  by  autooxidation) 
accompanied  by  TBARS  and  ethane  generation.  The  PBPD  module  described  distribution  and 
exhalation  of  ethane. 
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PBPK  Sub-Model  for  Internal  Dose  of  Pro-Oxidant  Chemical 

The  PBPK  sub-model  for  internal  dose  of  pro-oxidant  chemicals  was  constructed  as  a 
traditional,  flow-limited  PBPK  mathematical  description  of  volatile  compounds  (Gargas  et  al., 
1986;  Yang  and  Andersen,  1994).  A  separate  input  compartment,  parallel  to  the  gastrointestinal 
tract,  was  added  to  this  sub-model  to  estimate  intraperitoneal  exposure.  The  rate  of  change  in  the 
amount  of  chemical  absorbed  from  the  peritoneal  cavity  (RAjp  [mg/h])  was  described  as  a 
product  of  the  rate  constant  of  absorption  (KAip  [1/h])  and  the  mass  of  chemical  remaining  in 
the  peritoneal  cavity  (MRjp  [mg]): 

RAip  =  KAip  *  MRip 

and  t 

MRip  =  Dip +  jRMRip*dt 

where  Dip  [mg/animal]  is  the  actual  dose  of  the  chemical  injected  i.p.,  RMRjp  [mg/h]  is  the  rate 
of  change  of  the  chemical  remaining  in  the  peritoneal  cavity  (disappearance),  and  t  [h]  is  time. 

It  was  estimated  that  only  a  small  fraction  of  the  lipophilic  pro-oxidant  chemical  may  diffuse 
directly  to  the  abdominal  mesenteric  fat,  whereas  about  99%  of  the  chemical  is  being  absorbed  to 
venous  circulation  (blood  absorption  ratio,  Bab  =  0.99  [ratio])  and  eventually  drained  to  the  portal 
blood  and  delivered  to  the  liver.  Therefore,  the  product  of  Bab  *  RAjp  was  added  to  the  liver 
compartment. 

The  PBPK  sub-model  parameters  for  CCI4,  based  on  Paustenbach  et  al.  (1988)  and  Gallo  et  al. 
(1993),  were  scaled  allometrically  to  B6C3F1  mice  and  calibrated  with  data  from  the  literature 
(Seckel  and  Byczkowski,  1996).  Similarly,  the  PBPK  sub-model  parameters  for  TCE,  based  on 
Fisher  et  al.  (1991),  were  scaled  allometrically  to  B6C3F1  mice  and  calibrated  with  data  from  the 
literature  (Das  et  al.,  1994).  The  PBPK  sub-model  parameters  for  BrCCl3  and  TBOOH  were  not 
verified  experimentally. 

Classic  PK Module  for  Intraperitoneal  Dosing  of  Pro-Oxidant  Chemicals 

For  many  chemical  compounds  the  PBPK  sub-model  parameters  are  not  immediately 
available,  whereas,  their  classic  pharmacokinetic  micro-  and  macro-constants  may  occasionally  be 
found  in  the  literature.  To  utilize  this  kind  of  data,  an  additional  classic  PK  module  was  included  in 
the  pharmacokinetic  sub-model  and  connected  directly  to  the  liver  compartment.  The  classic  PK 
module  was  based  on  a  multiexponential  equation  for  two-compartment  system: 
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Cl  =  Clo  *  kl,0  *  (exp(-(B  *  t)-exp(-a  *  t))/(a  -  P) 

where  Cl  is  local  concentration  of  pro-oxidant  chemical  [mg/kg];  Clo  is  estimate  of  initial 
concentration  of  the  chemical  [mg/kg];  kj  O  is  pharmacokinetic  transfer  micro-constant  [1/h];  a  is 
pharmacokinetic  macro-constant  [1/hr];  and  p  is  pharmacokinetic  macro-constant  [1/h].  The 
pharmacokinetic  constants  were  recalculated  from  the  literature  (e.g.,  for  CCI4  as  presented  by 
Seckel  and  Byc2kowski,  1996). 

BBPD  Module  for  Activation  of  Pro-Oxidant  Chemical  and  Free  Radical  Concentration 

Production  of  free  radicals  and  the  local  concentration  of  pro-oxidant  chemical-derived  free 
radicals  were  estimated  by  the  square  root  algorithm  presented  previously  (Byczkowski  and 
Flemming,  1996)  and  verified  with  TCE  (Byczkowski  et  al.,  1996;  Channel  et  al.,  1997).  For  the 
bioactivation  and  free  radical  quenching  reactions: 

ki  kt 

Pro-oxidant  chemical  ->  FRa(j  +  FRacj  -»  Nonradical  products 

it  was  assumed  that:  dFR/dt  =  kj  *  Clm  -  kt  *  FRacj  *  FRa(j  =  0. 

Therefore: 

FRad  =  Vkj  *  Clm^ 

where  FRa(j  is  a  steady  state  concentration  of  pro-oxidant  chemical-derived  free  radicals  [pmol/0.1 
g  liver];  Clm  is  local  pro-oxidant  chemical  concentration  [pmol/0.1  g  liver];  kj  is  a  rate  constant  of 
free  radical  formation  from  the  pro-oxidant  chemical  [1/h];  and  kt  is  the  lumped  rate  constant  of 
free  radical  recombination  and  quenching  by  the  biological  system  [1/h]. 

Algorithm  for  CYP  degradation  and  suicidal  inhibition  by  activated  pro-oxidant  chemical  was 
based  on  deterministic  one-hit  mechanism  according  to  the  reaction  (presented  previously  by 
Byczkowski  and  Flemming,  1996): 


ACR 

FRad  +  CYP  ->  INACTIVE  CYP 

it  was  assumed  that  different  kinds  of  CYP  taking  part  in  the  bioactivation  of  pro-oxidant  chemical 

have  uniform  sensitivity  to  FRad- 

Therefore: 

ACtrem  =  AC  *  exp(-ACR*  FRad  *  t) 
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where  ACtrem  is  concentration  of  active  CYP  remaining  over  time  [pmol/0.1  g  liver];  AC  is  the 
initial  concentration  of  active  CYP  [pmol/0.1  g  liver];  ACR  is  the  rate  constant  of  CYP  inactivation 
by  free  radicals  [1/h];  t  is  time  of  incubation  with  free  radicals  [h];  and  FRa(j  is  a  steady-state 
concentration  of  pro-oxidant  chemical-derived  free  radicals  [pmol/0.1  g  liver]. 

BBPD  Sub-Model  for  Lipid  Peroxidation 

The  BBPD  sub-model  estimated  activity  of  lipid  peroxidation  expressed  as  TBARS  production 
and  ethane  generation  (Figure  7).  The  BBPD  module  for  lipid  peroxidation  in  the  liver  (calibrated 
with  TBARS  production),  confirmed  with  experimental  data  for  TBOOH  and  BrCCl3  in  precision- 
cut  mouse  liver  slices,  was  published  elsewhere  by  Byczkowski  et  al.  (1996).  A  simplified  scheme 
of  this  module  is  shown  in  Figure  7  (depicted  by  ovals).  The  source  codes  of  *.CSL  and  *.CMD 
files  for  this  module  are  archived  in  PBPK-L  Public  Domain  Source  Library  and  are  accessible 
through  the  World  Wide  Web  at  the  following  URL:  http://www.navy.al.wpafb.af.mil/new.htm 

The  PBPD  module  for  ethane  exhalation  (depicted  by  rectangles  in  Figure  7)  was  calibrated 
with  experimental  data  from  our  laboratory  (presented  previously  by  Seckel  and  Byczkowski, 
1996).  The  PBPD  module  estimated  rates  of  metabolism,  distribution,  and  exhalation  of  ethane 
generated  in  the  liver.  On  a  molecular  basis,  about  0.1%  of  hydroperoxides  derived  from  natural 
lipids  will  decompose  to  yield  ethane  (efficiency  of  ethane  generation  from  fat,  EFfe  <  0.001 
[molar  ratio];  Gardner,  1989;  Janero,  1990). 

BBDR  Sub-Model  for  Cellular  Target  Inhibition 

An  inhibition  of  activities  of  cellular  targets,  caused  by  free  radicals,  was  estimated  by  a 
biologically  based  dose-response  (BBDR)  sub-model  composed  of  two  modules,  deterministic 
and  stochastic  (presented  previously  by  Byczkowski  and  Flemming,  1996).  This  sub-model  was 
governed  by  dose-dependent  algorithms  with  time  of  exposure  to  free  radicals  tp  [h]  fixed  as  a 
fraction  (Ft  [ratio])  of  the  time  needed  to  reach  the  maximum  effect  (Tme  [h]).  For  independent 
variable,  the  initial  local  concentrations  of  pro-oxidant  chemical  Co  [pM]  were  increased  stepwise 
between  the  estimated  maximum  “no  effect”  dose  (Cmjjj  [pM])  and  minimum  100%  effect  dose 
(Cmax  [pM])  with  the  interval  Ca  [pM]  resulting  in  a  number  of  iterations  “i”.  The  initial  local 
concentrations  Co  were  captured  as  an  array  DOSCj .  The  local,  steady  state  concentrations  of  pro¬ 
oxidant  chemical-derived  free  radicals  FRacj]y]  [pM]  (FRadM  =  10  *  FRacj),  at  fixed  time  of 
exposure  tp  [h],  were  captured  as  an  array  DOSj  along  with  initial  local  concentrations  of  pro- 
oxidant  chemical  Co  [pM]. 
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PBPD  Sub-Model  for  Lipid  Peroxidation 


cone  CEP 


Figure  7.  An  interlinked  biologically  based  pharmacodynamic  (PBPD)  sub-model  describing  lipid 
peroxidation  combined  the  BBPD  module  for  TBARS  production  (ovals)  with  the  PBPD 
module  for  ethane  exhalation  (rectangles). 

CONE  -  concentration  of  inhaled  ethane  [ppm];  QP  -  alveolar  ventilation  rate  [L/h];  CEV  -  mixed 
venous  blood  ethane  concentration  [mg/L];  CEP  -  concentration  of  exhaled  ethane  [ppm];  QC  - 
cardiac  output  [L/h];  CEA  -  concentration  of  ethane  in  arterial  blood  [mg/L];  QR  -  blood  flow  to 
rapidly  perfused  tissues  [L/h];  QS  -  blood  flow  to  slowly  perfused  tissues  [L/h];  QF  -  blood  flow  to 
fat  tissue  [L/h];  QL  -  blood  flow  to  liver  tissue  [L/h];  QG  -  blood  flow  through  the  portal  vein  [L/h]; 
CEVG  -  concentration  of  ethane  in  the  portal  vein  [mg/L];  AEM  -  amount  of  metabolized  ethane 
[mg];  CEVL  -  concentration  of  ethane  in  venous  blood  leaving  the  liver  [mg/L];  PXREM  - 
accumulated  remaining  hydroperoxides  [  pm  o  1/0.1  g];  TBARS  -  thiobarbitunc  acid  reactive 
substances  [pmol/0.1  g];  TPX  -  accumulated  total  hydroperoxides  [pmol/0.1  g];  ANOX  -  vitamin 
E-type  antioxidants  [pmol/0. 1  g];  AUTOX  -  hydroperoxides  produced  by  autooxidation  [umo  1/0.1 
g];  PUF  -  polyunsaturated  fat  [pmol/0.1  g];  ACTIND  -  activated,  free  radical  form  of  chemical 
inducer  [pmol/0.1  g];  IND1  -  internal  dose  of  chemical  inducer  1  [pmol/0.1  g];  ACR  -  activator 
(CYP)  loss  rate  [l/h];  CL  -  delivered  dose  of  pro-oxidant  chemical  [mg/kg]. 
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Deterministic  Module 

The  deterministic  module  estimated  inhibition  of  uniform  cellular  targets  (In  [ratio])  by  pro¬ 
oxidant  chemical-derived  free  radicals  (FRadM)  as  an  exponential  decay  function: 

kd 

FRadM  +  CELLULAR  TARGET  ->  INACTIVE  TARGET 
it  was  assumed  that  homogenous  CELLULAR  TARGETS  have  uniform  sensitivity  to  FRad- 

In  =  l0  *  exp(-kd  *  FRadM  *  tp) 

where  In  is  a  remaining  activity,  expressed  as  a  fraction  of  remaining  active  cellular  targets,  relative 
to  the  amount  before  inhibition  [ratio];  I()  is  the  initial  concentration  of  active  cellular  targets, 
assumed  to  be  100%  (Io  =  1-  [percentage/100]);  kd  is  the  rate  constant  of  cellular  target  inactivation 
by  free  radicals  [100%/pM/h];  tp  is  time  of  exposure  to  free  radicals  [h];  and  FRadM  are  the  local, 
steady  state  concentrations  of  pro-oxidant  chemical-derived  free  radicals  [pM],  The  values  of  the 
remaining  relative  activity  In  were  captured  as  an  array  Inh: 

Inh  =  Pnb  In2>  —  Ini]- 


Stochastic  Module 

The  stochastic  module  estimated  inhibition  of  non-uniform  cellular  targets  (PROB  [ratio])  by 
exposure  for  time  tp  [h]  to  pro-oxidant  chemical-derived  free  radicals  (array  DOS)  as  a  1  minus 
time-weighted  fraction  of  a  cumulative  Gaussian  distribution  function: 

FRadM  +  MULTIPLE  TARGETS  ->  RANGE  OF  INHIBITORY  RESPONSES 
it  was  assumed  that  non-uniform  MULTIPLE  TARGETS  give  normal  distribution  of  the 
INHIBITORY  RESPONSES  to  FRad- 

PROB  =  1  -  Ft  *  I  1/V2rc  *  exp(-z2/2)  *  dz 
-  00 

where  x  =  (DOS  -  M)/SD;  PROB  is  probability  of  cellular  targets  to  remain  active,  relative  to  the 
amount  before  inhibition  (expressed  as  a  fraction  of  remaining  active  cellular  targets,  relative  to  the 
amount  before  inhibition  [ratio]);  DOS  is  array  of  free  radical  concentration  values  [pM],  generated 
during  exposure  to  the  range  of  concentrations  of  pro-oxidant  chemical  (between  maximum  no 
effect”  and  minimum  “100%  effect”  doses);  Ft  is  as  a  fraction  [ratio]  of  the  time  needed  to  reach 
the  maximum  effect  (Ft  =  tp/Tme  );  M  is  mean  of  the  cumulative  Gaussian  distribution  of  free 
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radical  concentration  values  [pM];  SD  is  standard  deviation  of  the  cumulative  Gaussian  distribution 
of  free  radical  concentration  values  [pM];  and  z  is  a  variable  of  integration. 

PARAMETRIZATION  AND  CALIBRATION  OF  SUB-MODELS  WITH  DATA 
PBPK  Sub-Model  for  Internal  Dose  of  Pro-Oxidant  Chemical 

Chemical-dependent  parameters  for  TCE  PBPK  module  in  B6C3F1  mice  (confirmed  with 
data  from  Fisher  et  al.,  1991)  were  presented  previously  by  Das  et  al.  (1994),  and  those 
parameters  for  CCI4  (confirmed  with  data  from  Sanzgiri  et  al.,  1995;  Gallo  et  al.,  1993;  Gargas  et 
al.,  1986;  and  Paustenbach  et  al.,  1986;  1988),  scaled  to  B6C3F1  mice,  were  presented 
previously  by  Seckel  and  Byczkowski  (1996).  These  parameters,  optimized  with  SIMUSOLV®, 
are  listed  in  Table  1 .  The  animal-specific  PBPK  modeling  parameters  were  as  recommended  by 
ILSI,  RSI  (1994)  for  mice. 


TABLE  1 .  CHEMICAL-SPECIFIC  PHARMACOKINETIC  PARAMETERS  FOR  TCE  AND 
CCI4  IN  B6C3F1  MICE 


Parameter  •  Description 

TCE 

CCI4 

PBC 

Partition  coefficients  [ratio] 

Blood/air 

13.4e 

4.52a 

PLC 

Liver/blood 

2.03e 

3.14a 

PFC 

Fat/blood 

41.3e 

79.4a 

PRC 

Rapidly  perfused  tissue/blood 

2.03e 

3.14a 

PSC 

Slowly  perfused  tissue/blood 

1.0e 

2.43b 

MW 

Molecular  weight  [g/mol] 

131.5 

153.82 

VMAXC 

Metabolism  constants 

Maximum  velocity  [mg/hr/kg] 

33 .0e 

0.65C 

KM 

Michaelis-Menten  constant  [mg/L] 

0.25e 

0.25c 

KFC 

1st  order  rate  constant  [1 /hr/kg] 

2.4e 

0.0C 

KAIP 

Absorption  rate  [1/hr] 

First  order  i.p.  uptake 

1.0d 

1.45d 

klO 

Pharmacokinetic  transfer  constants  fitted  [1/hr] 
Micro-constant 

0.03d 

0.3d 

a 

Macro-constant 

0.01d 

1.5d 

3 

Macro-constant 

2.0d 

1.6d 

a  from  Gargas  (1988). 
b  from  Evans  et  al.  (1994). 
c  from  Gargas  et  al.  (1986) 
d  from  Seckel  and  Byczkowski  (1996). 
e  from  Das  et  al.  (1994). 
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Classic  PK  Module  for  Intraperitoneal  Dosing  of  Pro-Oxidant  Chemicals 

The  pharmacokinetic  micro-  and  macro-constants  [1/h]  were  fitted  to  data  from  our 
laboratory  obtained  with  B6C3F1  mice  treated  i.p.  with  TCE  and  CCI4  (Seckel  and  Byczkowski, 

1996) .  The  pharmacokinetic  transfer  constants  for  classic  PK  module  are  listed  in  Table  1. 

BBPD  Module  for  Activation  of  Pro-Oxidant  Chemical  and  Free  Radical  Concentration 

Two  algorithms  describing  the  relationship  between  steady-state  concentration  of  free  radicals 
and  local  TCE  concentration  were  tested,  the  square  root  algorithm  (assuming  that  two  free  radicals 
are  formed  from  one  molecule  of  TCE,  as  described  above)  and  a  linear  algorithm  (Byczkowski  et 
al.,  1996;  Channel  et  al.,  1997).  The  linear  algorithm  described  the  reaction  in  which  one  free 
radical  is  formed  from  one  molecule  of  TCE: 

ki  kt 

TCE  -»  FRacj  -»  Nonradical  products 

Assuming  steady  state  concentrations  of  free  radicals:  dFR/dt  =  kj  *  Clm  -  kt  *  FRad  = 

FRad  =  kj  *  Cutf/kt 

where  FRad  is  a  steady-state  concentration  of  TCE-derived  free  radicals  [pmol/0.1  g  liver];  Clm  is 
local  TCE  concentration  [pmol/0.1  g  liver];  kj  is  a  rate  constant  of  free  radical  formation  from  TCE 
[1/h];  and  kt  is  the  lumped  rate  constant  of  free  radical  recombination  and  quenching  by  the 
biological  system  [1/h]. 

Quantitative  measurements  of  FRa(j  in  vitro  using  an  EPR-spin  trapping  method  failed  to 
confirm  either  algorithm  (Figure  8).  Despite  a  relatively  large  variability  in  time-dependent 
measurements  of  TBARS  produced  in  liver  slices  incubated  with  TCE  (Figure  9),  the  dose- 
dependent  data  were  much  better  fitted  with  the  square  root  algorithm  (Figure  10,  curve  B)  than 
with  the  linear  algorithm  (Figure  10,  curve  A).  Thus,  the  quantitative  measurements  of  TBARS  in 
vitro  confirmed  an  adequate  description  of  the  relationship  between  concentration  of  free  radicals 
and  local  TCE  concentration  by  the  square  root  algorithm  (Byczkowski  et  al.,  1996;  Channel  et  al., 

1997) . 

BBPD  Sub-Model  for  Lipid  Peroxidation 

As  the  square  root  algorithm  and  the  time-dependent  activator  degradation  equation  were 
introduced  into  the  BBPD  module  for  lipid  peroxidation  in  the  liver  (originally  calibrated  as  a 
linear  algorithm  with  the  experimental  data  from  precision-cut  mouse  liver  slices  with  TBOOH  and 
B1CCI3,  Byczkowski  et  al.,  1996),  the  modeling  parameters  had  to  be  recalculated  and  optimized 
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with  SIMUSOLV®  software.  All  optimized  parameters  added  or  changed  from  those  originally 
published  by  Byczkowski  et  al.  (1996)  are  listed  in  Table  2.  The  chemical-specific  parameters 
(factors  ACTDGF  and  PTIND)  were  estimated  and  optimized  for  a  range  of  concentrations  for 
TCE  (0.4  -5.6  p.mol/0.1  g  liver)  and  CCI4  (0.5  -  7.6  pmol/0.1  g  liver). 


BBPD  Module:  Effect  of  TCE  on  Generation  of  Free  Radicals 


01234 


TCE  [jjmol/0.1  g  liver] 

Figure  8.  Calibration  of  the  algorithm  describing  concentration  of  free  radicals  under  steady-state 
conditions  (FRAD  [pmol/g])  with  experimental  data  from  Steel-Goodwin  et  al.  (1995) 
for  free  radical  generation  by  different  concentrations  of  TCE  [pmol/0.1  g  liver]  in 
mouse  liver  slices  using  an  EPR/spin-trapping  method. 

FRAD  -concentration  of  PBN-reactive  free  radicals  [pmol/g  liver].  The  squares  represent  actual 
average  experimental  data  points  (after  subtraction  of  physiological  background  levels  of  free 
radicals  produced  in  the  absence  of  TCE).  The  continuous  lines  are  computer-generated  simulations 
involving:  A  -  linear  algorithm;  B  -  square  root  algorithm. 

Since  the  original  BBPD  sub-model  for  lipid  peroxidation  was  calibrated  with  BrCCl3 
(Tappel  et  al.,  1989;  Byczkowski  et  al.,  1996),  we  have  used  the  same  well  known  pro-oxidant 
chemical  to  check  if  it  will  lead  to  a  measurable  lipid  peroxidation  at  comparable  local 
concentrations  both  in  vitro  and  in  vivo. 
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BBPD  Module:  Effects  of  TCE  on  Lipid  Peroxidation 


Figure  9.  Results  of  time-dependent  simulation  of  lipid  peroxidation  in  mouse  liver  slices  induced 
by  1  mM  TCE.  TBARS  -  thiobarbituric  acid  reactive  substance  x  10*3  [mmol/0.1  g  liver]. 
Small  squares  depict  average  experimental  data  from  our  laboratory  (n=4)  described  by 
Byczkowski  et  al.  (1996).  The  continuous  lines  depict  computer  simulations  with  BBPD  sub-model 
involving  the  square  root  algorithm  and  parameters  optimized  by  SIMUSOLV®  software  (listed  in 
Table  2),  amount  of  TBARS  at  time=0  was  subtracted  from  the  data. 


TABLE  2.  PHARMACODYNAMIC  PARAMETERS  FOR  TCE  AND  CCI4  IN  B6C3F1  MICE 


Parameter 

Description 

Optimized  numerical  value 

Factors  [1/p.mol] 

ACTDGF1 

Activator  degradation  factor  1  (TCE) 

0.0014 

ACTDGF2 

Activator  degradation  factor  2  (CCI4) 

1.75 

PTIND1 

Potency  of  inducer  1  (TCE) 

250. 

PTIND2 

Potency  of  inducer  2  (CCI4) 

4408. 

Rate  constants  [1/h] 

AUTOXF 

Autooxidation  rate 

0.00013 

PXREDF 

Hydroperoxide  reduction  rate 

0.17 

ACR 

Activator  degradation  rate 

0.025 

INDLF 

Inducer  loss  rate 

0.0001 

22 


BBPD  Module:  Effects  of  TCE  on  Lipid  Peroxidation 


Figure  10.  The  results  of  dose-dependent  simulations  of  lipid  peroxidation  in  mouse  liver  slices  for 
0.5  h  induced  by  different  concentrations  of  TCE  [pmol/0.1  g  liver]. 

TBARS  -  thiobarbituric  acid  reactive  substance  x  10‘3  [mmole/0.1  g  liver].  Small  squares  depict 
average  experimental  data  from  our  laboratory  (n=4).  Lines  depict  computer  simulations  with 
BBPD  sub-model  involving:  A  -  linear  algorithm  (PT1ND  1=6.9);  B  -  square  root  algorithm 
(PTIND  1=250).  The  other  parameters  were  optimized  by  SIMUSOLV®  software.  Amounts  of 
TBARS  in  untreated  controls  were  subtracted  from  the  data  (Byczkowski  et  al.,  1996). 

Bromotrichloromethane  stimulated  TBARS  generation  by  precision-cut  mouse  liver  slices 
in  vitro  (Figure  1 1).  The  lowest  measurable  effect  was  noticed  at  0.1  mM  BrCCl3  in  the  medium. 
Considering  its  partition  coefficients  (medium/air  =  2,  and  liver/air  =  29),  the  100  pM  BrCCl3 
could  produce  an  initial  concentration  of  0.145  p.mol/0.1  g  liver.  The  dose-response  curves  had  a 
characteristic  sigmoidal  shape  with  a  plateau  reached  at  about  1  mM  BrCCl3,  when  incubated  for 
0.5  h,  and  at  about  0.5  mM  BrCCl3  when  incubated  for  1  h.  A  prolonged  incubation  (above  1  h) 
with  higher  than  0.5  mM  concentrations  of  BrCCl3  failed  to  increase  further  the  rate  of  TBARS 
generation  by  liver  slices  (results  not  shown  here). 

Figure  12  shows  the  results  of  measurement  of  ethane  exhalation  by  mice  in  a  closed  gas 
chamber.  The  treatment  with  BrCCl3  at  a  dose  above  0.025  g/kg  resulted  in  increased  ethane 
exhalation  measured  lh  after  the  exposure.  Twice  as  high  dose  was  required  to  produce  a 
significant  increase  in  ethane  exhalation  after  0.5  h  from  the  exposure.  Again,  the  dose-response 
curves  had  a  sigmoidal  shape. 
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Effect  of  BrCCh  on  TBARS  Production 


Figure  11.  Effects  of  different  doses  of  BrCCl3  (concentration  in  the  medium  [mM])  on 

production  of  thiobarbituric  acid  reactive  substances  (TBARS  -  thiobarbituric  acid 
reactive  substance  x  10‘3  [mmole/0.1  g  liver])  by  mouse  liver  slices  incubated  for 
either  0.5  or  1  h  in  the  presence  of  inducer. 

The  respective  background  values  were  subtracted  from  data  points.  Background  production  of 
TBARS  in  controls  were  respectively:  0.0226  pmole/g  liver  (±  0.0025  S.D.,  n=4)  at  0.5  h,  and 
0.0303  pmole/g  liver  (±0.0035  S.D.,  n=4)  at  1  h.  Data  points  were  significantly  different  from 
the  corresponding  controls  (at  p<0.05,  n=4)  by  Student’s  t-test. 


Figure  13  shows  actual  time  courses  of  stimulated  ethane  exhalation  by  mice  treated  with 
four  different  doses  of  BrCCl3  (0.025,  0.05,  0.1,  and  1.0  g/kg).  The  reliable  ethane  quantification 
threshold  was  0.025  ppm  There  was  no  significant  difference  between  ethane  exhalation 
(followed  for  up  to  2  h)  in  mice  before  and  after  i.p.  injection  of  0.2  mL  of  mineral  oil  only  (by 
one-way  Anova  at  p>0.05). 
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Effect  of  BrCCIa  on  Ethane  Exhalation 


0  0.5  1 

BrCCI3  [grtcgfl 

Figure  12.  Effects  of  different  doses  of  BrCCl3  (injected  i.p.  [g/kg])  on  ethane  exhalation  by 
five  mice  during  either  0.5  or  1  h  from  the  exposure  to  inducer  (Ethane  -  concentra¬ 
tion  in  the  chamber  [ppm]). 

The  respective  physiological  ethane  exhalation  values  (background)  from  the  same  group  of  five 
mice,  measured  before  treatment,  were  subtracted  from  the  data  points  (respectively,  0.0435  ppm 
[±0.022  S.D.,  n=5]  at  0.5  h,  and  0.0447  ppm  [±0.020  S.D.,  n=4]  at  1  h). 

Verification  of  BBPD  Sub-Model for  Lipid  Peroxidation  in  vitro 

The  BBPD  module  of  lipid  peroxidation  sub-model  was  further  verified  with  in  vitro  effects  of 
CCI4  on  the  generation  of  TBARS  by  precision-cut  mouse  liver  slices  (Figure  14).  The  parameters 
were  optimized  with  SIMUSOLV®  to  satisfy  all  set  (time-  and  dose-dependent)  of  experimental 
data. 

Verification  of  BBPD  Sub-Model for  Lipid  Peroxidation  in  vivo 

Based  on  the  “upside  down”  PBPK  model  for  ethane  distribution  (input  in  the  liver,  output 
in  the  lung),  a  PBPD  module  for  ethane  exhalation  was  constructed,  assuming  that  the  production 
of  reactive  free  radical  metabolites  of  xenobiotics  takes  place  in  the  liver  (Figure  7).  At  first,  the 
module  was  calibrated  with  data  from  mice  in  vivo,  inhaling  a  known  concentration  of  ethane  in 
a  closed  gas  chamber  (Figure  15  presented  previously  by  Seckel  and  Byczkowski,  1996).  These 
calibrations  allowed  us  to  verify  partition  coefficients  of  ethane  and  to  estimate  metabolism 
constants  (Vmexc,  Kem,  and  Kefc).  Pharmacokinetic  parameters  for  ethane  exhalation  module 
are  listed  in  Table  3. 
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Time-Dependent  Effects  of  BrCCIa  on  Ethane  Exhalation 


Figure  13.  Time-course  of  effects  of  four  different  doses  of  BrCCl3  (injected  i.p.  [g/kg])  on 
ethane  exhalation  (concentration  in  the  chamber  [ppm])  by  five  mice. 

The  respective  physiological  ethane  exhalation  values  (background)  from  the  same  group  of  five 
mice,  measured  before  treatment,  were  subtracted  from  the  data  points.  The  horizontal  line 
depicts  a  reliable  detection  level  of  ethane  by  the  method  used  (detectibility  threshold  +  2  SD). 
The  ethane  exhalation  time  courses  after  exposure  to  0.05  (after  1  h  but  not  0.5  h),  0.1  and  1.0  g 
BrCCl3/kg  were  significantly  different  from  the  corresponding  control  curves  before  treatment 
(by  one-way  Anova  at  p<0.05).  Other  details  are  the  same  as  in  Figure  12. 

Next,  the  PBPD  sub-model,  including  the  module  for  ethane  exhalation,  was  tested  versus 
experimental  data  from  mice  treated  in  vivo  with  TCE  (Figure  16).  The  square  root  algorithm  was 
included  in  the  BBPD  module  for  free  radical  concentration  and  the  PBPD  sub-model  simulated 
time-dependent  effects  of  different  doses  of  TCE  on  ethane  exhalation  in  B6C3F1  mice  in  vivo 
(Figure  16).  The  reliable  ethane  quantification  threshold  was  0.025  ppm. 
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BBPD  Module:  Effects  of  CCI4  on  Lipid  Peroxidation 


Figure  14.  Time-dependent  effects  of  four  concentrations  of  CCI4  (D-0.1,  11-0.5,  H-1.0,  and  B-1.5 
mM)  on  TBARS  generation  by  precision-cut  mouse  liver  slices  in  vitro.  Lines  are 
computer-generated  simulations  by  the  lipid  peroxidation  sub-model  (optimized 
parameters:  PTIND1=4408;  ACR=0.025;  other  parameters  the  same  as  in  Table  2). 


TABLE  3.  PHARMACOKINETIC  PARAMETERS  FOR  ETHANE  EXHALATION  IN 
B6C3F1  MICE 


Parameter  Description  Numerical  value 


PLA 

Partition  coefficients  [ratio] 

Liver/air 

0.828 

PGA 

Gut/air 

0.996 

PFA 

Fat/air 

2.444 

PSA 

Slowly  perfused  tissue/air 

0.979 

PRA 

Richly  perfused  tissue/air 

0.996 

PEB 

Blood/air 

1.305 

MW 

Molecular  weight  [g/mol] 

30.0 

VMEXC 

Optimized  metabolism  constants 

Maximum  velocity  [mg/hr/kg] 

0.286 

KEM 

Michaelis-Menten  constant  [mg/L] 

0.51 

KEFC 

1st  order  rate  constant  [1 /hr/kg] 

2.786 

EFFE 

Efficiency  of  ethane  generation  [molar  ratio] 

0.001 

Data  presented  previously  by  Seckel  and  Byczkowski  (1996). 
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PBPD:  Closed  Chamber  Ethane  Gas  Uptake 


ID 
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Figure  15.  The  results  of  PBPD  sub-model  simulations  of  our  experimental  data  from  mice  inhaling 
ethane  in  a  closed  chamber  (initial  concentration  of  ethane  1  ppm). 

CEP  -  Concentration  of  ethane  in  a  closed  chamber  [ppm];  Time  -  [h].  Symbols  depict  experimental 
data  collected  from  the  closed  chamber  (four  animals  per  time-point);  continuous  line  is  the  PBPD 
model  simulation  (Seckel  and  Byczkowski,  1996). 

At  this  step,  the  yield  of  ethane  generation  from  lipid  hydroperoxides  was  estimated  (EFfe  =  0.001 
[molar  ratio]  )  by  comparing  the  simulated  molar  amounts  of  ethane  in  vivo  (Figure  16)  with 
TBARS  in  vitro  (Figure  10),  produced  in  response  to  the  same  total  cumulated  dose  (area  under 
the  concentration  curve,  AUC)  of  TCE  in  the  liver. 

Finally,  the  calibrated  PBPD  sub-model  was  verified  with  our  data  for  ethane  exhalation, 
induced  by  four  different  doses  of  carbon  tetrachloride  (CCLj-specific  parameters  were  taken  from 
our  in  vitro  calibration.  Figure  14).  At  this  verification  step,  the  parameters  estimated  during  the  in 
vivo  calibration  with  TCE  were  not  further  adjusted.  Figure  17  shows  the  time-dependent  effects  of 
four  different  doses  of  CCI4  on  ethane  exhalation  along  with  the  computer-generated  simulations 
by  the  BBPD  sub-model.  The  dose-dependent  effects  of  CCI4  on  ethane  exhalation  and  the 
computer-generated  simulations  at  three  different  times  from  treatment  are  presented  in  Figure  18. 
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PBPD:  Effects  of  TCE  on  Ethane  Exhalation  in  Mice 

TCE 


Figure  16.  The  results  of  PBPD  model  simulations  of  our  experimental  data  from  mice  treated  i.p. 

with  three  different  doses  of  TCE  (□-0.26,11-1.0,  and  B-2.6  g/kg)  and  exhaling  ethane 
in  a  closed  chamber  (Byczkowski  et  al.,  1996). 

CEP  -  concentration  of  ethane  in  a  closed  chamber  [ppm];  Time  -  [h].  Symbols  depict  experimental 
data  collected  from  the  closed  chamber  (five  animals  per  time-point)  with  the  respective 
physiological  ethane  exhalation  values  (background)  from  the  same  group  of  five  mice, 
measured  before  treatment,  subtracted  from  the  data  points.  Lines  are  the  PBPD  model 
simulations.  The  horizontal  line  depicts  a  reliable  ethane  quantification  threshold  by  the  method 
used  (detectibility  threshold  +  2  S.D.).  Inset  shows  PBPK  sub-model  simulations  of  local  TCE 
concentrations  in  the  liver.  CL  -  local  concentration  [mg/kg  liver];  T  -  time  [h]. 


BBDR  Sub-Model  for  Cellular  Target  Inhibition 

Two  modules,  composing  the  BBDR  sub-model  for  cellular  target  inhibition  by  pro-oxidant 
chemical-derived  free  radicals,  were  based  on  either  deterministic  or  stochastic  equation, 
respectively,  and  were  calibrated  with  in  vitro  data  from  the  literature  (Byczkowski  and 
Flemming,  1996). 

Deterministic  Module 

The  deterministic  BBDR  module  was  calibrated  with  experimental  data  of  Vroegop  et  al., 
(1995)  for  inhibition  of  amino  acid  and  glucose  transporters  in  N  18  neuronal  hybridoma  cell  line  in 
culture,  incubated  with  cumene  hydroperoxide  (Figure  19  A)  or  hydrogen  peroxide  (Figure  19  B), 
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PBPD:  Effects  of  CCL  on  Ethane  Exhalation  by  Mice  In  Vivo 


Time  [h]  Time  [h] 

Figure  17.  Time-dependent  effects  of  four  different  doses  of  CCI4  i.p.  on  the  concentration  of 

ethane  exhaled  by  mice  in  a  closed  gas  chamber  in  vivo  (A.  IU-0.075,  B.O-O.15,  Cfi- 
0.3,  and  D.  B-1.5  g/kg). 

Lines  are  computer-generated  simulations  by  the  ethane  exhalation  BBPD  sub-model  (parameters 
that  describe  lipid  peroxidation  were  the  same  as  in  Table  3).  The  respective  physiological  ethane 
exhalation  values  (background)  from  the  same  groups  of  five  mice,  measured  before  treatment, 
were  subtracted  from  the  data  points.  The  reliable  ethane  quantification  threshold  by  the  method 
used  was  0.025  p.p.m.  (detectibility  threshold  +  2  SD). 

respectively,  Heffetz  et  al.  (1990),  and  Hecht  and  Zick  (1992),  for  inhibition  of  protein  tyrosine 
phosphatase  (PTyrPase)  in  rat  hepatoma  cells  in  culture  incubated  with  hydrogen  peroxide  (Figure 
19  C)  or  vanadate  (Figure  19  D),  respectively.  The  rate  constant  for  quenching  of  free  radicals  by 
all  biological  systems  investigated  (kt  [1/pM/h])  was  estimated  as  200.0  [1/pM/h]  and  was  fixed 
during  all  simulations.  Rate  constants  of  free  radical  formation  (kj  [1/pM/h])  were  chemical- 
specific  and  rate  constants  of  cellular  target  inhibition  [1/p.M/h]  were  biological  system-specific, 
and  they  were  optimized  with  SIMUSOLV®  software.  Initial  concentrations  of  active  targets  were 
assumed  to  be  100%  in  uninhibited  biological  systems  before  incubation  with  pro-oxidant 
chemicals  (Iq  =  1.0  [%/100]).  Timing  parameters  [h]  were  dependent  on  experimental  setup.  The 


final  simulation  parameters  are  listed  in  Table  4,  as  presented  previously  by  Byczkowski  and 
Flemming  (1996). 

BBPD  Model:  Dose-Dependent  Effects  of  CCI4on  Ethane 


Dose  of  CCI4[g/kg] 

Figure  18.  Dose-dependent  effects  of  CCI4  i.p.  on  the  concentration  of  ethane  exhaled  by  mice  in  a 
closed  gas  chamber  in  vivo,  measured  at  three  different  times  from  treatment  ([3-0.33,8!- 
0.75,  and  1-1.0  h). 

Lines  (A,  B,  and  C)  are  computer-generated  simulations  by  the  ethane  exhalation  BBPD  sub-model 
(parameters  that  describe  lipid  peroxidation  were  the  same  as  in  Table  3).  The  respective 
physiological  ethane  exhalation  values  (background)  from  the  same  groups  of  five  mice, 
measured  before  treatment,  were  subtracted  from  the  data  points.  The  reliable  ethane 
quantification  threshold  by  the  method  used  was  0.025  ppm  (detectibility  threshold  +  2  SD). 

Stochastic  Module 

The  stochastic  BBDR  module  was  calibrated  with  experimental  data  of  Vroegop  et  al.  (1995) 
for  inhibition  of  amino  acid  transporter  and  mitochondrial  activity  in  N  18  neuronal  hybridoma  cell 
line  in  culture,  incubated  with  6-OH  dopamine  (Figure  20  A,  C)  or  hydrogen  peroxide  (Figure  20  B, 
D),  respectively.  Optimized  simulation  parameters  are  listed  in  Table  4,  as  presented  previously  by 
Byczkowski  and  Flemming  (1996). 

Some  data  from  the  literature  covered  a  dose-response  characteristics  of  the  wide  range  of  pro¬ 
oxidant  chemical  concentrations,  especially  when  the  response  kinetics  was  slow.  For  such  data 
sets,  to  discriminate  between  the  targeted  and  random  hit  mechanisms  using  the  BBDR  sub-model, 


it  was  necessaiy  to  express  the  pro-oxidant  chemical  concentrations  on  a  logarithmic  scale  rather 
than  on  a  linear  one  (Figure  21). 


BBDR:  Effect  of  Cum.OOH  on  AA  Transporter  BBDR:  Effect  of  H2Q2  on  Glucose  Transporter 


Figure  19.  Calibration  of  the  deterministic  BBDR  module  with  in  vitro  data  for  inhibition  of 
cellular  targets  by  free  radicals  generated  by  pro-oxidant  chemicals.  Lines  are  computer¬ 
generated  simulations  by  the  deterministic  BBDR  module.  Black  squares  depict  data 
from  the  literature.  Insets  show,  for  comparison,  simulations  by  the  stochastic  BBDR 
module  (curves  depicted  by  empty  squares)  with  the  same  data  (black  squares). 

Optimized  parameters  are  listed  in  Table  4. 

A.  Effects  of  different  concentrations  of  cumene  hydroperoxide  (Cum.OOH  [pM])  on  activity  of 
amino  acid  (AA)  transporter  in  N  18  neuronal  hybridoma  cell  line  in  culture  (INH  [ratio])  after  1  h 
(plus  1  h  preincubation  with  Cum.OOH;  data  points  from  Vroegop  et  al.,  1995).  B.  Effects  of 
different  concentrations  of  hydrogen  peroxide  (H2O2  [pM])  on  activity  of  glucose  transporter 
(Glue)  in  N  18  neuronal  hybridoma  cell  line  in  culture  (INH  [ratio]),  after  1  h  (plus  1  h 
preincubation  with  H2O2;  data  points  from  Vroegop  et  al.,  1995).  C.  Effects  of  different 
concentrations  of  hydrogen  peroxide  (H2O2  [pM])  on  activity  of  protein  tyrosine  phosphatase 
(PTyrPase)  in  rat  hepatoma  cells  in  culture  (INH  [ratio])  after  5  min  (plus  20  min  preincubation 
with  H2O2;  data  points  from  Heffetz  et  al.,  1990).  D.  Effects  of  different  concentrations  of 
vanadate  (Vanad  [p.M])  on  activity  of  protein  tyrosine  phosphatase  (PTyrPase)  in  rat  hepatoma  cells 
in  culture  (INH  [ratio]),  after  8  min  (plus  0.5  h  preincubation  with  vanadate;  data  from  Hecht  and 
Zick,  1992). 
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TABLE  4.  PARAMETERS  FOR  THE  BIOLOGICALLY  BASED  PHARMACODYNAMIC 
MODEL  DESCRIBING  THE  INHIBITION  OF  CELLULAR  TARGETS  BY 
CHEMICALLY  GENERATED  FREE  RADICALS 


Parameter  Description 

Numerical  value 

ki 

Rate  constants  of  FR  formation  [1/pM/h] 
from  Cum.OOH 

100.0a 

ki 

from  60HD 

200.0a 

ki 

from  H2O2,  vanadate  and  pervanadate 

18.0a,f 

ki 

from  TCE 

-O 

O 

O 

O 

O'N 

kt 

Rate  constant  of  FR  quenching  [1/pM/h] 
fixed  for  all  biological  systems 

200.0C 

kd 

Rate  constants  of  cellular  target  inhibition  [1/pM/h] 

AA  transporter 

0.075a 

kd 

mitochondria 

0.053 

kd 

Glue,  transporter 

O.ia 

kd 

PTyrPase 

i.of 

10 

Initial  concentration  of  active  targets  [%/100] 

assumed  value  for  all  biological  systems 

1.0 

tp 

Time  of  pre-incubation  with  FR  [h] 

Fig  19  A-B 

l.Od 

tp 

Fig.  8  B 

0.333e 

tp 

Fig.  19  D 

0.58 

tp 

Fig.  19  C 

0.333^ 

TSTOP 

Length  of  experiment  [h] 

Fig.  19A-B 

2.0d 

TSTOP 

Fig.  8  B 

0.333C 

TSTOP 

Fig.  19  D 

0.6338 

TSTOP 

Fig.  19  C 

0.42^ 

Abbreviations:  Cum.OOH  -  cumyl  hydroperoxide;  6-OHD  -  6-hydroxy  dopamine;  TCE  - 
trichloroethylene;  FR  -  free  radicals;  AA  -  amino  acids;  Glue.  -  glucose;  PTyrPase  -  protein  tyrosine 
phosphatase. 

a  fitted  to  Vroegop  et  al.  (1995). 
b  fitted  to  Steel-Goodwin  et  al.  (1995). 

c  estimated  from  Vroegop  et  al.  (1995),  Steel-Goodwin  et  al.  (1995),  Heffez  et  al.  (1990),  and 
Hecht  and  Zick  (1992). 
d  experimental  from  Vroegop  et  al.  (1995). 
e  experimental  from  Steel-Goodwin  et  al.  (1995). 
f  fitted  to  Heffez  et  al.  (1990)  and  Hecht  and  Zick  (1992). 
g  experimental  from  Hecht  and  Zick  (1992). 
h  experimental  from  Heffez  et  al.  (1990). 
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BBDR:  Effect  of  6-OH  Dopamine  on  AA  Transporter  BBDR:  Effect  of  H2Q2  on  AA  Transporter 


Figure  20.  Calibration  of  the  stochastic  BBDR  module  with  in  vitro  data  for  inhibition  of  cellular 
targets  by  free  radicals  generated  by  pro-oxidant  chemicals.  Curves  depicted  by  empty 
squares  are  computer-generated  simulations  by  the  stochastic  BBDR  module.  Black 
squares  depict  data  from  the  literature.  Insets  show,  for  comparison,  simulations  by  the 
deterministic  BBDR  module  (lines)  with  the  same  data  (black  squares).  Optimized 
parameters  are  listed  in  Table  4. 

A.  Effects  of  different  concentrations  of  6-hydroxy  dopamine  (6-OH  Dopamine  [pM])  on  activity 
of  amino  acid  (AA)  transporter  in  N  18  neuronal  hybridoma  cell  line  in  culture  (PROB  [ratio]), 
after  1  hour  (plus  1  h  preincubation  with  6-OH  dopamine;  data  points  from  Vroegop  et  al.,  1995). 

B.  Effects  of  different  concentrations  of  hydrogen  peroxide  (H2O2  [pM])  on  activity  of  amino  acid 
transporter  (AA)  in  N  18  neuronal  hybridoma  cell  line  in  culture  (PROB  [ratio]),  after  1  hour  (plus 
1  h  preincubation  with  H2O2;  data  points  from  Vroegop  et  al.,  1995).  C.  Effects  of  different 
concentrations  of  on  activity  of  mitochondrial  reduction  of  MTT  (Mitochondria)  in  N  18  neuronal 
hybridoma  cell  line  in  culture  (PROB  [ratio]),  after  1  hour  (plus  1  hr  preincubation  with  60HD; 
data  points  from  Vroegop  et  al.,  1995).  D.  Effects  of  different  concentrations  of  hydrogen  peroxide 
(H2O2  [pM])  on  mitochondrial  reduction  of  MTT  (Mito)  in  N  18  neuronal  hybridoma  cell  line  in 
culture  (PROB  [ratio]),  after  1  hour  (plus  1  hr  preincubation  with  H2O2;  data  points  from  Vroegop 
etal.,  1995). 


For  some  data  of  Vroegop  et  al.  (1995),  the  simulations  with  BBDR  sub-model  could  not 
distinguish  between  the  targeted  or  random  hit  mechanisms,  and  both  deterministic  and  stochastic 
BBDR  modules  gave  an  equal  “goodness  of  fit”  (Figure  22). 
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BBDR:  Effect  of  PerVanadate  on  PTyrPase 


BBDR:  Effect  of  Cum.OOH  on  Mitochondria 


Figure  21.  Calibration  of  the  BBPD  sub-model  with  in  vitro  data  for  inhibition  of  cellular  targets 
by  free  radicals  generated  by  pro-oxidant  chemicals.  Lines  are  computer-generated 
simulations  by  the  deterministic  BBDR  module  (A)  whereas  curves  depicted  by  empty 
squares  are  computer-generated  simulations  by  the  stochastic  BBDR  module  (B).  Black 
squares  depict  data  from  the  literature.  Insets  show,  for  comparison,  simulations  by 
either  stochastic  (A)  or  deterministic  (B)  BBDR  modules,  with  the  same  data  (black 
squares).  Optimized  parameters  are  listed  in  Table  4. 

A.  Effects  of  different  concentrations  of  pervanadate  (Pervan.  [Log  pM])  on  the  activity  of  protein 
tyrosine  phosphatase  (PTyrPase)  in  rat  hepatoma  cells  in  culture  (INH  [ratio]),  after  15  min  (plus 
0.5  h  preincubation  with  Pervan.;  data  points  from  Heffez  et  al.,  1990).  B.  Effects  of  different 
concentrations  of  cumene  hydroperoxide  (Cum.OOH  [Log  pM])  on  activity  of  mitochondrial 
reduction  of  MTT  (Mitochondria)  in  N  18  neuronal  hybridoma  cell  line  in  culture  (PROB  [ratio]), 
after  1  h  (plus  1  h  preincubation  with  60HD;  data  points  from  Vroegop  et  al.,  1995).  Optimized 
parameters  are  listed  in  Table  4. 

BBDR:  Effect  of  Cum.OOH  on  Glucose  Transporter 


Figure  22.  Calibration  of  the  BBPD  sub-model  with  effects  of  different  concentrations  of  cumene 
hydroperoxide  (PCONC  [pM  Cum.OOH])  on  activity  of  glucose  transporter  (GLU)  in 
N  18  neuronal  hybridoma  cell  line  in  culture  (INH  [ratio]),  after  1  h  (plus  1  h  pre¬ 
incubation  with  Cum.OOH;  data  points  from  Vroegop  et  al.,  1995). 

The  continuous  line  is  a  computer-generated  simulation  by  the  deterministic  BBDR  module.  Black 
squares  depict  data  from  the  literature.  Insets  show,  for  comparison,  a  simulation  by  the 
stochastic  BBDR  module  (the  curve  depicted  by  empty  squares)  with  the  same  data  (black 
squares).  Optimized  parameters  are  listed  in  Table  4. 
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DISCUSSION 


PBPD  MODEL 

Several  mathematical  models  of  lipid  peroxidation  and  free  radical  reactions  leading  to  the 
oxidative  stress  were  developed  and  described  in  the  literature  (Babbs  and  Steiner,  1990;  Suzuki 
and  Ford,  1994;  Antunes  et  al.,  1994;  Vroegop  et  al.,  1995).  However,  none  of  these  models  could 
be  applied  for  dose-response  characterization  of  pro-oxidant  chemicals  in  vivo.  The  purpose  of  this 
modeling  effort  was  to  provide  a  quantitative  tool,  based  on  biochemical  mechanisms,  capable  of 
predicting  the  biological  effects  of  pro-oxidant  chemicals  (Figure  23). 


PBPD  Model:  Interlinked  BBPD/PBPK  Description  of  Chemically 
-Induced  Oxidative  Stress 


INPUT  SUB-MODEL  FUNCTION 


OUTPUT 


O  Pro-Oxidant 

q  Chemical _ 

^  Exposure 


PBPD 


BIOMARKER 

e.g.,  Ethane 
exhalation 


Local  concentrations  of  active 
metabolite  in  target  tissue 


Dosimetry  of  free  radicals  generated 
in  target  tissue  by  pro-oxidant 


Co 


FR 


Testing  for  quantitative  relation 
between  dose  and  response 

Testing  for  mode  of  action 
deterministic  vs.  stochastic 
and  predict  dose-dependent  distri¬ 
bution  of  remaining  active  targets 


BIOMARKER 

INH  or 
PROB 


Quantification  of  biological  RELATIVE 

response  depending  on  remaining  ACTIVITY 
or  unaffected  activity  of  cellular  targets 


ESTIMATE  OF 
BIOLOGICAL  EFFECT 


RESPONSE 


Figure  23.  A  scheme  of  PBPD  model  of  chemically  induced  oxidative  stress.  For  verification  of  the 
quantitative  PBPD  time-  and  dose-response  model  for  pro-oxidant  chemicals,  a  continu¬ 
ously  responding  end-point  was  measured  as  a  biomarker  (lipid  peroxidation  which 
generates  TBARS  and  ethane). 

PBPK  -  physiologically  based  pharmacokinetic  sub-model;  BBPD  -  biologically  based 
pharmacodynamic  sub-model;  PBPD  -  physiologically  based  pharmacodynamic  module;  BBDR  - 
biologically  based  dose-response  sub-model. 

The  resultant  PBPD  model  of  chemically  induced  oxidative  stress,  consisting  of  three  sub¬ 
models  (PBPK/PK,  BBPD,  BBDR)  and  a  PBPD  module  for  ethane  exhalation,  may  be  used  to 
predict  time-  and  dose-response  of  biological  systems,  both  in  vitro  and  in  vivo,  to  pro-oxidant 
chemicals,  such  as  volatile  chlorinated  hydrocarbons  (e.g.,  TCE,  CCI4,  or  BrCCl3). 
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Several  crucial  assumptions  were  made  during  the  construction  of  this  PBPD  model.  The 
PBPK  sub-model  was  built  assuming  the  blood-flow  limited  delivery  of  pro-oxidant  chemicals 
(Evans  and  Andersen,  1995;  Kedderis,  1996).  This  was  justified  in  the  case  of  volatile  chlorinated 
hydrocarbon  solvents  (Fisher  et  al.,  1991;  Gallo  et  al.,  1993;  Gargas  et  al.,  1986;  Paustenbach  et  al., 
1988),  and  if  needed,  the  sub-model  may  be  expanded  with  algorithms  involving  “diffusion 
coefficients”  for  diffusion-limited  chemicals  (e.g.,  xenobiotics  similar  to  dioxin;  Kohn  et  al.,  1993). 
The  PBPK  sub-model  may  be  fitted  with  other  sets  of  chemical-specific  parameters  and  may  be 
scaled  allometrically  to  simulate  local  delivery  of  the  pro-oxidant  chemical  in  different  animal 
species  (Yang  and  Andersen,  1994).  The  PK  module  can  accept  pharmacokinetic  macro-  and 
micro-constants  for  two-compartment  classical  model. 

The  BBPD  sub-model  was  built  on  a  template  of  a  previously  published  mathematical  model 
for  chemically  induced  lipid  peroxidation  in  precision-cut  liver  slices  (Byczkowski  et  al.,  1996), 
assuming  a  steady-state  concentration  of  pro-oxidant  chemical-derived  free  radicals  in  the  liver 
over  time  at  each  dose  level.  This  assumption  does  not  hold  during  the  initial  phase  of  free  radical 
generation/mixing  and  during  the  terminal  phase  of  free  radical  metabolism/quenching  (Figure  24 
A)  and,  therefore,  it  could  not  be  used  for  rapid  pulses  of  free  radicals.  On  the  other  hand,  under 
experimental  conditions,  the  free  radical  reactions  and  diffusion  rates  are  very  fast  (Antunes  et  al., 
1994),  so  the  equilibrium  in  the  biological  system  is  achieved  almost  instantaneously  and,  thus,  the 
steady-state  approximation  gave  satisfactory  estimates  of  free  radical  biological  effects  in  cellular 
systems  (McKenna  et  al.,  1991;  Vroegop  et  al.,  1995).  Consequently,  the  local  dose  of  free  radicals 
in  the  liver  was  considered  as  an  array  (DOS)  of  sustained  steady-state  concentrations  over  time  at 
each  pro-oxidant  chemical  delivered  dose  (Figure  24  B)  ignoring  both  the  initial  and  the  terminal 
phase  (assumed  a  set  of  “zero  order”  functions  of  FR  overtime;  Figure  24). 
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Steady  State  Concentration  of  Pro-Oxidant  Chemical-Derived 

Free  Radicals  An  Array  of  Free  Radical  Steady  State  Concentrations 


Figure  24.  Representations  of:  A.  a  steady  state  concentration  of  pro-oxidant  chemical-derived 
free  radicals  (FRad).  B.  an  array  (DOS)  of  the  steady  state  concentrations  of  pro-oxidant  chemical- 
derived  free  radicals  (FRad)- 

The  PBPD  module  was  built  as  an  “upside  down”  PBPK  model  for  a  volatile  compound 
(ethane),  assuming  its  production  in  the  liver  only,  and  the  blood-flow  limited  redistribution  (Seckel 
and  Byczkowski,  1996).  Obviously,  other  metabolically  competent  organs  may  also  generate  free 
radicals  from  pro-oxidant  chemicals  and  contribute  to  the  ethane  production,  but  since  the  liver  is 
the  most  active  metabolically  (Kulkami  and  Byczkowski,  1994a),  its  contribution  is  overwhelming. 

Once  the  maximum  no-effect  concentration  of  pro-oxidant  chemical  in  the  liver  (Cmjji)  and 
minimum  100%  effective  concentration  (Cmax)  have  been  estimated  (along  with  the  time  needed 
to  accomplish  the  maximum  effect,  Tme)  using  the  PBPK  and  BBPD  sub-models,  the  BBDR  sub¬ 
model  was  used  to  simulate  and  predict  the  dose-dependent  inhibitions  of  cellular  targets.  Because 
in  the  BBDR  sub-model,  delivered  doses  of  pro-oxidant  chemical  in  the  liver  (expressed  as  initial 
concentrations,  Co)  were  set  as  a  discrete  variable  (with  an  interval  C^),  the  corresponding  steady- 
state  concentrations  of  free  radicals  over  time  were  returned  by  the  model  as  an  array  (DOS;  Figure 
24  B).  Considering  the  two  possible  modes  of  action  of  free  radicals  on  cellular  targets, 
deterministic  or  stochastic,  the  frequency  of  hitting  targets  depended  on  specificity  and 
homogeneity  of  interactions  (Figure  25).  Thus,  in  the  deterministic  mode,  a  specific  “shooting”  of 
uniformly  susceptible  targets  by  free  radicals  in  a  homogenous  phase  gave  a  number  of  interactions 
(effective  hits)  proportional  to  the  number  of  free  radicals  (Figure  25  A). 
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Free  Radical  Interactions  with  Cellular  Targets  Interactions  of  Free  Radical  Arrays  with  Cellular  Targets 


Figure  25.  Representations  of  interaction  between:  A.  steady  state  concentrations  of  free  radicals 
(FRad)  and  homogenous  cellular  targets.  B.  arrays  (DOS)  of  the  steady-state  concen¬ 
trations  of  free  radicals  and  heterogenous  cellular  targets. 

On  the  other  hand,  in  the  stochastic  mode,  a  random  “shooting”  of  differently  susceptible 
targets  by  free  radicals  in  non-homogenous  phases,  gave  a  bell-shape  distribution  curve  of  the 
number  of  interactions  (effective  hits)  over  the  arrays  of  free  radicals  (Figure  25  B).  Assuming  that 
this  is  a  normal  distribution,  an  algorithm  for  a  1  minus  cumulative  Gaussian  distribution  function 
was  used  in  the  stochastic  module  to  predict  the  fraction  of  remaining  active  cellular  targets,  after 
Vroegop  et  al.  (1995).  These  assumptions  made  it  possible  to  distinguish  between  a  specific 
inhibition  of  enzymatic  or  transporting  cellular  activities  (e.g.,  protein  tyrosine  phosphatase  activity 
by  H2C>2-derived  free  radicals;  Figure  19  C)  and  a  non-specific  inhibition  (e.g.,  mitochondrial 
reduction  of  MTT  activity  by  H202-derived  free  radicals;  Figure  20  D),  based  on  the  shape  of  a 
dose-response  curve  fitted  to  the  experimental  data. 

OXIDATIVE  STRESS  AND  LIPID  PEROXIDATION 

Pro-oxidant  chemicals  are  bioactivated  by  CYP  (mainly  in  the  liver)  to  metabolites  producing 
free  radicals  and  initiating  lipid  peroxidation  (Larson  and  Bull,  1992;  Cojocel  et  al.,  1989;  Rosen 
and  Rauckman,  1982;  Poyer  et  al.,  1978;  Burdino  et  al.,  1973).  Our  EPR  spin-trapping  study  of 
mouse  liver  treated  in  vitro  with  pro-oxidant  chemicals  confirmed  the  increased  production  of  free 
radicals  (Steel-Goodwin  et  al.,  1996),  but  an  attempt  to  use  the  EPR  signal  as  a  quantitative  end 
point  to  calibrate  the  production  of  free  radicals  from  TCE  turned  to  be  unsuccessful  (Figure  8).  On 
the  other  hand,  TBARS  production  in  the  liver  was  found  to  be  a  convenient  quantitative  end  point, 
following  the  free  radical  generation  (Figure  10)  and  lipid  peroxidation,  and  thus,  it  was  used  as  a 
biomarker  for  PBPD  model  calibrations  in  vitro  (Figures  9, 1 1,  and  14). 


39 


Riely  et  al.,  (1974)  have  demonstrated  that  chlorinated  hydrocarbons,  such  as  CCI4,  increase 
ethane  exhalation  by  mice  due  to  lipid  peroxidation  in  vivo.  Ethane  exhalation  by  animals  treated 
with  CCI4  was  enhanced  by  pretreatment  with  phenobarbital  and  suppressed  by  a-tocopherol 
(vitamin  E),  and  it  was  suggested  that  the  formation  of  carbon  trichloromethyl  free  radical 
(CC13*)  by  CYP  in  the  liver  was  responsible  for  initiation  of  lipid  peroxidation  (Riely  et  al., 
1974).  Bromotrichloromethane,  another  pro-oxidant  chemical  from  which  CYP  generates  CC13*, 
increased  ethane  exhalation  in  rats,  especially  in  vitamin  E-  and  selenium-deficient  group  (Burk 
and  Lane,  1979).  However,  the  same  authors  postulated  also  that  CCI4-  and  BrCCl3 -induced 
lipid  peroxidation  does  not  necessarily  correlate  with  liver  necrosis,  and  may  be  involved  at  an 
early  stage  of  hepatotoxicity  (Burk  and  Lane,  1979).  Since  then,  the  ethane  exhalation  test  was 
proved  to  be  a  reliable,  non-invasive  index  of  oxidative  stress  in  experimental  animals  and 
humans  (Jeejeebhoy,  1991;  Kneepkens  et  al.,  1994),  and  thus,  we  used  it  as  a  biomarker  for 
PBPD  model  calibrations  in  vivo  (Figures  12,  13, 16,  17,  and  18). 

On  the  other  hand,  despite  causing  lipid  peroxidation  as  a  relatively  early  event,  oxidative 
stress  may  actually  stimulate  cellular  proliferation,  induce  apoptosis,  and  at  a  very  high  dose  of 
pro-oxidant  chemical,  it  may  cause  necrosis  (Byczkowski  and  Kulkami,  1996).  Some  of  these 
biological  effects  seemed  to  be  mediated  by  latching  into  the  cellular  signal  transduction  process 
(Byczkowski  and  Channel,  1996).  Because  chemically  induced  oxidative  stress  depends  on  a 
dose  of  pro-oxidant  chemical,  it  seemed  essential  in  our  study  to  follow  effects  of  a  range  of 
doses  of  pro-oxidant  chemicals  on  lipid  peroxidation  (Figures  13,  16,  and  17)  rather  than  to 
check  only  a  time-course  at  one  effective  dose,  as  some  other  authors  have  done  (Riely  et  al., 
1974;  Burk  and  Lane,  1979). 

Comparison  of  Pro-Oxidant  Chemical  Concentrations  in  vitro  vs.  in  vivo : 

Since  the  original  mathematical  model  for  lipid  peroxidation  was  calibrated  with  BrCCl3  in 
vitro  (Tappel  et  al.,  1989),  we  used  the  same  compound  to  compare  a  feasibility  of  measuring  the 
biomarkers  in  our  biological  systems  both  in  vitro  (Figure  1 1)  and  in  vivo  (Figure  12),  with  a 
similar  range  of  the  pro-oxidant  chemical  dosage.  Assuming  that  about  99%  of  the  1  g/kg  of 
BrCCl3  injected  i.p.,  at  the  highest  dose  level,  will  eventually  pass  through  the  liver  in  vivo,  the 
total  cumulated  internal  dose  (AUC)  was  estimated  as  about  2  mg/0.1  g  liver  over  4  h  (or  10 
(j.mol/0.1  g  liver/4  h),  equivalent  to  about  2.5  pmol/0.1  g  liver/h.  This  amount  would  be 
compatible  with  the  initial  concentration  of  0.145  to  2.175  pmol/O.l  g  liver  used  in  vitro  (Figure 
11).  However,  even  though  B1CCI3  partitions  preferentially  in  the  blood,  its  peak  concentration 
in  the  liver  venous  blood  after  the  highest  dose  used  reached  only  about  0.2  mg/L  or  lpM 
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concentration.  Considering  the  liver/blood  partitioning  of  B1CCI3  (0.7),  this  could  produce,  at 
most,  the  peak  concentration  of  70  pmol  per  0.1  g  liver  in  vivo.  Actual  PBPK  sub-model 
simulations  (not  shown  here)  returned  the  highest  sustained  local  concentrations  in  mouse  of  65 
pmol/0.1  g  liver.  However,  since  the  PBPK  sub-model  was  not  validated  with  experimental  data 
for  liver  concentrations  in  vitro  vs.  in  vivo,  these  estimates  remain  tentative.  Assuming  that  the 
free  radical  damage  to  a  biological  system  is  additive,  the  comparison  of  AUC  for  internal  dose 
in  vivo  with  the  initial  pro-oxidant  concentration  in  vitro  seems  to  be  justified  for  the  assessment 
of  biological  effects  (such  as  lipid  peroxidation). 

Effects  of  Pro-Oxidant  Chemicals  on  Lipid  Peroxidation 

From  our  in  vitro  experiments,  presented  in  Figure  1 1,  it  appeared  that  lipid  peroxidation  is 
both  time-  and  dose-dependent  phenomenon;  TBARS  generation  increased  non-linearly  with 
both  increasing  concentration  of  BrCCl3  and  increasing  time  of  incubation.  However,  liver  slices 
incubated  with  the  highest  B1CCI3  concentration  (1.5  mM)  for  1  h  generated  slightly  less 
TBARS  than  with  the  lower,  1  mM  concentration  of  BrCCl3  (Figure  11).  This  decrease  is 
unlikely  to  be  caused  by  the  necrotic  action  of  a  high  dose  of  BrCCl3  acting  for  a  prolonged 
time,  since  the  liver  slices  were  screened  for  signs  of  necrosis  after  the  incubation.  More 
probable  is  a  mechanism  of  CYP  suicide-inhibition  caused  by  the  accumulated  damage  to  the 
enzyme  by  CC13*  and/or  Br*  free  radicals. 

Similarly,  BrCCl3  administered  in  vivo  showed  non-linear  time-  and  dose-dependent 
characteristics  of  lipid  peroxidation  stimulation,  measured  by  ethane  exhalation.  This  appeared 
from  the  data  presented  in  Figures  12  and  13  that  the  maximum  no-effect  dose  of  B1CCI3  of 
about  0.05  g/kg,  when  measured  half  an  hour  after  the  exposure,  dropped  to  almost  0.025  g/kg 
after  one  hour  from  the  exposure.  On  the  other  hand,  the  maximum  stimulation  of  ethane 
exhalation,  measured  1  h  after  the  exposure,  was  reached  at  0.1  g  of  BrCCl3/kg,  while  ten  times 
as  high  a  dose  of  BrCCl3  was  still  not  saturable  when  measured  0.5  h  after  the  exposure  (Figure 
12).  Even  more  dramatic  differences  were  noted  2  h  after  the  exposure  (Figure  13);  at  lg 
BrCC^/kg,  the  ethane  exhalation  was  still  at  the  plateau  level  (0.33  p.p.m.),  while  ten  times 
lower  a  dose  (0.1  g  BrCCl3/kg)  caused  twice  as  high  stimulation  of  ethane  exhalation  (0.68 
p.p.m.). 

These  results  suggest  that  the  extent  of  chemically  induced  oxidative  stress  and  effects  on 
lipid  peroxidation  are  linked  with  both  time  and  dose  of  BrCCl3  by  a  non-linear  function.  The 
non-linearity  may  result  from  limited  supply  of  antioxidants  (threshold)  and  from  accumulated 
over  time  free  radical  insult  to  the  CYP  enzymatic  system  (suicidal  inhibition).  Therefore, 
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considering  biological  effects  and  modeling  in  silico  the  dose-response  for  pro-oxidant 
chemicals,  both  independent  variables  (time  and  dose)  were  taken  into  account  to  assure  realistic 
predictions  and  the  same  should  be  done  in  future  risk  characterizations. 

Similar  characterizations  of  the  dose-response  for  TCE  (Figures  10  and  16)  led  us  to  the 
estimate  of  an  effective  dose-response  range  of  the  local  liver  concentrations  between  0.5  mM 
and  50  mM  (or  0.05  to  5.0  pmol/0.1  g  liver),  and  a  maximum  “no  observable  effect”  dose  of 
TCE  in  vivo  (for  up  to  1  h)  above  0.26  g/kg.  These  estimates  suggested  that  TCE  is  about  an 
order  of  magnitude  less  potent  a  pro-oxidant  than  B1CCI3  (based  on  a  mg  of  mass  comparison). 
On  the  other  hand,  CCI4  had  a  pro-oxidant  potency  similar  to  that  of  BrCCl3  (Figures  14  and 
17). 

Parameters  of  lipid  peroxidation  and  the  feasibility  of  square  root  algorithm,  which  links 
concentration  of  pro-oxidant  with  production  of  free  radicals,  were  determined  in  our  experiments 
with  mouse  liver  slices  treated  in  vitro  with  TCE,  CCI4,  and  B1CCI3.  The  calibrated  algorithms  and 
the  parameters  estimated  in  vitro  were  included  in  the  PBPD  model  simulations  of  time-  and  dose- 
dependent  effects  of  different  doses  of  pro-oxidant  chemicals  on  ethane  exhalation  in  B6C3F1  mice 
in  vivo  (Figures  16  and  18  show  examples  of  TCE  and  CCI4  effects  vs  model  predictions).  It  seems 
that  the  PBPD  model  described  kinetics  and  dynamics  of  chemically  induced  lipid  peroxidation 
relatively  well,  at  least  for  TCE  and  CCI4. 

EFFECTS  OF  FREE  RADICALS  ON  CELLULAR  TARGETS 

Cellular  effects  of  oxidative  stress  and  free  radical  insult  to  cellular  targets  are  well  described 
in  the  literature  (Byczkowski  and  Gessner,  1988;  Kulkami  and  Byczkowski,  1994a;  1994b; 
Byczkowski  and  Channel,  1996;  Byczkowski  and  Kulkami,  1996).  Depending  on  the  biochemical 
mechanism,  the  effects  may  include  inhibition  of  enzymatic  activities,  inhibition  of  cellular 
transporter  activities,  uncoupling  of  oxidative  phosphorylation  and  inhibition  of  the  respiratory 
chain,  interference  with  cellular  membrane  receptors  and  signal  transduction  pathway, 
deregulation  of  transcription  and  gene  expression,  and  damage  to  vital  macromolecules  (lipids, 
proteins,  and  DNA).  Cellular  mechanisms  and  problems  associated  with  risk  characterization  of 
pro-oxidant  chemicals  were  reviewed  in  detail  in  the  previous  report  by  Byczkowski  and 
Flemming  (1995).  Essentially,  the  interaction  of  pro-oxidant  derived  free  radicals  with  cellular 
targets  may  be  either  specific  (uniformly  susceptible  targets,  suspended  in  the  same  phase  as  free 
radicals)  or  random  (differently  susceptible  targets,  or  suspended  in  non-homogenous  phases). 
The  specific  interactions  were  adequately  described  by  deterministic  module  of  our  BBDR  sub- 
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model  (Figures  19  and  21),  whereas  random  interactions  were  relatively  well  simulated  by  the 
stochastic  module  (Figures  20  and  21). 

Setting  a  Hypothesis  with  the  PBPD  Model 

There  is  growing  body  of  evidence  that  pro-oxidant  chemicals  and  oxidative  stress  can 
interfere  with  cellular  signal  transduction  pathway  (Byczkowski  and  Channel,  1996).  Chen  and 
Chan  (1993)  using  3T3-L1  cells,  cultured  in  a  serum-free  medium,  demonstrated  that  pro-oxidant 
chemicals  (e.g.,  orthovanadate  or  DMNQ)  increased  [3H]thymidine  incorporation  to  DNA  and 
enhanced  expression  of  the  c -fos  gene.  The  authors  suggested  that  pro-oxidant  compounds, 
causing  oxidative  stress,  increased  tyrosine  protein  phosphorylation  early  in  the  signal 
transduction  cascade,  leading  to  augmentation  of  cell  proliferation  (Chen  et  al.,  1990;  Chen  and 
Chan,  1993).  This  effect  may  be  accounted  for  by  a  selective  effect  on  redox-sensitive  protein 
tyrosine  phosphatase  (PTyrPase).  Typically,  PTyrPases  contain  a  nucleophilic  cysteinyl  residue  in 
their  catalytic  center  (Stone  and  Dixon,  1994)  and  their  enzymatic  activities  are  rapidly  inhibited 
by  small  disulfides  (Ziegler,  1985).  It  seems  that  the  cysteinyl  residue  must  be  kept  in  the  reduced 
-SH  form,  therefore,  thiol-directed  reagents  that  oxidize  it  cause  inhibition  of  PTyrPases  (Fischer 
et  al.,  1991). 

Accordingly,  Mendelson  et  al.,  (1996)  demonstrated  that  oxidative  stress  caused  by  CCI4  in 
the  liver,  activated  the  JNK  family  of  protein  kinases  and  increased  the  AP-1  transcription  factor 
binding  to  DNA.  Because  mitogen  activated  kinases  (MAPK)  must  be  kept  phosphorylated  to 
transmit  this  effect,  it  is  possible,  that  inhibition  of  PTyrPase  may  stop  deactivation  of  phosphatase 
MKP-1  which  has  high  affinity  for  phosphorylated  p38  as  a  substrate.  This,  in  turn,  may  cause  a 
rapid  dephosphorylation  of  p38,  which  was  down-regulated  following  the  oxidative  stress 
(Mendelson  etal.,  1996). 

Since  a  similar  subversion  of  signal  transduction  was  suggested  for  TCE  (Maronpot  et  al., 
1995)  but  was  never  proved,  it  is  possible  that  TCE-derived  free  radicals  can  affect  intracellular 
redox  potential  and  may  lead  to  the  inhibition  of  PTyrPase  activity  by  causing  oxidation  of 
essential  cysteinyl  thiols.  Our  EPR  spin  trapping  experiments  confirmed  the  presence  of  free 
radicals  in  TCE-treated  liver,  but  failed  to  reveal  their  identity  (Steel-Goodwin  et  al.,  1995;  1996). 
Assuming  that  the  TCE-derived  free  radicals  are  well  water-soluble  (Gonthier  and  Barret,  1989; 
Mason,  1992;  Ni  et  al.,  1994),  they  should  interact  with  PTyrPase  suspended  in  the  same  phase  as 
free  radicals,  in  a  specific  way  (the  cysteinyl  -SH  residues  represent  uniformly  susceptible 
targets).  If  this  is  the  case,  within  the  relevant  TCE  local  concentrations  in  the  liver,  an 
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exponential  inhibition  of  PTyrPase  should  be  observed  with  a  50%  inhibitory  concentration 
around  0.5  xnM  (0.05  pmol  TCE/0.1  g  liver),  as  predicted  by  the  BBDR  sub-model.  However,  if 
only  lipid-soluble  secondary  and  tertiary,  lipid  peroxide-derived  free  radicals  are  available  in  the 
liver,  they  should  result  in  a  non-specific  random  inhibition  of  PTyrPase,  giving  a  sigmoidal,  dose- 
response  with  a  50%  inhibitory  concentration  around  22.0  mM  (2.2  pmol  TCE/0.1  g  liver;  Figure 
26). 


BBDR:  Simulated  Effects  of  TCE  on  PTyrPase  Activity 


TCE  (yumo  I  f  0.1  9  liver) 


Figure  26.  The  results  of  BBDR  sub-model  simulations  of  the  expected  inhibition  of  protein 
tyrosine  phosphatase  activity  in  the  liver  of  mice  in  vivo  after  treatment  with  the 
estimated,  effective  pro-oxidant  range  of  TCE  doses  (resulting  in  the  local  concen¬ 
trations  between  Cmin  =  0.05  and  Cmax  =  5.0  pmol  TCE/0.1  g  liver;  other  parameters 
are  listed  in  Table  4). 

TCE  -  local  concentration  of  TCE  in  the  liver  [pmol/0.1  g  liver].  A.  ENH  -  remaining  percentage  of 
uninhibited  activity  by  the  deterministic  BBDR  module  (x  102);  I50A  -  50%  inhibitory 
concentration  (0.05  pmol  TCE/0.1  g  liver).  B.  PROB  -  remaining  percentage  of  uninhibited  activity 
by  the  stochastic  BBDR  module  (x  102);  I50B  -  50%  inhibitory  concentration  (2.2  pmol  TCE/0.1  g 
liver). 

Deterministic  and  stochastic  PBPD  model  theoretical  simulations  (unconfirmed 
experimentally)  of  the  remaining  percentage  of  the  enzymatic  activity  of  protein  tyrosine 
phosphorylase  in  the  liver  inhibited  by  TCE-derived  free  radicals  specifically  (Figure  26  A)  and 
nonspecifically  (Figure  26  B)  showed  quite  different  dose-dependent  characteristics  and  different 
by  two  orders  of  magnitude  I50  values,  within  the  relevant  range  of  TCE  concentrations.  Which 
dose-response  profile  will  actually  occur  remains  to  be  seen,  but  it  was  possible  with  the  aim  of  the 
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PBPD  model  to  simulate  the  predicted  outcome  within  the  relevant,  internal  pro-oxidant  dosage  of 
TCE,  and  set  a  verifiable  working  hypothesis  for  the  future  experimental  research. 

CONCLUSION 

The  resultant  PBPD  hybrid  model  may  be  used  for  pharmacodynamic  description  of 
chemically-induced  oxidative  stress  in  mice,  for  planning  the  future  experiments  and  setting  the 
verifiable  working  hypotheses,  and  potentially  may  be  useful  for  a  risk  characterization. 
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APPENDIX:  Source  Codes  of  PBPD  Model  Written  in  ACSL® 


•  PBPD  Sub-Model 

Sub-Model  for  Pro-Oxidant  Induced  Lipid  Peroxidation  *.CSL  FILE 

’This  PBPD  sub-model  predicts  ethane  exhalation  by  mice  in  a  closed  * 
'chamber  after  dosing  with  CC14,  TCE  and  BrCC13  based  on  TBARS' 

'model  published  by  Byczkowski  et  al.,  1996;  interlinked  with  ' 
'PBPK/TCPK  sub-model  Das  et  al.,  1994  and  Seckel  and  Byczkowski' 

'1996.  Final  version  1/03/97  ' 

PROGRAM:  PBPD  FOR  P-OX-INDUCED  OXIDATIVE  STRESS 


INITIAL 

LOGICAL  CC  $ ' Flag  set  to  .TRUE,  for  closed  chamber  runs 


'  Miscellaneous  commands 

CONSTANT  CC  =  .true.  $’ Default  to  close  chamber 

CONSTANT  PBPK  =  1.  $’  Default  to  interlink  with  PBPK  model 

'  to  use  two-compartment  CPK  switch  PBPK=0. 

CONSTANT  CPK=  0.  $’ Switch  CPK=1  to  use  two- compartment  PK 

’PEROXIDATION  PARAMETERS  FOR  LIVER  [per  0.1  g  of  liver] 


CONSTANT 

PUF 

7.0 

CONSTANT 

PXZLUF= 

12.0 

CONSTANT 

PXZHUF= 

24.0 

CONSTANT 

AN  0X1  = 

0.0037 

CONSTANT 

ACT1  = 

0.0003 

CONSTANT 

GSH 

0.6 

CONSTANT 

PXTTBA= 

0.1 

CONSTANT 

LPUF  = 

3. 

CONSTANT 

HPUF  = 

4.  * 

CONSTANT 

ANOX2  = 

0. 

CONSTANT 

ANOX3  = 

0. 

CONSTANT 

EFANOl 

= 

1. 

CONSTANT 

EFAN02 

= 

0. 

CONSTANT 

EFAN03 

= 

0. 

CONSTANT 

ANOXUF 

= 

2. 

CONSTANT 

ACT  2 

= 

0. 

CONSTANT 

ACACT1 

= 

1. 

CONSTANT 

ACACT2 

- 

0. 

CONSTANT 

ACTDGF 

= 

1.75 

CONSTANT 

IND2 

= 

0. 

CONSTANT 

PTIND1 

= 

4408. 

CONSTANT 

PTIND2 

= 

0. 

CONSTANT 

INDLF 

= 

0.0001 

CONSTANT 

PXRATE 

= 

0.00029 

CONSTANT 

AUTOXF 

= 

0.00013 

CONSTANT 

GPENZA 

= 

1. 

CONSTANT 

PXREDF 

= 

0.17 

CONSTANT 

PHYSPX 

= 

0. 

CONSTANT 

ANREG 

= 

0.0007 

CONSTANT 

ANOXR 

- 

0.001 

CONSTANT 

BCKGD 

= 

0.0 

CONSTANT 

ACR 

= 

0.025 

$ ’  [micromol/0 . 1  g  liver] 

$ ’peroxidizability  of  PUFA  L[l/h] 

$ ’peroxidizability  of  PUFA  H[l/h] 

$ ' Vit . E  antiox [mcmol/0 . Ig  liv] 

$ ' activatorl [mcmol  cytP450/0 . Igliv] 

$’ glutathione  [mcmol/0. 1  g  liver] 

$’ yield  of  TBARS/mol  hydroperoxide 
$' LA-derivative  PUFA [mcmol /0 . 1  g] 

$ ’HA-derivative  PUFA [mcmol/0 . 1  g] 
$’non-Vit.E  ant i ox [mcmol/0 . 1  g  liv] 

$’ added  antioxidant [mcmol/0 . lg  liv] 

$’ effectiveness  of  Vit.E 
$ 1  effectiveness  of  non-Vit.E  a-o 
$’ effectiveness  of  added  a-o 
$’ antioxidant  use  factor  [1/mcmol] 

$’ activator  2  [mcmol/0. 1  g  liver] 

$’ activity  of  activator  1 
$’ activity  of  activator  2 
$’ activator  degradation  fctr  [1/mcmol] 

$ ’ inducer2 [mcM  chemical/0 . lg  liv] 

$* potency  of  inducer  1  [1/mcmol] 

$' potency  of  inducer  2  [1/mcmol] 

$' inducer  loss  factor  [1/h] 

$ ’peroxidation  rate  [mcmol/mcmol] 

$  1  autooxidation  factor  [1/h] 

$ ’ glutathione  peroxidase [1/mcmol] 

$ ’hydroperoxide  reduction  factor  [1/h] 

$ ’physiological  levels  of  Hperoxides 
$’ antioxidant  regenerated  in  situ 
$* antioxidant  regeneration  rate  constant 
$ ’ TBARS  in  control  [mcm/0.1  g] 

$' activator  degradation  rate  constant 


’TIMING  COMMANDS 


CONSTANT  POINTS  =  200. 

CINT=TSTOP/POINTS 

CONSTANT  TCHNG  =4.0  $' Length  of  inhalation  exposure  to  P-Ox  [h] 

CONSTANT  TINF  =  0.01  $'Length  of  IV  infusion  with  P-Ox  [h] 

'repeated  gavage  dosing 

constant  DAYS=0. 08333  $' Duration  of  simulation  [DAYS]  if  pdays=0. 
constant  pdays=0.  $'No-gavage  days  in  cycle 

'PARAMETERS  FOR  MOUSE  IN  VIVO 

CONSTANT  KLC  =  .050  $' First  order  chamber  loss  [Lin  fraction/h] 

CONSTANT  KS  =  100000.  $' Suppression  rate  constant  for  metabolism 

CONSTANT  AVBW  =  38.0  $' Average  BW  of  mice=Tot  W/n  [g] 

BW  =  AVBW/1000  $' Average  body  weight  per  mouse  [kg] 

CONSTANT  QPC  =  30.00  $' Alveolar  ventilation  rate  [1/h] 

CONSTANT  QCC  =16.5  $' Cardiac  output  [1/h] 

QLC  =  0.24  -  QGC  • 

CONSTANT  QGC  =  .175  $' Fractional  blood  flow  to  gut 

CONSTANT  QFC  =  .05  $' Fractional  blood  flow  to  fat 

CONSTANT  QSC  =  .238  $'  Fractional  blood  flow  to  slow 

CONSTANT  QRC  =  .472  $' Fractional  blood  flow  to  rapid 

CONSTANT  VLC  =  .05  $' Fraction  liver  tissue 

CONSTANT  VGC  =  .033  $' Fraction  gut  tissue 

CONSTANT  VSC  =  .558  $' Fraction  slow  tissue 

CONSTANT  VRC  =  .031  $' Fraction  rapid  tissue 

CONSTANT  .  VFC  =  .1  $' Fraction  fat  tissue 

'  PARAMETERS  FOR  ETHANE 

CONSTANT  PLA  =  0.828  $' Liver/air  partition  coefficient 

CONSTANT  PGA  =  0.996  $' Gut/air  partition  coefficient 

CONSTANT  PFA  =  2.444  $'Fat/air  partition  coefficient 

CONSTANT  PSA  =  0.979  $' Slowly  perfused  tissue/air  partition 

CONSTANT  PRA  =  0.996  $' Richly  perfused  tissue/air  partition 

CONSTANT  PEB  =  1.305  $' Blood/air  partition  coefficient 

PEL=PLA/PEB  $ 'Liver/blood  partition  coefficient 
PEG=PGA/PEB  $ 'Gut/blood  partition  coefficient 
PEF=PFA/PEB  $ ' Fat/blood  partition  coefficient 
PES=PSA/PEB  $' Slow/blood  partition  coefficient 
PER=PRA/PEB  $' Rich/blood  partition  coefficient 

CONSTANT  BACKE=0 .  $' Background  ethane  concentration  [ppm] 

CONSTANT  MEW  =  30.  $ 'Molecular  weight  Et  [g/mol] 

I 

CONSTANT  VMEXC=0.286  $' Maximum  velocity  of  metabolism  [mg/h-lkg] 

CONSTANT  KEM  =  0.51  $ 'Michaelis-Menten  constant  [mg/L] 

CONSTANT  KEFC=  2.786  $' First  order  metabolism  rate  constant [1/h-l kg] 

CONSTANT  NRATS=  5.  $' Number  of  mice  [for  closed  chamber] 

CONSTANT  VCHC  =  0.745  $ 'Volume  of  closed  chamber  [L] 

CONSTANT  SODA  =0.005  $' Volume  of  soda  lime  [L] 

'  PARAMETERS  FOR  P-Ox 

CONSTANT  BAB  =  0.99  $' Fraction  absorbed  from  i.p.  to  portal  blood 

CONSTANT  FAB  =  0.0001  $'Fraction  absorbed  from  i.p.  directly  to  fat 
'  Note:  (BAB  +  FAB)  <  1. 

'Carbon  Tetrachloride  -  specific  constants  ' 

CONSTANT  PLCC  =  3.14  $' Liver /blood  partition  coefficient 


CONSTANT  PFCC  =79.4  $'Fat/blood  partition  coefficient 
CONSTANT  PRCC  =  3.14  $' Richly  perfused  tissue/blood  partition 

CONSTANT  PBCC  =  4.52  $' Blood/air  partition  coefficient 

CONSTANT  PSCC  =  2.43  $' Slowly  perfused  tissue/blood  partition 

CONSTANT  MWCC  =153.82  $ 'Molecular  weight  P-Ox  [g/mol] 

CONSTANT  DOSE  =0.  $'Dose  of  active  P-Ox  [mg/kg] 

'Estimate  of  initial  concentration  of  chemical  in  the  liver  [mg/L] 
CLO  =  DOSEIP/VLC 


CONSTANT 

KAIP  =1.45 

$’ First 

.  order  i.p.  uptake  [1/h] 

CONSTANT 

DOSEIP  =  1500. 

$'i.p. 

dose  [mg/kg] 

DIP  =  DOSEIP*BW 

$'i.p. 

dose  per  animal  [mg] 

CONSTANT 

PD0SEC=  0. 

,  $ 'Oral 

dose  [mg/kg] 

DOSEC  =  PD0SEC*BW 

$'p.o. 

dose  per  animal  [mg] 

CONSTANT 

1 

KAC  =  1. 

$ 'Oral 

uptake  rate  [1/hr] 

CONSTANT 

IVDOSC  =  0. 

$ ' i . v. 

dose  [mg/kg] 

CONSTANT 

CONCC  =0. 

$' Inhaled  concentration  [ppm] 

'Metabolic  constants  for  CC14  ' 

CONSTANT  VMAXCC=0 . 65  $ 'Maximum  velocity  of  metabolism  [mg/h-lkg] 
CONSTANT  KMC  =0.25  $ 'Michaelis-Menten  constant  CC14  [mg/L] 

CONSTANT  KFCC=  0.0  $'First  order  metabolism  rate  constant  [1/h-lkg] 

'Pharmacokinetic  transfer  constants  fitted 

CONSTANT  klO  =  0.3  $' Pharmacokinetic  transfer  micro-constant  [1/h] 

CONSTANT  alpha=  1.5  $' Pharmacokinetic  macro-constant  [1/h] 

CONSTANT  beta=  1.6  $' Pharmacokinetic  macro-constant  [1/h] 

CONSTANT  EFFE  =0.001  $' Efficiency  of  ethane  generation  [molar  ratio] 

'  Inhalation  P-Ox  Exposure  definition 

IF  (CC)  RATS  =  NRATS  $' Closed  chamber  simulation 

IF  (CC)  KL  =  KLC 

IF  (.NOT.CC)  RATS  =  0.  $'Open  chamber  simulation 

IF  (.NOT.CC)  KL  =  0. 

'Turn  off  chamber  losses  so  concentration  remains  constant 

AI0C  =  CONCC*VCH*MWCC/24450.  $' Initial  amount  of  P-Ox  in  chamber  [mg] 

'  PARAMETERS  FOR  CLOSED  CHAMBER 

CONE=  BACKE 

VCH  =  VCHC- (RATS*BW) -SODA  $'Net  chamber  volume  [L] 

AEI0  =  CONE*VCH*MEW/24450.  $' Initial  amount  of  E  in  chamber  [mg] 

'  INITIALIZATION  ' 

'RESETS  INITIAL  CONDITIONS  BEFORE  PEROXIDATION  ' 

CONSTANT  ACTLOS=0.  $' activator  loss  ' 

CONSTANT  AUTOX  =0.  $' autooxidation  ' 

CONSTANT  PXREDG=0 .  $' Hydroperoxides  red.  by  GSH  peroxidase' 

CONSTANT  PXLUF  =0.00199  $ ' L-Hydroperoxides  formed  ' 

CONSTANT  PXHUF  =0.0053  $ ' H-Hydroperoxides  formed  ' 

CONSTANT  PXREM  =0.  $' accumulated  remaining  Hydroperoxides  ' 

CONSTANT  TPX  =0.  $' accumulated  total  Hydroperoxides  formed' 

CONSTANT  ILR  =0.  $’rate  of  inducer  loss  =0.  ' 

CONSTANT  CL  =  0.  $'no  inducer 

'Reset  output  arrays  at  initialization  ' 

ETH  =0. 

CONC=  0. 
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Scaled  parameters  for  mouse 


QC  =  QCC*BW**0 . 74 
QP  =  QPC*BW**0 . 74 
QL  =  QLC*QC 
QG  =  QGC*QC 
QF  =  QFC*QC 
QS  =  0 . 24*QC-QF 
QR  =  0 . 7  6*QC-QL-QG 
QLB=  QL+QG 
VL  =  VLC*BW 
VG  =  VGC*BW 
VF  =  VFC*BW 
VS  =  0 . 82*BW-VF 
VR  =  0. 101*BW-VL-VG 
VMEX  =  VMEXC*  BW*  *0.7 


KEF  =  KEFC/BW**0 . 3 
VEK  =  VMEXC /KEM 
VMAXC4  =  VMAXCC*BW**0 . 7 
KFC4  =  KFCC/BW**0 . 3 


INTEGER  DAY 

'  Repeated  gavage  dosing  with  P-Ox 

tstop=  (days+PDAYS) *24 . 

CINT=  tstop/points 

DAY=- 1 .  $ ' TO  START  GAVAGE  ON  MONDAY  -1,  TUES  0,  WEDN  1,  ETC 

END  $'End  of  Initial 

DYNAMIC 

»  GAV  =  FEED  MICE  p.o.  with  P-Ox  YES=1,  N0=0 . 

DISCRETE  CAT1 

INTERVAL  CAT  =24.  $' EXECUTE  CAT1  EVERY  24  hr 

DAY=DAY+1 

IF (DAY. GT. DAYS)  GOTO  OUT 
IF (MOD (DAY, 7) .GE.5)  GOTO  OUT 

GAV  =  1.  $' GAVAGE  =  YES 

SCHEDULE  CAT 2  .AT.  T+0.01  $' SCHEDULE  END  OF  GAVAGE 

OUT..  CONTINUE 
END  $ ' END  OF  CAT1 

DISCRETE  CAT 2 

GAV  =  0. 

END  $ ' END  OF  CAT 2 

ALGORITHM  IALG  =  2 
DERIVATIVE 


<<<<MODULE  LOCAL  DOSE  OF  P-Ox»» 

Estimate  of  actual  cone,  of  chemical  in  the  liver  CL(C) [mg/L]  by  two- 
compartment  open-system  classic  pharmacodynamic  CPK  sub-model  as  CL. 
Alternatively,  CL(C)  is  calculated  continuously  by  PBPK  sub-model  as 
CLC' 


««SUB-MODEL  CPK»» 


$’ GAVAGE  =  NO 
$'Gear  stiff  method 


CL  =  CL0*kl0* (exp (-beta*t) -exp (-alpha*t))/ (alpha  -  beta) 


*  ««SUB-MODEL  PBPK  FOR  P-Ox»» 

'Includes  code  for  suppression  of  metabolism 

IVRC  =  IVDOSC*BW/TINF 

'  i.p.  dosing  with  P-Ox 

RMRIP  =  -KAIP*MRIP  $'Rate  of  change  of  amount  in  i.p.  cavity  [mg/h] 
MRIP  =  INTEG (RMRIP, DIP)  $ 'Amount  of  toxicant  in  i.p.  cavity  [mg] 

'  AIP  =  amount  of  P-Ox  absorbed  from  i.p. 

RAIP  =  KAIP*MRIP  $'Rate  absorption  from  i .p.  cavity ■ [mg/h] 

AIP  =  DIP  -  MRIP  $ 'Amount  of  toxicant  absorbed  from  i.p.  [mg] 

'  CIC  =  Concentration  in  inhaled  air  [mg/L] 

CIZONE  =  RSW( (T.LT.TCHNG) .OR.CC, 1. ,  0. ) 

RAIC  =  RATS*QP* (CAC/PBCC-CIC)  -  (KLC*AIC) 

AIC  =  INTEG (RAIC, AI0C) 

CIC  =  AIC/VCH*CIZONE 
■  CP  =  CIC*24450./MWCC 

'  Repeated  gavage  dosing  with  P-Ox 

'  MRC  =  Amount  remaining  in  stomach  [mg] 

RMRC  =  gav*dosec/tinf  -  raoc 
MRC  =  integ (rmrc, 0 . ) 

'  AOC  =  total  mass  input  from  stomach 

RAOC  =  kac*mrc 
AOC  =  integ (raoc, 0 . ) 

'  AGC  =  Amount  in  gut  [mg] 

RAGC  =  QG* (CAC-CVGC)  $ ' This  is  in  addition  to  RAOC 
AGC  =  INTEG (RAGC, 0. ) 

CVGC  =  AGC/ (VG*PRCC) 

CGC  =  AGC/VG 

'  CAC  =  Concentration  in  arterial  blood  [mg/L] 

CAC  =  (QC*CVC+QP*CIC) / (QC+ (QP/PBCC) ) 

AUCBC  =  INTEG (CAC, 0. ) 

'  AXC  =  Amount  exhaled  [mg] 

CXC  =  CAC/PBCC 

CXPPMC  =  (0.7*CXC+0.3*CIC) *24450.  /MWCC 
RAXC  =  QP*CXC 
AXC  =  INTEG (RAXC,  0.  ) 

'  ASC  =  Amount  in  slowly  perfused  tissues  [mg] 

RASC  =  QS* (CAC-CVSC) 

ASC  =  INTEG (RASC, 0. ) 

CVSC  =  ASC/ (VS*PSCC) 

CSC  =  ASC/VS 

'  ARC  =  Amount  in  rapidly  perfused  tissues  [mg] 

RARC  =  QR* (CAC-CVRC) 

ARC  =  INTEG ( RARC, 0. ) 

CVRC  =  ARC/ (VR*PRCC) 

CRC  =  ARC/VR 

'  AFC  =  Amount  in  fat  tissue  [mg] 
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RAFC  =  QF* (CAC-CVFC)  +  (FAB*RAIP) 

AFC  =  INTEG (RAFC, 0 . ) 

CVFC  =  AFC/ (VF*PFCC) 

CFC  =  AFC/VF 

ALC  =  Amount  in  liver  tissue  [mg] 

RALC  =  QL* (CAC-CVLC) +QG* (CVGC-CVLC) +raoc-RAMC+ (BAB*RAIP) 
ALC  =  INTEG (RALC, 0. ) 

CVLC  =  ALC/ (VL*PLCC) 

CLC  =  ALC/VL 

AMC  =  Amount  metabolized  [mg] 

RAMC  =  (VMAXC4*CVLC) / (KMC+CVLC)  +  KFC4*CVLC*VL  - 
AMC  =  INTEG (RAMC,  0.) 
ambc=  arnc/bw 


'  IVC  =  Intravenous  infusion  rate  [mg/hr] 

IVC  =  IVRC* ( 1 . -step (tinf ) ) 

'  CVC  =  Mixed  venous  blood  concentration  P-Ox  [mg/L] 

CVC= (QF*CVFC+ (QL+QG) *CVLC+QS*CVSC+QR*CVRC+IVC+RAIP* (1- (BAB+FAB) ) ) /QC 

'  CTMASS  =  mass  balance  of  P-Ox  [mg] 

CTMASS  =  AFC+ALC+ASC+ARC+AMC+AXC+AGC 

cbal  =  aoc-ctmass  $’gavage,  repeated;  mass  bal' 

cball=  mrc-ctmass 

'  CDOSEX  =  Net  amount  P-Ox  absorbed  [mg] 

CDOSEX  =  AIC+AOC+IVC*TINF+AIP 
BWCC  =  CDOSEX/bw 

MOLCC=  BWcc/MWCC  $'Milimoles  P-Ox  absorbed  [mmoles/kg] 


iaa/^end  OF  PBPK  SUB-MODEL  FOR 


P-Ox 


AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA? 


- link  - 

Link  with  local  concentration  of  P-Ox:  ' 

Switch  interlinking  with  PBPK  or  CPK  yields  normalized  concentration 
in  liver,  compatible  with  slices  [mcromole/0 . 1  g] 

IND1  =  (CL*CPK  +  CLC*PBPK) / (MWCC*10)  $'MWCC  *  10  =  1538.2 


A=CL/(MWCC*10) 

B=CLC/ (MWCC*10) 

AUCA  =  INTEG  (A,  0) 

AUCB  =  INTEG (B,0) 

PROCE  DURAL 

if  (IND1.LE.0.)  IND1=0.0  $ 1  Prevents  from  attempt  SQRT  negative  value1 

END  $  1  End  of  procedural 


<<<<MODULE  LIPID  PEROXIDATION»» 


PROCEDURAL 

IF  (IND1.GE.CL0)  IND1  =  CL0/(MWCC  *  10) 
IF  (LPUFRE.LE.O)  LPUFRE  =  0. 

IF  (HPUFRE.LE. 0)  HPUFRE  =  0. 

IF  ( ACTLOS . GE . ACTEF)  ACTLOS  =  ACTEF 
IF  (ANOXRE . LE . 1 . e-10 )  ANOXRE  =  l.e-10 
IF  (GSHREM.LE. 0)  GSHREM  =  0. 

IF  ( INDLOS . GE . INDEF)  INDLOS  =  INDEF 
IF  (PXREDA.GE.TPX)  PXREDA  =  TPX 
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END  $'End  of  procedural  ' 

i****************BBPB  SUB-MODEL  FOR  LIPID  PEROXIDATION**********' 

f _ _  _ _ _ 1 

*#s  Corresponding  to  equations  in  Byczkowski  et  al.  (1996)  1 

'This  part  has  been  calibrated  in  vitro  in  mouse  liver  slices  * 

I _ —  — _ T 

'1.  Remaining  polyunsaturated  fatty  acids  [micromol/0 . lgliv]  ' 

LPUFRE  =  LPUF  -  PXLUFA  -  AUTOXA/2 
HPUFRE  =  HPUF  -  PXHUFA  -  AUTOXA/2 

'2.  Effective  activator  [micromoe/0 . 1  g  liver]  ' 

ACTEF  =  ACT1*ACACT1  +  ACT2*ACACT2 

'3.  Activator  loss  [micromole/0.1  g  liver]  ' 

ACTLOS  =  ACTEF*ACTDGF*TPX 

'4.  Remaining  activator  [micromole/0.1  g  liver]  ' 

ACTREM  =  (ACTEF  -  ACTLOS) * exp (-ACR*INDREM*t) 

'5.  Effective  inducer  [micromole/0.1  g  liver]  ' 

INDEF  =  SQRT (IND1*PTIND1)  +  IND2*PTIND2 

'6.  Remaining  inducer  [micromole/0. lg  liv]  ' 

INDREM  =  INDEF  -  INDLOS 

'6. a.  Inducer  loss  rate  [micromol/0 . 1  g  liver/hr]  ' 

ILR  =  INDEF* INDLF 

'7.  Activated  inducer  [micromol/0 . 1  g  liver]  ' 

ACTIND  =  INDREM*ACTREM 

'8.  Effective  antioxidant  [micromol/0 . lg  liver]  ' 

ANOXEF  =  AN0X1*EFAN01  +  AN0X2*EFAN02  +  AN0X3*EFAN03 

'9.  Remaining  antioxidant  [mcromol/0.1  g]  ' 

ANOXRE  =  ANOXEF  -  ANOXEF*TPX*ANOXUF  +  ANREG*exp (-ANOXR*t) 

' 10. Hydroperoxides  formed  by  action  of  activated  inducer  on 
'  PUFA  [micromol/0 . 1  g  liver/hr]  ' 

PXLUF  =  LPUFRE*PXZLUF*ACTIND*PXRATE/ANOXRE 
PXHUF  =  HPUFRE*PXZHUF*ACTIND*PXRATE/ANOXRE 

' 12 .Autooxidation  [micromole/0.1  g  liver/hr]  ' 

AUTOX  =  (LPUFRE  +  HPUFRE) *AUTOXF*TPX/ANOXRE 

' 14 .Accumulated  total  hydroperoxides  formed  [mcmol/0.1  g  liver]' 
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TPX  =  AUTOXA  +  PXLUFA  +  PXHUFA  +  PHYSPX 

'15. Remaining  glutathione  [micromol/0 . 1  g  liver]  1 

GSHREM  =  GSH  -  PXREDA 

’ 16. Hydroperoxides  reduced  by  glutathione  peroxidase  1 

1  [micromol/0 . 1  g  liver/hr]  ' 

PXREDG  -  PXREM* GPEN ZA* GSHREM* PXREDF 

' 18 .Accumulated  remaining  hydroperoxides  [micromole/0.1  g  liv.]T 


PXREM  —  TPX  -  PXREDA 

*19. Amount  of  TBARS  from  accumulated  remaining  hydroperoxides  T 
1  [micromol/0 . 1  g  liver]  1 

TBARS  =  PXREM*  PXTTBA  +  BCKGD 

’  Inducer  lost  over  time  [micromol/0 . 1  g  liver]  * 


INDLOS  =  INTEG { ILR, 0 . ) 

1 11  .Accumulated  hydroperoxides  formed  by  action  of  activated 
1  inducer  on  PUFA  [micromole/0.1  g  liver] 

PXLUFA  =  INTEG (PXLUF,  0.) 

PXHUFA  =  INTEG ( PXHUF ,  0.) 

1 13. Accumulated  autooxidation  [micromol/0 . 1  g  liver] 


AUTOXA  =  INTEG  (AUTOX,  0.) 

’ 17  .Accumulated  hydroperoxides  reduced  by  glutathione 
f  [micromole/0.1  g  liver] 


PXREDA  =  INTEG ( PXREDG,  0.) 

AAAAAEnd  of  BBPD  Sub-model  for  lipid  peroxidationAAAAAAAAAAAAAA ' 


’ - - - LINK - T 

1  Link  with  Accumulated  Remaining  Hydroperoxides  produced  in  the  liver  1 

1  Assumed  lipid  peroxidation  in  the  liver  only,  evoked  by  P~Ox 
concentration  * 

1  Lipid  peroxidation  is  driven  by  hydroperoxides  generated  by  free  1 

’radicals  depending  on  SQRT  of  local  concentration  of  P-Ox  in  liver 

ETH  =  PXREM* EFFE  $T Molar  Amount  of  ethane  produced  [mcM/0.1  g  liver]1 

’Rate  of  evolution  of  ethane  produced  in  the  liver  [mg  E/hr /mouse]  ? 


REOX= (ETH*VL*300 ) / (t+le-12) 

’REOX  =  f (SQRT(RPOX) ) ’ 

’  300=MEW/1000/0 . 0001  over  time  prevented  to  start  from  dividing  by  0 


<<<<MODULE  ETHANE  EXHALATION»» 
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'This  part  has  been  calibrated  in  mice  in  vivo  ' 

'****PBPD  SUB-MODEL  FOR  ETHANE***************************************** 
'  CEI  =  Concentration  in  inhaled  air  (mg/L) 

REAI  =  RATS*QP* (CEA/PEB-CEI ) - (KL*AEI) 

AEI  =  INTEG (REAI, AEIO)  $  'Chamber 

CEI  =  AEI/VCH  $  'With  N  mice 

CEP  =  CEI*24450 . /MEW  $' concentration  in  closed  chamber 


'  CEA  =  Concentration  in  arterial  blood  [mg/L]  ' 

CEA  =  (QC*CEV+QP*CEI ) / (QC+ (QP/PEB) ) 

'  AEX  =  Amount  exhaled  per  mouse  [mg]  ' 

CEX  =  CEA/PEB 

CXEPM  =  (0 . 7*CEX+0 . 3*CEI ) *24450 . /MEW 
REAX  =  QP*CEX 
AEX  =  INTEG (REAX, 0 . ) 

'  AES  =  Amount  in  slowly  perfused  tissues  per  mouse  [mg]  ' 

REAS  =  QS* (CEA-CEVS) 

AES  =  INTEG (REAS,  0.  ) 

CEVS  =  AES/ (VS*PES) 

CES  =  AES/VS 

'  AER  =  Amount  in  rapidly  perfused  tissues  per  mouse  [mg]  ' 

REAR  =  QR* (CEA-CEVR) 

AER  =  INTEG (REAR,  0.  ) 

CEVR  =  AER/ (VR*PER) 

CER  =  AER/VR 

'  AEF  =  Amount  in  fat  tissue  per  mouse  [mg]  ' 

REAF  =  QF* (CEA-CEVF) 

AEF  =  INTEG (REAF,  0 . ) 

CEVF  =  AEF/ (VF*PEF) 

CEF  =  AEF/VF 

'  AEG  =  Amount  in  gut  tissue  per  mouse  [mg]  ' 

REAG  =  QG* (CEA-CEVG) 

AEG  =  INTEG (REAG,  0.  ) 

CEVG  =  AEG/ (VG*PEG) 

CEG  =  AEG/VG 

'  AEL  =  Amount  in  liver  tissue  per  mouse  [mg]  ' 

REAL  =  QL*  (CEA-CEVL)  +QG*  (CEVG-CEVL)  -REAM+REOX 
AEL  =  INTEG (REAL,  0 . ) 

CEVL  =  AEL/ (VL* PEL) 

CEL  =  AEL/VL 

'  AEM  =  Amount  metabolized  per  mouse  w/suppression  (KS)  [mg] 

REAM  =  (VMEX*CEVL) / (KEM+CEVL* (1+CEVL/KS) )  +  KEF*CEVL*VL  $ ' [mg/h]  ' 

AEM  =  INTEG (REAM, 0 . )  $' Amount  [mg]  ' 

'  CEV  =  Mixed  venous  blood  concentration  per  mouse  [mg/L]  1 

CEV  =  (QF*CEVF+  (QL+QG) *CEVL+  QS*CEVS+  QR*CEVR) /QC 

'AMOUNT  INHALED  PER  MOUSE 
REINH  =  QP*CEI 
AEINH  =  INTEG (REINH, 0) 
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'TMESS  =  MASS  BALANCE  PER  MOOSE  ’ 

TMESS  =  ( AE  S  + AER+ AE  F-KAEM+AEL+ AEX+ AEG ) 

BEL  =  AEINH  -  TMESS 

i  aaaaEND  0F  SOB-MODEL  FOR  ETHANEAAAAAAAAAAAAA AAAAAAAAAAAAAAAAAAAAAAAAAA ’ 


TERMT (T.GE.TSTOP) 


$ 'Termination  at  TSTOP 


END  $'End  of  derivative  ' 

END  $'End  of  dynamic  ' 

I  AAAAAAAAA/\AAAAAAAAA/\/\AAAAAAAAAAAAAAA/\AAAAAAAA.y\AAAAAAA/\AAAAAAAA/\l 


TERMINAL  $' Resets  parameters  to  initial  values1 

ACTREM  -  ACTEF 
INDREM  -  INDEF 
ANOXRE  =  ANOXEF 

'Save  arrays  of  dependent  variable* 

ETH  =  CEP 
CONC  -  Indl 

END  $  1  End  of  terminal  1 


AAAAEND  OF  PBPD  SUB-MODEL  FOR  P-OX-INDUCED  LIPID  PEROXIDATIONAAA 1 


END  $ ' End  of  program  ' 


•  Command  files  for  CCI4  and  TCE 

PBPD  Sub-Model  *.CMD  FILE  for  CC14 

'ETHANE.CMD' 

'GAS  UPTAKE/EXHALATION  DATA  FOR  ETHANE  IN  MICE' 
PREPAR  T,  'ALL' 

SET  GRDCPL=.F.  $' Turns  off  grid  lines' 

PROCED  CONDI T 

'TABLE  3:  CONDITIONS  FOR  MICE' 

SET  KLC=0 . 05,  KS=100000 . , POINTS=300 . 

SET  PLA=0 . 8285,  PGA=0.996,  PFA=2.444,  PRA=0.996 
SET  PSA=0 .979,  PEB=1.305,  EFFE=0.001 
SET  MEW=30 . ,  AVBW=31. 

SET  VMEXC=0 .286,  KEM=0.51,  KEFC=2.786 
SET  NRATS=5,  VCHC=0.745,  SODA=0 . 005 
SET  QPC=30. 0,  QCC=16 . 5 

DISPLAY  QPC,  QCC,  VMEXC, KEM, KEFC, PEB, PLA,  PGA,  PFA,  PSA 
END 

PROCED  FIG15 

'Ethane  Uptake  from  Closed  Chamber  1  ppm' 

SET  TITLE  =  'ETHANE  UPTAKE' 

SET  BACKE=1 .227,  DOSEIP=0. 

SET  KLC= .05,  KS=100000 . , 

SET  PLA=0 .8285,  PGA=0.996,  PFA=2.444,  PRA=0.996 

SET  PSA=0. 979,  PEB=1.305 

SET  MEW=30 . ,  DAYS=0 . 125 

SET  VMEXC=0 .286,  KEM=0.51,  KEFC=2.786 

SET  NRATS=4 ,  VCHC=0.745,  SODA=0 . 005 

SET  QPC=30 . 0 ,  QCC=16 .5,  AVBW=31. 


DATA 

T 

0.0 


CEP 


INITIAL 
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0.08300 

0.988875 

0.17000 

0.966695 

0.25000 

0.961003 

0.33000 

0.903836 

0.50000 

0.822943 

0.75000 

0.743472 

1.00000 

0.946910 

1.25000 

0.599756 

1.50000 

0.564610 

1.75000 

0.530173 

2.00000 

0.387864 

2.25000 

0.443700 

2.50000 

0.384679 

2.75000 

0.348181 

3.00000 

0.311571 

END  $ 'END 
START 

OF  DATA’ 

PLOT  CEP 
END 

PROCED  CONDCCL 

'  TABLE  1&2:  CHEMICAL-SPECIFIC  CONDITIONS  FOR  LIPID  PEROXIDATION  BY  CC14 T 
SET  PTIND1=44 08 . ,  ACR=0.025, 

SET  kl0=0 . 3, alpha=l . 5,beta=l . 6, KAIP=1 .45 
SET  ACTDGF=1 .  75 ,  INDLF=0 . 0001, MWCC-153 . 8 
SET  ANREG=0 .0007,  ANOXR=0.001 
END 

PROCED  FIG17A 

1  Ethane  Exhalation  after  CC14  0.075  g/kg' 

SET  TITLE  =  *  ETHANE  EXHALATION  AFTER  CC14 1 

PREPAR  t,  *  ALL  * , CL 

SET  DAYS=0. 083333,  NRWITG- . F. 

SET  BACKE=0 . , POINTS=300, AVBW=32 . 82 
SET  DOSEIP-75 . 


1  [h] 

DATA 

T 

0. 

0.0833 

0.166667 

0.25 

0.333333 

0.416667 

0.5 

0.75 

1. 

1.25 

1.5 

1.75 

2. 

END  $  T  END  OF 
START 


[ppm]  1 

CEP 

0. 

0.027758 

0.081802 

0.057029 

0.030957 

0.019462 

0.04063 

0.046445 

0.035008 

0.097736 

0.120978 

0.121889 

0.184791 

DATA1 


PLOT  CEP,  'lo'-O. ,  1  hi f  =0.25,  ,‘xhi,=2.2 
END  $  *  End  of  file* 


PROCED  FIG17B 

1  Ethane  Exhalation  after  CC14  0.15  g/kgl 
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SET  TITLE  =  'ETHANE  EXHALATION  AFTER  CC14 * 

PREPAR  t, 'ALL', CL 

SET  DAYS=0. 083333,  NRWITG=. F. , 

SET  BACKE=0 . , POINTS=200, AVBW=32 .56 
SET  DOSEIP=150. 


'  [h] 

DATA 

T 

0. 

0.0833 

0.166667 

0.25 

0.333333 

0.416667 

0.5 

0.75 

1. 

1.25 

1.5 

1.75 

2. 


[ppm]  ' 

CEP 

0. 

0. 

0. 

0.010312 

0.043063 

0.071015 

0.020494 

0.107916 

0.089671 

0.076225 

0.099269 

0.076916 

0.042414 


END  $ ' END  OF  DATA' 
START 


PLOT  CEP, 'lo'=0. , 'hi '=0.25, 'xhi'=2.2 
END  $ ' End  of  file' 


PROCED  FIG17C 

'Ethane  Exhalation  after  CC14  0.3  g/kg' 

SET  TITLE  =  'ETHANE  EXHALATION  AFTER  CC14  ' 

PREPAR  t, 'ALL', CL 

SET  DAYS=0. 083333,  NRWITG=. F . , 

SET  BACKE=0 . , POINTS=200, AVBW=38  . 

SET  DOSEIP=300 . 


1  [h] 

DATA 

[ppm]  1 

T 

CEP 

0. 

0. 

0.0833 

0. 

0 . 166667 

0.038776 

0.25 

0.072132 

0.333333 

0.070974 

0.416667 

0.123789 

0.5 

0.096211 

0.75 

0.109947 

1. 

0.097474 

1.25 

0.092724 

1.5 

0.073632 

1.75 

0.110104 

2. 

0.095158 

END  $ ' END  OF  DATA' 

START 

PLOT  CEP, 'lo'=0. , 'hi '=0 . 25, 'xhi'=2.2 
END  $ ' End  of  file' 

PROCED  FIG17D 

'Ethane  Exhalation  after  CC14  1.5  g/kg' 

SET  TITLE  =  'ETHANE  EXHALATION  AFTER  CC14 ’ 
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PREPAR  t,  * ALL 1 , CL 

SET  DAYS=0 .083333,  NRWITG=.F., 

SET  BACKE=0 . , POINTS-300, AVBW=31 . 8 
SET  DOSEIP=1500. 


1  [h] 

DATA 

T 

0. 

0.0833 

0.166667 

0.25 

0.333333 

0.416667 

0.5 

0.75 

1. 

1.25 

1.5 

1.75 

2. 


[ppm]  f 

CEP 

0. 

0. 

0. 

0.065464 

0.090988 

0.027442 

0.059182 

0.170162 

0.194608 

0.148473 

0.175496 

0.132515 

0.182836 


END  $  '  END  OF  DATA 
START 


PLOT  CEP, ' lo '=0 . , ’hi 1 =0 .25, 'xhi'=2.2 
END  $  '  End  of  file1 


PROCED  FIG18 

'Ethane  Exhalation  after  CC14  1.5,  0.3,  0.15,  0.075  g/kg' 
SET  TITLE  =  'ETHANE  EXHALATION  AFTER  CC14 ' 

PREPAR  /clear  t,CEP 

SET  DAYS=0. 083333,  NRWITG= . T . , FTSPLT= . T . 

SET  BACKE=0 . ,  POINTS-300. 

SET  AVBW=35. 

SET  NRATS=5 , 

SET  DOSEIP-75. 

START 


SET  DOSEIP=150 . 
START 


SET  DOSEIP=300 . 
START 

SET  DOSEIP—1500 . 


1  [h] 

DATA 

T 

0. 

0.0833 

0.166667 

0.25 

0.333333 

0.416667 

0.5 

0.75 

1. 

1.25 

1.5 

1.75 


[ppm]  1 

CEP 

0. 

0.  . 

0. 

0.065464 

0.090988 

0.027442 

0.059182 

0.170162 

0.194608 

0.148473 

0.175496 

0.132515 


INITIAL 
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INITIAL 


2. 

0. 

0.0833 

0.166667 

0.25 

0.333333 

0.416667 

0.5 

0.75 

1. 

1.25 

1.5 

1.75 

2. 

0. 

0.0833 
0 . 166667 
0.25 

0.333333 

0.416667 

0.5 

0.75 

1. 

1.25 

1.5 

1.75 

2. 

0. 

0.0833 

0.166667 

0.25 

0.333333 

0.416667 

0.5 

0.75 

1. 

1.25 

1.5 

1.75 

2. 

END  $  ?  END 
START 


0.182836 

0. 

0. 

0.038776 

0.072132 

0.070974 

0.123789 

0.096211 

0.109947 

0.097474 

0.092724 

0.073632 

0.110104 

0.095158 

0. 

0. 

0. 

0.010312 

0.043063 

0.071015 

0.020494 

0.107916 

0.089671 

0.076225 

0.099269 

0.076916 

0.042414 

0. 

0.027758 
0.081802 
0.057029 
0.030957 
0.019462 
0.04063 
0.046445 
0.035008 
0.097736 
0.120978 
0.121889 
0.184791 
F  DATA1 


INITIAL 


INITIAL 


PLOT  CEP,  1 lo'=0. ,  *  hi  '“0 . 25,  ,xhif==2.2 
END  $'End  of  filef 


PBPD  Sub-Model  *.CMD  FILE  for  TCE 

PROCED  FIG16 

SET  TITLE='MICE  ETHANE  AFTER  TCE : 2 . 6, 1 , 0 . 26g/kg ' 

SET  TSTOP=2 . ,  NRWITG=.t. , 

PREPAR  t,  'all' 

SET  kl0=0 . 03,  alpha=0 . 01,  beta=2.0,  EFFE=0.001, 

SET  PBCC=13 . 4 ,  PLCC=2 .03,  PFCC=4 1 . 3,  PRCC=2 . 03,  PSCC=1. 
SET  MWCC=131 . 5,  VMAXCC=33.,  KMC=0.25,  KFCC=2 . 4 
SET  ACTDGF=0. 0014,  PTIND1=250,  KAIP=1.,  DOSE=2600., 
START 
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SET  DOSE=1000 . 
START 


SET  DOSE=260. 

1  th]  [ppm]  1 
DATA 

T  CEP 

0.  0.000001  INITIAL 

0.25  0.036148 

0.333  0.027971 

0.75  0.038132 

1.  0.042805 

1.25  0.045230 

1.5  0.093413 

1.75  0.075427 

2.  0.110495 

0.  0.000001  INITIAL 

0.1667  0.032907 

0.25  0.016053 

0.3333  0.026683 

1.5  0.026005 

1.75  0.056478 

2.  0.060730 

0.  0.000001  INITIAL 

0.75  0.034561 

END  $,end  of  data1 

START 

PLOT  CEP,  CLC,  'xhi' =2 . 

END  $Tend  of  file* 

***★★★***★**★**★********★**■*■■*•*-*  'k'k'k'k'k'k'k'k'k'k'k'k-k-k-k'k-k'tr-k-k-k-k'k-k'k'k'k'k'k'k-k  +  'k-te'k-k'k'k'k-k'k 

•  BBDR  Sub-Model 

Deterministic  and  stochastic  modules  *.CSL  FILE 

PROGRAM:  FREE  RADICAL  DOSE  RESPONSE 

’A  program  that  calculates  amount  of  free  radicals  from  initial  ’ 

’local  concentration  of  P-Ox  and  estimates  their  cellular  effect  1 
’based  on  Vroegop  et  al.  1995.  Final  version  for  SIMUSOLV  1/10/97’ 

==  OPEN (UNIT=41, STATUS= ’ NEW* , FILE= ’ GRAPHG . FIG 1 } 

WRITE (41, 10) 

10..  FORMAT  (’  PROC  NDATA  ’  /  ’  DATA’/  ’  PCONC  PROB  ’  ) 

INITIAL  $’ Pre-execution  section  of  program  ’ 

CONSTANT  ki  =0.001  $’rate  constant  of  free  radical  formation  1 

CONSTANT  kt  =0.002  $’rate  constant  of  free  radical  recombination’ 

’  ki  kt  ’ 

’C  - >  FR  +  FR  - >  Nonradical  products  ’ 

’assuming  early  phase  C=C0,  and  the  steady  state  condition:  ’ 

’ dFR/dt  =  ki*C0  -  kt*FR*FR  =  0.  (const,  ki,  kt  [l/uM*h]  )  ’ 

CONSTANT  kd  =  0.075  $’rate  constant  of  receptor  inactivation  ’ 

’  kd 

’FR  +  RECEPTOR  >  RECEPTOR  Inactive 

’assuming  first  order  process  for  receptor  inactivation  and  ’ 

’uniform  sensitivity  fo  FR  with  rate  constant  kd  [100%/uM*h]  ’ 

CONSTANT  10  =  1.  $’ Initial  concentration  of  receptors  [%*E-2] ’ 
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CONSTANT  tp  =  1.  $ ' Time  of  receptors  exposure  to  FR  [h] 

CONSTANT  a  =  0.075  $' Receptor  population  response  constant 

CONSTANT  TSTOP  =  2.  $' Length  of  experiment  [h] 

CONSTANT  TME=  1.  $'Time  to  maximum  effect  [h] 

CONSTANT  POINTS=  1.  $' Number  of  communication  intervals 

Ft  =  tp/TME  $' Fraction  of  max  time  [ratio] 
'Parameters  for  dose-response  simulation 

CONSTANT  CMIN  =  100.  $' Starting  concentration  [uM] 

CONSTANT  CMAX  =  1000.  $' Final  concentration  [uM] 

CONSTANT  CDELT  =  100.  $' Concentration  interval  [uM] 

'Miscellaneous  parameters 
INTEGER  I , ND 

CINT  =  TSTOP/POINTS  $' Interval  for  saving  data  to  * . RRR  file 
REAL  DOS (1000) ,  PROB  (1000) , DOSC (1000) 

'arrays  of  data  saved  for  plotting  PROBability  vs  DOSe 
'Initialize  variables  for  dose-response  calculation 
PCONC=0 .  $  FRad=0 .  $  Inh=1.0 

CO  =  CMIN-CDELT 

ND  =  INT ( (CMAX  -  CMIN) /CDELT)  +  1 
1  =  1 

'Start  of  dose-response  loop 

RESTRT. .C0=C0+cdelt 

END  $'End  of  INITIAL  section 


DYNAMIC  $' Beginning  of  execution  section  of  program 


FR  =  SQRT (ki*C0/kt) 


$' Concentration  of  free  radicals  [uM] 


,AAA^AAAAAAAAAAAAAAAAADETERMINISTIC  MODULE! a aaaaaaaaaaaaaaaaaaaaaa 

In  =  10  *  exp (-kd*FR*tp)  $ 'Amount  of  active  receptors  remaining 

'after  exposure  to  FR  for  time  =  tp 

TERMT(T.GE. TSTOP)  $’Stop  simulation  when  T  >=  TSTOP 


END  $'End  of  dynamic  section 

TERMINAL  $' Post-execution  section  of  program 


PCONC  =  CO  $'Save  current  concentration  [uM] 

FRad  =  FR  $'Save  cone,  of  free  radicals  [uM]  ' 

Inh  =  In  $ 'Relative  amount  of  active  receptors  ' 

DOS (I)  =  FRad  $'Save  I-th  FR  concentration  to  array  ' 

DOSC ( I ) =PCONC  $'Save  I-th  chemical  cone,  to  array  ' 

1  =  1+1 

CALL  LOGD ( .FALSE. ) 

IF  (C0.lt . cmax)  goto  restrt  $' Restart  to  initial  unless  done  ' 

,aaaaaaaaaaaaaaaaaaaaaaaaST0CHASTIC  MODULEaa aaaaaaaaaaaaaaaaaaaaaa 

x ************** ******START  FORTRAN  SUBROUTINE********************' 


PROCEDURAL (PROB  =  ND, DOS, Ft) 

CALL  OMPHI (ND, DOS, PROB, Ft) 

END  $'End  of  procedural' 
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PROCEDURAL (DOSC, PROB) 

DO  20  K  =  1,  ND 

WRITE (41, 30)  DOSC (K) , PROB (K) 

20. . CONTINUE 

30..  FORMAT (  F9 . 3, 2X, F6 . 4 ) 

WRITE (41, 40) 

40..  FORMAT (’■  END') 

END 

END  $'End  of  terminal  section' 

END  $'End  of  program' 


SUBROUTINE  OMPHI (ND, Dl, PROB, Ft) 

C 

C  OMPHI  FINDS  1  -  CUMMULATIVE  GAUSSIAN  DISTRIBUTION. 

C 


INTEGER  I,ND,NDMAX 
PARAMETER ( NDMAX=1 0  0  0 ) 

REAL  M,  SD,  D,  SI,  S2,N,  Ft,  PROB  (NDMAX)  ,  Dl  (NDMAX) 


c 

VARIABLE 

TYPE 

DESCRIPTION 

c 

I 

INTEGER 

INDEX  VARIABLE 

c 

ND 

INTEGER 

NUMBER  OF  DOSES 

c 

NDMAX 

INTEGER 

MAX  NUMBER  OF  DOSES 

c 

Ft 

REAL 

FRACTION  OF  TIME  TO  MAX  EFFECT 

c 

M 

REAL 

MEAN  OF  DOSE 

c 

SD 

REAL 

STANDARD  DEVIATION  OF 

DOSE 

c 

D 

REAL 

DOSE  VALUE 

c 

Dl 

REAL 

ARRAY  OF  DOSE  VALUES 

c 

PROB 

REAL 

ARRAY  OF  RESULTS  OF  1 

-  CUM.  GAUSS 

C 

C 

C 


DIST 


CALCULATES  THE  MEAN  AND  STANDARD  DEVIATION  OF  DOSE. 

51  =  0.0 

52  =  0.0 

DO  10  I  =  1,  ND 
D  =  D1(I) 

51  =  SI  +  D 

52  =  S2  +  D*D 


10  CONTINUE 

N  =  FLOAT (ND) 

M  =  Sl/N 

SD  =  SQRT ( (S2  -  S1*S1/N) / (N  -  1)) 

C 

C  FINDS  AND  OUTPUTS  THE  PROBABILITY  OF  CUMMULATIVE  GAUSSIAN. 

C 


DO  20  I  =  1, ND 
D  =  Dl (I) 

PROB ( I )  =  1.0  -  Ft  *  PHI (M, SD, D) 
20  CONTINUE 
RETURN 
END 

FUNCTION  PHI (M, SD, D) 


C 

C  PHI  FINDS  THE  CUMMULATIVE  GAUSSIAN  DISTRIBUTION  USING  THE 

C  ERROR  FUNCTION  AND  THE  COMPLEMENTARY  ERROR  FUNCTION. 

C 

REAL  M, SD, D, X 
C 

C  VARIABLE  TYPE  DESCRIPTION 

C  M  REAL  INPUT  MEAN  OF  DOSE(D) 

C  SD  REAL  INPUT  STANDARD  DEVIATION  OF  DOSE(D) 
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D  REAL  INPUT  DOSE 

X  REAL  INTERNAL  VARIABLE 

X  —  (D  -  M)/SD 
X  =  X/ SQRT (2.0) 

IF  (X  .GE.  0.0)  THEN 

PHI  =  (1.0  +  ERF (X)0  /2 . 0 
ELSE 

PHI  =  ERFC(-X) /2.0 
ENDIF 
RETURN 
END 

FUNCTION  erf (x) 

ERF  FINDS  THE  ERROR  FUNCTION. 

USES  gammp 

REAL  erf, x, gammp 
if (x. It . 0 . ) then 

erf=-gammp { .5, x**2) 
else 

erf=gainmp  ( .  5,x**2) 
endif 
return 
END 

FUNCTION  gammp (a, x) 

GAMMP  FINDS  THE  INCOMPLETE  GAMMA  FUNCTION  P(a,x). 

REAL  a,  gammp,  x 

USES  gcf , gser 

REAL  gammcf ,  gamser ,  gin 

if (x.lt .0. .or.a.le.0. ) pause  'bad  arguments  in  gammp ' 
if (x.lt. a+1 . ) then 

call  gser (gamser, a, x, gin) 
gammp=gamser 
else 

call  gcf (gammcf , a, x, gin) 
gammp=l . -gammcf 
endif 
return 
END 

SUBROUTINE  gcf ( gammcf ,  a, x, gin) 

GCF  RETURNS  THE  INCOMPLETE  GAMMA  FUNCTION  Q(a,x)  EVALUATED 
BY  ITS  CONTINUED  FRACTION  REPRESENTATION  AS  GAMMCF.  IT 
ALSO  RETURNS  InGAMMA  ( a )  AS  GLN  t 

INTEGER  ITMAX 

REAL  a,  gammcf , gin, x, EPS, FPMIN 

PARAMETER  ( ITMAX=100 , EPS=3 . e-7 , FPMIN=1 . e-30 ) 

USES  gammln 

INTEGER  i 

REAL  an, b, c, d, del, h, gammln 

gln=gammln  (a) 

b=x+l.-a 
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c=l . / FPMIN 
d=l . /b 
h=d 

do  11  i— 1, ITMAX 
an=-i* (i-a) 
b=b+2 . 
d=an*d+b 

if (abs(d) . It . FPMIN) d=FPMIN 
c=b+an/c 

if (abs(c) . It. FPMIN) c=FPMIN 

d=l./d 

del=d*c 

h-h*del 

if (abs (del-1.) .It. EPS) goto  1 
11  continue 

pause  'a  too  large,  ITMAX  too  small  in  gcf’ 

1  gammcf=exp (-x+a*log (x)  -gin)  *h 

return 
END 

SUBROUTINE  gser  (gamser , a, x, gin) 

GSER  RETURNS  THE  INCOMPLETE  GAMMA  FUNCTION  P(a,x)  EVALUATED 
BY  ITS  SERIES  REPRESENTATION  AS  GAMSER.  IT  ALSO  RETURNS 
InGAMMA  ( a )  AS  GLN. 

INTEGER  ITMAX, n 

REAL  a,  gamser,  gin, x,  EPS,  ap,  del,  sum,  gammln 
PARAMETER  ( ITMAX=100,  EPS=3 .  e-7 ) 

USES  gammln 

gln=gammln ( a ) 
if (x. le. 0 . ) then 

if (x. It. 0. ) pause  'x  <  0  in  gser' 
gamser=0. 
return 
endif 
ap=a 

sum=l . /a 
del=sum 

do  11  n=l, ITMAX 
ap=ap+l. 
del=del*x/ap 
sum=sum+del 

if (abs (del) . It . abs ( sum) *EPS) goto  1 
continue 

pause  'a  too  large,  ITMAX  too  small  in  gser* 
gamser=sum*exp (-x+a*log (x) -gin) 
return 
END 

FUNCTION  gammln (xx) 

GAMMLN  RETURNS  THE  VALUE  In  [GAMMA  (xx)  ]  FOR  xx  >  0 . 

REAL  gammln,  xx 
INTEGER  j 

DOUBLE  PRECISION  ser , stp, tmp, x, y, cof ( 6) 

SAVE  cof, stp 

DATA  cof , stp/7  6 . 1800917294714 6d0, -86. 505 32 032941677 dO, 

*24 . 014 09824 083 091d0, -1 . 231739572450155d0, 

* . 1208  6509738  6617 9d-2, 
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5395239384 953d-5, 2. 506628274 6310005d0/ 
x=xx 
y=x 

tmp=x+5 . 5d0 

tmp= (x+0 . 5d0 ) *log (tmp) -tmp 
ser=l . 000000000190015d0 
do  11  j=l,6 
y=y+l . dO 

ser=ser+cof ( j ) /y 
11  continue 

gammln=tmp+log (stp*ser/x) 

return 

END 

FUNCTION  erf c (x) 

C 

C  ERFC  FINDS  COMPLIMENTARY  ERROR  FUNCTION. 

C 

REAL  erfc,x 
C 

C  USES  gairunp,  gammq 

C 

REAL  gairunp,  gammq 
if  (x.  It .  0 . )  then 

erfc=l.+gammp{  . 5,x**2) 
else 

erf c=gairanq  ( . 5, x**2 ) 
endif 
return 
END 

FUNCTION  gammq (a, x) 

C 

C  GAMMQ  RETURNS  THE  INCOMPLETE  GAMMA  FUNCTION  Q(a,x)  =  1  - 

C  P (a,x) . 

C 

REAL  a,  gammq,  x 
C 

C  USES  gcf,gser 

C 

REAL  gammcf , gamser, gin 

if (x. It . 0 . . or . a . le . 0 . ) pause  'bad  arguments  in  gammq* 
if (x.lt .a+1. ) then 

call  gser (gamser, a, x, gin) 
gammq=l . - gamser 
else 

call  gcf (gammcf , a, x, gin) 
gammq=gammcf 
endif 
return 
END 

**************** ************************************************************ 

•  Command  files  and  graphics  for  deterministic  and  stochastic  modules 

BBDR  Sub-Model  *.CMD  FILES 


i  aaa/va/s/n/vaaa^a/sa/vaa^qATA  FROM  LITERATURE 


AAAAAAAAAAAAAAAAAAAAA  f 


' - Steel-Goodwin  et  al.,  (1995) - 

PROCED  FIG8 

'kinetics  of  FRAD  production  from  TCE  in  precision  cut' 
'mouse  liver  slices'' 

PREPAR  ’CLEAR’ , Pconc,  FRad 
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SET  TITLE  =  ' BBPD  Module:  Effect  of  TCE  on  Generation  of  Free  Radicals' 

SET  ki=900, kt=200, cmin=l . , cdelt=l . , cmax=40 . , tp=0 . 333 

SET  TSTOP=0 . 333 , NRWITG= . T . ,  TME=0 . 333 

'data  normalized  to  physiological  backgrounds . ' 

'uM/g  uM/g' 

DATA 

Pconc  FRad 

0.  0. 

9.85  3.5 

18.9  9.4 

32.8  11.9 

37.15  15. 

END  $'end  of  data' 

START 

PLOT  FRAD, 'TAG'  =  '-  FR  CONCENTR.  (um/g)  ' ,  'XTAG'  =  '  -  TCE  (um/g)  ' 

END  $'end  of  Linda' 

' - Vroegop,  et  al.  (1995) - ' 

PROCED  FIG19A 

'effect  of  cumene  hydroperoxide  on  amino  acid' 

'transport  1  =  19000  cpm' 

SET  TITLE  =  ' BBDR:  Effect  of  Cum.OOH  on  AA  Transporter’ 

PREPAR  ' CLEAR', Pconc, Inh 

SET  cmin=10.,cdelt=10.,cmax=1000.,TP=l.,TSTOP=2. 

SET  NRWITG= . T . , ki=100, kt=200, kd=0 . 075,  TME=1 . 


'  uM 

%  Control' 

DATA 

Pconc 

Inh 

0. 

1. 

1. 

0.84 

3.3 

0.91 

10. 

0.95 

33. 

0.71 

50. 

0.69 

80. 

0.57 

130. 

0.52 

200. 

0.47 

330. 

0.37 

1000. 

0.22 

END 

$  fend 

of 

START 

PLOT 

Inh, 1 lo 1 =0 . 

,  'TAG 

T  _  1 

' XTAG f  -  1 

-  CUM.OOH 

END 

$ '  end 

of 

data' 

-  ACT.  AA  TRANSPORTERS’  ... 
(uM)  ' 

CHDOSEAA' 


PROCED  FIG19B 

'effect  of  H202  on  glucose' 

'transporter  1  =  25250  cpm’ 

SET  TITLE  =  'BBDR:  Effect  of  H202  on  Glucose  Transporter' 
PREPAR  'CLEAR', Pconc, Inh 

SET  cmin=100 . , cdelt=100. , cmax=10000. ,TP=1 . ,TSTOP=2 . 

SET  NRWITG= . T . , ki=18 , kt=200 , kd=0 . 1,  TME=1 . 


'uM 

%  Control 

DATA 

Pconc 

Inh 

0. 

1. 

3.3 

0.92 

10. 

0.95 
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33.  0.89 

100.  0.76 

330.  0.57 

1000.  0.48 

3300.  0.22 

10000.  0.05 

33000.  0.04 

END  $'end  of  data' 

START 

PLOT  Inh, 'lo'=0. ,  ' TAG ' = '  -  ACT.  GLUC  TRANSPORTERS'  .. 

'xhi ’=10000, 'XTAG'  =  '  -  H202  (uM) ' 

END  $ 'end  of  H202DG' 

' - Heffetz  et  al.  (1990) - - 

PROCED  FIG19C 

'effect  of  H202  on  protein  Tyr  phosphatase' 

'Tyr  Pase  activity  measured  as  [32P]  remaining' 

' in  [32P] poly- (Glu, Tyr) ' 

SET  TITLE  =  ' BBDR:  Effect  of  H202  on  PTyrPase' 

PREPAR  ' CLEAR ’ , Pconc, Inh 

SET  cmin=0. 01, cdelt=10 . , cmax=500 . , TP=0 . 3333, TSTOP=0 . 42 
SET  NRWITG=. T . , ki=18 , kt=200,  kd=l . ,  TME=0 . 3333 

'mM  %  Control' 

DATA 

Pconc  Inh 

0.  1. 

10.  0.4953 

20.  0.7143 

25.  0.5333 

30.  0.5143 

50.  0.3143 

100.  0.1430 

200.  0.1143 

300.  0.0476 

500.  0.00952 

END  $'end  of  data’ 

START 

PLOT  Inh,  ' lo ' =0 . , ' TAG '  = '  -  ACT.  PTPase'... 

'XTAG'  =  '  -  H202  (uM) ' 

END  $ ' end  of  H202PTP ' 

' - Hecht  and  Zick  (1992) - 

PROCED  FIG19D 

'effect  of  vanadate  on  protein  Tyr  phosphatase' 

'Tyr  Pase  activity  measured  as  [32P]  remaining' 

'in  [32P] poly- (Glu, Tyr) ’ 

SET  TITLE  =  'BBDR:  Effect  of  Vanad.  on  PTyrPase' 

PREPAR  'CLEAR' , Pconc, Inh 

SET  cmin=0. 1, cdelt=10 . , cmax=1000 . , TP=0 . 5, TSTOP=0 . 633 
SET  NRWITG=. T . , ki=18 , kt=200 , kd=l . ,  TME=0 . 7 


*  uM 

DATA 

%  Control*  . 

Pconc 

Inh 

0. 

1. 

1. 

0.529 

10. 

0.456 

100. 

0.191 

1000. 

0.294 

END 

$  *  end 
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START 

PLOT  Inh,  ' lo ' =0 . , ' xlo ' =0 . , ' TAG'  = '  -  ACT.  PTPase'... 

' XTAG '  =  '  -  VANAD  (uM) ' 

END  $ ' end  of  VANADTP' 

' - Vroegop,  et  al.  (1995) - 

PROCED  FIG20A 

'effect  of  6 -OH  dopamine  on  amino  acid  transport  1  =  5600  cpm' 
PREPAR  'CLEAR', Pconc, Inh 

SET  TITLE  =  ' BBDR:  Effects  of  6-OH  Dopamine  on  AA  Transporter' 
SET  TP=1 . , TSTOP=2 . , ki=200 ,  kd=0 . 075 , cmin=10 . , cdelt=3 . ,  cmax=300 . 
SET  NRWITG= . T . , Kt=200, TME=1 . 

'uM  %  Control' 

DATA 

Pconc  Inh 

0.  1. 

10.  1. 

20.  0.916 

40.  1. 

60.  0.916 

80.  0.8779 

90.  0.8321 

100.  0.4351 

200.  0.2214 

300.  0.2137 

END  $'end  of  data' 

START 

PLOT  Inh, 'lo'=0.,  ' TAG ' = '  -  ACT.  AA  TRANSPORTERS'  ... 

'XTAG'  =  '  -  6-OH  DOPAMINE  (uM) ' 

END  $'end  of  60HDAA ' 

PROCED  FIG20B 

'effect  of  H202  on  amino  acid  transport  1  =  6350  cpm' 

SET  TITLE  =  'BBDR:  Effect  of  H202  on  AA  Transporter' 

PREPAR  ' CLEAR' , Pconc, Inh 

SET  cmin=10 . , cdelt=l . , cmax=100 . , TP=1 . , TSTOP=2 . 

SET  NRWITG=. T . , ki=18, kt=200, kd=0 .075,  TME=1 . 

'uM  %  Control' 

DATA 

Pconc  Inh 

0.  1. 

0.33  1.0 

10.  0.86 

33.  0.56 

55.  0.42 

100.  0.34 

END  $ ' end  of  data' 

START 

PLOT  Inh, 'lo'=0. ,  ' TAG ' = '  -ACT.  AA  TRANSPORTERS'  ... 

'XTAG'  =  '  -  H202  (uM) ' 

END  $'end  of  H202DAA' 

PROCED  FIG20C 

'effect  of  6-OH  dopamine  on  mitochondria  staining  1  =  0.51  OD' 
SET  TITLE  =  'BBDR:  Effect  of  6-OH  Dopamine  on  Mitochondria' 
PREPAR  ' CLEAR' , Pconc, Inh 

SET  kd=0 . 05, cmin=40 . , cdelt=2 . , cmax=200 . , TP=1 . ,  TSTOP=2 . 

SET  NRWITG=.T.,ki=200,kt=200,  TME=2.04 

' uM  %  Control ' 
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DATA 

Pconc 

Inh 

0. 

1. 

10. 

1.0 

20. 

1.0 

40. 

1.0 

60. 

0.982 

80. 

0.927 

90. 

0.918 

100. 

0.795 

200. 

0.498 

300. 

0.502 

500. 

0.466 

END 

START 

PLOT  Inh, 

1 lo ' =0 . 

'XTAG'  =  '  - 

END 


$'end  of  data' 

' TAG ' = '  -  ACT.  MITOCHONDRIA'  ... 
6-OH  DOPAMINE  (uM) ' 

$ ' end  of  60HDMIT ' 


PROCED  FIG20D 

'effect  of  H202  on  mitochondria  staining  1  =  0.49  OD' 
SET  TITLE  =  ' BBDR:  Effect  of  H202  on  Mitochondria' 
PREPAR  ' CLEAR ' , Pconc, Inh 

SET  cmin=10 . ,  cdelt=l . , cmax=100 . , TP=1 . , TSTOP=2 . 

SET  NRWITG=. T . , ki=18 , kt=200 , kd=0 . 05,  TME=1. 


'  uM  %  Control ' 


DATA 

Pconc 

Inh 

0. 

1. 

0.33 

0.97 

10. 

1. 

18. 

0.95 

33. 

0.69 

55. 

0.55 

100. 

0.50 

END 

START 

PLOT  Inh, 

1  lo 1 =0 

'XTAG'  = 

END 


$'end  of  data' 

' TAG ' = '  -  ACT. 
-  H202  (uM) ' 

$ ' end  of  H202DM' 


MITOCHONDRIA'  ... 


' - Heffetz  et  al.  (1990) - - - 

PROCED  FIG21A 

'effect  of  pervanadate  on  protein  Tyr  phosphatase' 

'in  the  presence  of  2  mM  H202' 

'Tyr  Pase  activity  measured  as  [32P]  remaining' 

'in  [32P] poly- (Glu, Tyr) ' 

SET  TITLE  =  'BBDR:  Effect  of  PerVanadate  on  PTyrPase’ 
PREPAR  ’ CLEAR' , Pconc, Inh 

SET  cmin=0 . 1,  cdelt=10 . ,  cmax=1000 . , TP=0 . 333, TSTOP=0 . 42 
SET  NRWITG=.T. , ki=18, kt=200, kd=l . , TME=0 . 333 


1  uM 
DATA 

%  Control 

Pconc 

Inh 

0.01 

1. 

1. 

0.9467 

10. 

0.6 

100. 

0.3933 

1000. 

0.1267 
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$,end  of  data1 


END 
START 

PLOT  Inh, 'lo'=0. , 'xlog' , 'xlo'=0. 01, 'TAG'='  -ACT.  PTPase'... 

' XTAG '  =  '  -  PERVAN  (uM) ' 

END  $ ' end  of  VANADTP' 

> - Vroegop,  et  al.  (1995) - ' 

PROCED  FIG21B 

'effect  of  CumOOH  on  mitochondria  staining  1  =  0.49  OD' 

SET  TITLE  =  ' BBDR :  Effect  of  Cum.OOH  on  Mitochondria' 

PREPAR  'CLEAR' , Pconc, Inh 

SET  cmin=10 . , cdelt=10 . , cmax=1000 . , TP=1 . ,  TSTOP=2 . 

SET  NRWITG=. T. , ki=100, kt=200, kd=0 . 05,  TME=4 . 

' uM  %  Control ' 

DATA 

Pconc  Inh 

0.  1. 

3.  1.0 

10.  1. 

33.  0.95 

100.  0.91 

330.  1.0 

1000.  0.79 

END  $'end  of  data' 

START 

PLOT  Inh, 'lo'=0.,  ' TAG ' = '  -ACT.  MITOCHONDRIA'  ... 

'xlog',  'xlo'=l, 'XTAG'  =  '  -  Log  CUM.OOH  (uM) ' 

END  .  $ 'end  of  CHDMIT ' 

PROCED  FIG22 

'effect  of  cumene  hydroperoxide  on  glucose' 

'transporter  1  =  18300  cpm' 

SET  TITLE  =  'BBDR:  Effect  of  Cum.OOH  on  Glucose  Transporter’ 
PREPAR  ' CLEAR ' , Pconc,  Inh 

SET  cmin=10 . , cdelt=2 . , cmax=200 . ,  TP=1 . , TSTOP=2 . 

SET  NRWITG=. T. , ki=100, kt=200, kd=0 . 1,  TME=1 . 

'uM  %  Control' 

DATA 

Pconc  Inh 

0.  1. 

4.5  0.97 

10.  0.90 

20.  0.82 

45.  0.62 

100.  0.46 

125.  0.38 

150.  0.30 

200.  0.19 

END  $'end  of  data' 

START 

PLOT  Inh, 'lo'=0.,  ' TAG ' = '  -ACT.  GLUC  TRANSPORTERS '  ... 

•XTAG'  =  '  -  CUM.OOH  (uM) ' 

END  $ ' end  of  CHDG' 

' - Theoretical  simulation - ' 

PROCED  FIG26A 

’theoretical  BBDR  deterministic  simulation  with  the  range  of  TCE ’ 
’ local  doses  corresponding  to  ethane  exhalation  experiments  ' 

SET  TITLE  =  ’BBDR:  Simulated  Effects  of  TCE  on  PTyrPase  Activity' 


PREPAR  'CLEAR* ,PCONC,  Inh 

SET  cmin=0. 5,  cdelt=0 .5,  cmax=50 . , TP=0 .  5, TSTOP=0 . 633, TME=0 . 5 
SET  NRWITG=.T.,  ki=900.,  kt=200.,  kd=l. 

START 

'SET  TME=1 . ' 

' START ’ 

' SET  TME=2 . ' 

'  START ' 

PLOT  Inh, 'char '='  ' , 'lo ' =0 . , ' tag'=' -  %  ACT.  PTPase'... 

'xhi'=50. , 'xtag'=' -  TCE  (um/g) ’ 

END  $ 'end  of  TCEPTP' 

*********************************************************************** 

BBDR  Stochastic  Array  Graphics  *.CSL  FILE 

PROGRAM:  PLOT  FOR  GRAFGAUS 

’Program  allows  to  plot  GRAPHG.CMD  files  created  by  GRAFGAUS  ’ 
'Final  version  8/28/95  ' 

INITIAL 

variable  N  =  0.  $' defines  independent  variable  N  =  PCONC  ' 

CONSTANT  Nmax=100.$ 'defines  the  maximum  PCONC  value  in  plot  ' 

CONSTANT  CINT=  1.  $' reports  at  every  whole  unit  of  PCONC  ' 

CONSTANT  FRAD=  0.  $ ' sham  value  for  another  independent  variable' 

END 

DYNAMIC 
PCONC  =  N 

PROB  =  0.  $'sham  values  for  dependent  variables  ' 

ACT  =  0. 

END 

TERMT  (N.GE.Nmax)  $'Stop  at  maximum  dose  ' 

END 

****************************************************************************** 


Output  of  stochastic  module,  created  by  *BBDR  Sub-Model: 
BBDR  Array  Graphics  *.CMD  FILE 


'FIG  19  A' 
PROC  NDATA 


prepar  pconc,prob 

set  title= ' BBDR:  Effect  of  Cum.OOH  on  AA  Transporter' 


set  Nmax=1000. 
DATA 

PCONC  PROB 

0.  1. 

10.000  0.9927 

20.000  0.9883 

30.000  0.9834 

40.000  0.9781 
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50.000 

0.9723 

60.000 

0.9659 

70.000 

0.9591 

80.000 

0.9517 

90.000 

0.9439 

100.000 

0.9356 

110.000 

0.9269 

120.000 

0.9177 

130.000 

0.9082 

140.000 

0.8982 

150.000 

0.8879 

160.000 

0.8772 

170.000 

0.8661 

180.000 

0.8548 

190.000 

0.8431 

200.000 

0.8312 

210.000 

0.8190 

220.000 

0.8066 

230.000 

0.7940 

240.000 

0.7813 

250.000 

0.7683 

260.000 

0.7552 

270.000 

0.7419 

280.000 

0.7286 

290.000 

0.7152 

300.000 

0.7017 

310.000 

0.6881 

320.000 

0.6745 

330.000 

0.6609 

340.000 

0.6472 

350.000 

0.6336 

360.000 

0.6200 

370.000 

0.6065 

380.000 

0.5930 

390.000 

0.5796 

400.000 

0.5662 

410.000 

0.5530 

420.000 

0.5398 

430.000 

0.5268 

440.000 

0.5139 

450.000 

0.5011 

460.000 

0.4885 

470.000 

0.4760 

480.000 

0.4636 

490.000 

0.4514 

500.000 

0.4394 

510.000 

0.4276 

520.000 

0.4159 

530.000 

0.4045 

540.000 

0.3932 

550.000 

0.3821 

560.000 

0.3712 

570.000 

0.3605 

580.000 

0.3500 
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590.000 

0.3397 

600.000 

0.3296 

610.000 

0.3197 

620.000 

0.3100 

630.000 

0.3005 

640.000 

0.2913 

650.000 

0.2822 

660.000 

0.2734 

670.000 

0.2647 

680.000 

0.2562 

690.000 

0.2480 

700.000 

0.2400 

710.000 

0.2321 

720.000 

0.2245 

730.000 

0.2170 

740.000 

0.2098 

750.000 

0.2027 

760.000 

0.1958 

770.000 

0.1891 

780.000 

0.1826 

790.000 

0.1763 

800.000 

0.1702 

810.000 

0.1642 

820.000 

0.1584 

830.000 

0.1528 

840.000 

0.1473 

850.000 

0.1420 

860.000 

0.1369 

870.000 

0.1319 

880.000 

0.1271 

890.000 

0.1224 

900.000 

0.1179 

910.000 

0.1135 

920.000 

0.1093 

930.000 

0.1052 

940.000 

0.1012 

950.000 

0.0974 

960.000 

0.0937 

970.000 

0.0901 

980.000 

0.0866 

990.000 

0.0833 

1000.000 

0.0800 

END  $ ' 

end  of  data 

START 

PLOT  PROB, ' lo '=0, ' hi '=1, 'xhi '=1000, 'xtag '='  -CUM.OOH  (uM) ' 
END  $ ' END  of  file' 

'FIG  19  B' 

PROC  NDATA 
prepar  pconc,prob 

set  title= ' BBDR:  Effect  of  H202  on  Glucose  Transporter' 
set  Nmax=10000. 

DATA 

PCONC  PROB 
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INITIAL 


0. 

1. 

100.000 

0.9927 

200.000 

0.9883 

300.000 

0.9834 

400.000 

0.9781 

500.000 

0.9723 

600.000 

0.9659 

700.000 

0.9591 

800.000 

0.9517 

900.000 

0.9439 

1000.000 

0.9356 

1100.000 

0.9269 

1200.000 

0.9177 

1300.000 

0.9082 

1400.000 

0.8982 

1500.000 

0.8879 

1600.000 

0.8772 

1700.000 

0.8661 

1800.000 

0.8548 

1900.000 

0.8431 

2000.000 

0.8312 

2100.000 

0.8190 

2200.000 

0.8066 

2300.000 

0.7940 

2400.000 

0.7813 

2500.000 

0.7683 

2600.000 

0.7552 

2700.000 

0.7419 

2800.000 

0.7286 

2900.000 

0.7152 

3000.000 

0.7016 

3100.000 

0.6881 

3200.000 

0.6745 

3300.000 

0.6609 

3400.000 

0.6472 

3500.000 

0.6336 

3600.000 

0.6200 

3700.000 

0.6065 

3800.000 

0.5930 

3900.000 

0.5796 

4000.000 

0.5662 

4100.000 

0.5530 

4200.000 

0.5398 

4300.000 

0.5268 

4400.000 

0.5139 

4500.000 

0.5011 

4600.000 

0.4885 

4700.000 

0.4760 

4800.000 

0.4636 

4900.000 

0.4514 

5000.000 

0.4394 

5100.000 

0.4276 

5200.000 

0.4159 

5300.000 

0.4045 
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5400.000 

0.3932 

5500.000 

0.3821 

5600.000 

0.3712 

5700.000 

0.3605 

5800.000 

0.3500 

5900.000 

0.3397 

6000.000 

0.3296 

6100.000 

0.3197 

6200.000 

0.3100 

6300.000 

0.3005 

6400.000 

0.2913 

6500.000 

0.2822 

6600.000 

0.2734 

6700.000 

0.2647 

6800.000 

0.2562 

6900.000 

0.2480 

7000.000 

0.2400 

7100.000 

0.2321 

7200.000 

0.2245 

7300.000 

0.2170 

7400.000 

0.2098 

7500.000 

0.2027 

7600.000 

0.1958 

7700.000 

0.1891 

7800.000 

0.1826 

7900.000 

0.1763 

8000.000 

0.1702 

8100.000 

0.1642 

8200.000 

0.1584 

8300.000 

0.1528 

8400.000 

0.1473 

8500.000 

0.1420 

8600.000 

0.1369 

8700.000 

0.1319 

8800.000 

0.1271 

8900.000 

0.1224 

9000.000 

0.1179 

9100.000 

0.1135 

9200.000 

0.1093 

9300.000 

0.1052 

9400.000 

0.1012 

9500.000 

0.0974 

9600.000 

0.0937 

9700.000 

0.0901 

9800.000 

0.0866 

9900.000 

0.0833 

10000.000 

0.0800 

END  $ ' End  of  data' 

START 

PLOT  PROB, ' lo'=0, ' hi'=l, ' xhi '=10000, . 

'xtag'='  -H202  (uM) ' 

END  $ ' END  of  file' 


82 


'FIG  19  C' 

PROC  NDATA 
prepar  pconc,prob 

set  title=’BBDR:  Effect  of  H202  on  PTyrPase' 
set  Nmax=500. 

DATA 


PCONC 

PROB 

0. 

1. 

0.010 

0.9959 

10.010 

0.9817 

20.010 

0.9680 

30.010 

0.9525 

40.010 

0.9353 

50.010 

0.9163 

60.010 

0.8958 

70.010 

0.8740 

80.010 

0.8509 

90.010 

0.8269 

100.010 

0.8020 

110.010 

0.7764 

120.010 

0.7503 

130.010 

0.7238 

140.010 

0.6970 

150.010 

0.6702 

160.010 

0.6434 

170.010 

0.6167 

180.010 

0.5902 

190.010 

0.5641 

200.010 

0.5384 

210.010 

0.5131 

220.010 

0.4884 

230.010 

0.4643 

240.010 

0.4409 

250.010 

0.4181 

260.010 

0.3960 

270.010 

0.3747 

280.010 

0.3541 

290.010 

0.3343 

300.010 

0.3153 

310.010 

0.2971 

320.010 

0.2796 

330.010 

0.2629 

340.010 

0.2470 

350.010 

0.2318 

360.010 

0.2174 

370.010 

0.2037 

380.010 

0.1906 

390.010 

0.1783 

400.010 

0.1666 

410.010 

0.1556 

420.010 

0.1452 

430.010 

0.1354 

440.010 

0.1261 

450.010 

0.1174 
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460.010  0.1093 

470.010  0.1016 

480.010  0.0944 

490.010  0.0876 

END  $ ' End  of  data' 

START 

PLOT  PROB, ' lo ' =0, ' hi '=1 . , ' xhi ' =500 , ' xtag '= '  -H202  (uM) ' 
END  $ ' End  of  file' 

'FIG  19  D' 

PROC  NDATA 
prepar  pconc,prob 

set  title='BBDR:  Effect  of  Vanad.  on  PTyrPase' 
set  Nmax=1000. 

DATA 

PCONC  PROB 

0.  1.  INITIAL 

0.100  0.9976 

10.100  0.9930 

20.100  0.9890 

30.100  0.9848 

40.100  0.9803 

50.100  0.9754 

60.100  0.9702 

70.100  0.9646 

80.100  0.9587 

90.100  0.9525 

100.100  0.9460 

110.100  0.9391 

120.100  0.9320 

130.100  0.9246 

140.100  0.9170 

150.100  0.9091 

160.100  0.9011 

170.100  0.8928 

180.100  0.8843 

190.100  0.8756 

200.100  0.8668 

210.100  0.8578 

220.100  0.8487 

230.100  0.8395 

240.100  0.8302 

250.100  0.8207 

260.100  0.8113 

270.100  0.8017 

280.100  0.7921 

290.100  0.7825 

300.100  0.7728 

310.100  0.7631 

320.100  0.7535 

330.100  0.7438 

340.100  0.7342 

350.100  0.7246 

360.100  0.7150 
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370.100 

0.7055 

380.100 

0.6960 

390.100 

0.6866 

400.100 

0.6773 

410.100 

0.6680 

420.100 

0.6588 

430.100 

0.6498 

440.100 

0.6408 

450.100 

0.6319 

460.100 

0.6232 

470.100 

0.6145 

480.100 

0.6060 

490.100 

0.5976 

500.100 

0.5893 

510.100 

0.5811 

520.100 

0.5731 

530.100 

0.5652 

540.100 

0.5574 

550.100 

0.5498 

560.100 

0.5423 

570.100 

0.5350 

580.100 

0.5278 

590.100 

0.5207 

600.100 

0.5138 

610.100 

0.5070 

620.100 

0.5004 

630.100 

0.4939 

640.100 

0.4875 

650.100 

0.4813 

660.100 

0.4753 

670.100 

0.4694 

680.100 

0.4636 

690.100 

0.4579 

700.100 

0.4524 

710.100 

0.4471 

720.100 

0.4419 

730.100 

0.4368 

740.100 

0.4318 

750.100 

0.4270 

760.100 

0.4223 

770.100 

0.4177 

780.100 

0.4132 

790.100 

0.4089 

800.100 

0.4047 

810.100 

0.4006 

820.100 

0.3967 

830.100  • 

0.3928 

840.100 

0.3891 

850.100 

0.3854 

860.100 

0.3819 

870.100 

0.3785 

880.100 

0.3752 

890.100 

0.3720 

900.100 

0.3689 

910.100  0.3658 

920.100  0.3629 

930.100  0.3601 

940.100  0.3574 

950.100  0.3547 

960.100  0.3521 

970.100  0.3497 

980.100  0.3473 

990.100  0.3450 

END  $ ' END  of  data' 

START 

PLOT  PROB, 'lo'=0, 'hi'=l, 'xtag’='  -VANAD  (uM) ' 

END  $ ' end  of  file' 

'FIG  20  Af 
PROC  NDATA 
Prepar  Pconc,Prob 

Set  title= ' BBDR:  Effect  of  6-OH  Dopamine  on  AA  Transporter' 
Set  Nmax=300. 


DATA 

PCONC 

PROB 

0.0 

1.  INITIAL 

10.000 

0.9849 

20.000 

0.9670 

30.000 

0.9437 

40.000 

0.9154 

50.000 

0.8825 

60.000 

0.8459 

70.000 

0.8060 

80.000 

0.7639 

90.000 

0.7200 

100.000 

0.6752 

110.000 

0.6301 

120.000 

0.5853 

130.000 

0.5412 

140.000 

0.4982 

150.000 

0.4568 

160.000 

0.4172 

170.000 

0.3796 

180.000 

0.3442 

190.000 

0.3109 

200.000 

0.2800 

210.000 

0.2513 

220.000 

0.2249 

230.000 

0.2007 

240.000 

0.1786 

250.000 

0.1585 

260.000 

0.1402 

270.000 

0.1238 

280.000 

0.1090 

290.000 

0.0958 

300.000 

0.0840 

END 

$'End  of  Data' 

START 
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PLOT  PROB, ' lo ' =0 , ' hi ' =1 . , ' tag ' = '  -Act.  AA  Transp ' . 
'xhi  '=300,  ' Xtag'  =  '  -  6-OH  DOPAMINE  (uM)  ' 

END  $ ' End  of  File' 

'FIG  20  B' 

PROC  NDATA 
prepar  pconc,prob 

set  title= ' BBDR:  Effect  of  H202  on  AA  Transporter’ 
set  nmax=100. 

DATA 

PCONC  PROB 

0.  1.  INITIAL 

10.000  0.9816 

11.000  0.9776 

12.000  0.9732 

13.000  0.9683 

14.000  0.9629 

15.000  0.9570 

16.000  0.9506 

17.000  0.9437 

18.000  0.9362 

19.000  0.9283 

20.000  0.9199 

21.000  0.9110 
22.000  0.9016 

23.000  0.8917 

24.000  0.8814 

25.000  0.8706 

26.000  0.8594 

27.000  0.8478 

28.000  0.8358 

29.000  0.8234 

30.000  0.8107 

31.000  0.7977 

32.000  0.7843 

33.000  0.7707 

34.000  0.7568 

35.000  0.7427 

36.000  0.7284 

37.000  0.7139 

38.000  0.6993 

39.000  0.6845 

40.000  0.6696 

41.000  0.6547 

42.000  0.6396 

43.000  0.6246 

44.000  0.6095 

45.000  0.5944 

46.000  0.5794 

47.000  0.5644 

48.000  0.5494 

49.000  0.5346 

50.000  0.5198 

51.000  0.5052 
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52.000 

0.4907 

53.000 

0.4764 

54.000 

0.4622 

55.000 

0.4482 

56.000 

0.4344 

57.000 

0.4208 

58.000 

0.4074 

59.000 

0.3943 

60.000 

0.3813 

61.000 

0.3686 

62.000 

0.3562 

63.000 

0.3440 

64.000 

0.3320 

65.000 

0.3203 

66.000 

0.3089 

67.000 

0.2977 

68.000 

0.2868 

69.000 

0.2762 

70.000 

0.2658 

71.000 

0.2557 

72.000 

0.2459 

73.000 

0.2364 

74.000 

0.2271 

75.000 

0.2181 

76.000 

0.2094 

77.000 

0.2009 

78.000 

0.1927 

79.000 

0.1848 

80.000 

0.1771 

81.000 

0.1696 

82.000 

0.1624 

83.000 

0.1555 

84.000 

0.1488 

85.000 

0.1423 

86.000 

0.1360 

87.000 

0.1300 

88.000 

0.1242 

89.000 

0.1186 

90.000 

0.1133 

91.000 

0.1081 

92.000 

0.1031 

93.000 

0.0983 

94.000 

0.0937 

95.000 

0.0893 

96.000 

0.0851 

97.000 

0.0810 

98.000 

0.0772 

99.000 

0.0734 

100.000 

0.0699 

END  $'End  of  data' 
START 

PLOT  PROB, 'lo'=0, 'hi' 
END  $ ' END  of  file' 


'FIG  20  C' 

PROC  NDATA 
prepar  pconc,  prob 

set  title= ' BBDR:  Effect  of  6-OH  Dopamine  on  Mitochondria' 
SET  NMAX=200 . 


DATA 


PCONC 

PROB 

0. 

1. 

40.000 

0.9826 

50.000 

0.9672 

60.000 

0.9455 

70.000 

0.9178 

80.000 

0.8852 

90.000 

0.8489 

100.000 

0.8106 

110.000 

0.7720 

120.000 

0.7344 

130.000 

0.6990 

140.000 

0.6667 

150.000 

0.6379 

160.000 

0. 6129 

170.000 

0.5917 

180.000 

0.5740 

190.000 

0.5595 

200.000 

0.5478 

END 

START 

PLOT  PROB, 'lo'=0. 
'xhi '=200. , 


END 


INITIAL 


, 'hi'=l., 'tag'='  -ACT  Mi to'... 
'xtag'='  -6-OH  DOPAMINE  (uM) ' 


'FIG  20  D' 

PROC  NDATA 
prepar  pconc, prob 

set  title= ' BBDR:  Effect  of  H202  on  Mitochondria' 


set  Nmax=100. 
DATA 
PCONC 


PROB 


0. 

1.  INITIAL 

10.000 

0.9608 

20.000 

0.8800 

30.000 

0.7663 

40.000 

0.6362 

50.000 

0.5061 

60.000 

0.3874 

70.000 

0.2867 

80.000 

0.2059 

90.000 

0.1439 

100.000 

0.0982 

END 

$ ' end  of  data ' 

START 

PLOT  PROB, ' lo'=0 . , 'hi'= 

END  $ 

'END  of  file' 

'xhi'=100, 'xtag'='  -H202 


(uM)  ' 
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'FIG  21  A' 

PROC  NDATA 
prepar  pconc,prob 

set  title= ' BBDR:  Effect  of  PerVanadate  on  PTyrPase' 

set  Nmax=1000 

DATA 

PCONC  PROB 

0.01  1.  INITIAL 

0.100  0.9966 

10.100  0.9902 

20.100  0.9846 

30.100  0.9788 

40.100  0.9724 

50.100  0.9656 

60.100  0.9583 

70.100  0.9505 

80.100  0.9422 

90.100  0.9335 

100.100  0.9244 

110.100  0.9148 

120.100  0.9048 

130.100  0.8945 

140.100  0.8838 

150.100  0.8728 

160.100  0.8615 

170.100  0.8499 

180.100  0.8380 

190.100  0.8258 

200.100  0.8135 

210.100  0.8009 

220.100  0.7882 

230.100  0.7753 

240.100  0.7622 

250.100  0.7490 

260.100  0.7358 

270.100  0.7224 

280.100  0.7090 

290.100  0.6955 

300.100  0.6819 

310.100  0.6684 

320.100  0.6549 

330.100  0.6413 

340.100  0.6278 

350.100  0.6144 

360.100  0.6010 

370.100  0.5877 

380.100  0.5744 

390.100  0.5612 

400.100  0.5482 

410.100  0.5352 

420.100  0.5224 

430.100  0.5097 

440.100  0.4971 

450.100  0.4847 
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460.100 

0.4724 

470.100 

0.4603 

480.100 

0.4484 

490.100 

0.4366 

500.100 

0.4250 

510.100 

0.4135 

520.100 

0.4023 

530.100 

0.3912 

540.100 

0.3804 

550.100 

0.3697 

560.100 

0.3592 

570.100 

0.3489 

580.100 

0.3389 

590.100 

0.3290 

600.100 

0.3193 

610.100 

0.3098 

620.100 

0.3005 

630.100 

0.2914 

640.100 

0.2826 

650.100 

0.2739 

660.100 

0.2654 

670.100 

0.2571 

680.100 

0.2490 

690.100 

0.2411 

700.100 

0.2334 

710.100 

0.2259 

720.100 

0.2186 

730.100 

0.2115 

740.100 

0.2045 

750.100 

0.1978 

760.100 

0.1912 

770.100 

0.1848 

780.100 

0.1785 

790.100 

0.1725 

800.100 

0 . 1666 

810.100 

0.1609 

820.100 

0.1553 

830.100 

0.1499 

840.100 

0.1447 

850.100 

0.1396 

860.100 

0.1347 

870.100 

0.1299 

880.100 

0.1253 

890.100 

0.1208 

900.100 

0.1164 

910.100 

0.1122 

920.100 

0.1081 

930.100 

0.1041 

940.100 

0.1003 

950.100 

0.0966 

960.100 

0.0930 

970.100 

0.0895 

980.100 

0.0862 

990.100 

0.0829 

END  $ ' END  of  data’ 

START 

PLOT  PROB, ' lo'=0 . , 'hi'=l, 'xlog', ’xlo'=0.01,  .  .  . 

' xtag ' = '  -Log  PerVan  (uM) ' 

END  $'End  of  file' 

'FIG  21  B' 

PROC  NDATA 
prepar  pconc,prob 

set  title= ’ BBDR:  Effect  of  Cum.OOH  on  Mitochondria’ 
set  Nmax=1000. 

DATA 

PCONC  PROB 

0.1  1.  INITIAL 

10.000  0.9982 

20.000  0.9971 

30.000  0.9959 

40.000  0.9945 

50.000  0.9931 

60.000  0.9915 

70.000  0.9898 

80.000  0.9879 

90.000  0.9860 

100.000  0.9839 

110.000  0.9817 

120.000  0.9794 

130.000  0.9770 

140.000  0.9745 

150.000  0.9720 

160.000  0.9693 

170.000  0.9665 

180.000  0.9637 

190.000  0.9608 

200.000  0.9578 

210.000  0.9548 

220.000  0.9517 

230.000  0.9485 

240.000  0.9453 

250.000  0.9421 

260.000  0.9388 

270.000  0.9355 

280.000  0.9321 

290.000  0.9288 

300.000  0.9254 

310.000  0.9220 

320.000  0.9186 

330.000  0.9152 

340.000  0.9118 

350.000  0.9084 

360.000  0.9050 

370.000  0.9016 

380.000  0.8983 

390.000  0.8949 

400.000  0.8916 
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410.000 

0.8882 

420.000 

0.8850 

430.000 

0.8817 

440.000 

0.8785 

450.000 

0.8753 

460.000 

0.8721 

470.000 

0.8690 

480.000 

0.8659 

490.000 

0.8629 

500.000 

0.8599 

510.000 

0.8569 

520.000 

0.8540 

530.000 

0.8511 

540.000 

0.8483 

550.000 

0.8455 

560.000 

0.8428 

570.000 

0.8401 

580.000 

0.8375 

590.000 

0.8349 

600.000 

0.8324 

610.000 

0.8299 

620.000 

0.8275 

630.000 

0.8251 

640.000 

0.8228 

650.000 

0.8206 

660.000 

0.8183 

670.000 

0.8162 

680.000 

0.8141 

690.000 

0.8120 

700.000 

0.8100 

710.000 

0.8080 

720.000 

0.8061 

730.000 

0.8043 

740.000 

0.8024 

750.000 

0.8007 

760.000 

0.7990 

770.000 

0.7973 

780.000 

0.7957 

790.000 

0.7941 

800.000 

0.7925 

810.000 

0.7911 

820.000 

0.7896 

830.000 

0.7882 

840.000 

0.7868 

850.000 

0.7855 

860.000 

0.7842 

870.000 

0.7830 

880.000 

0.7818 

890.000 

0.7806 

900.000 

0.7795 

910.000 

0.7784 

920.000 

0.7773 

930.000 

0.7763 

940.000 

0.7753 
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950.000  0.7743 

960.000  0.7734 

970.000  0.7725 

980.000  0.7717 

990.000  0.7708 

1000.000  0.7700 

END  $ ' END  of  data' 

START 

PLOT  PROB, ' lo ' =0 . , ' hi ' =1 . , 'xlog ' , 'xlo ' =1 . ,  .  . . 

'xtag'='  -Log  Cum.OOH  (uM)  ' 

END  $ ' End  of  file' 

'FIG  22' 

PROC  NDATA 
prepar  pconc,prob 

set  title= ' BBDR:  Effect  of  Cum.OOH  on  Glucose  Transporter' 
set  Nmax=200. 

DATA 

PCONC  PROB 

0.  1.  INITIAL 

10.000  0.9792 

20.000  0.9487 

30.000  0.9068 

40.000  0.8553 

50.000  0.7963 

60.000  0.7323 

70.000  0.6656 

80.000  0.5984 

90.000  0.5325 

100.000  0.4692 

110.000  0.4098 

120.000  0.3548 

130.000  0.3048 

140.000  0.2598 

150.000  0.2200 

160.000  0.1850 

170.000  0.1546 

180.000  0.1284 

190.000  0.1061 
200.000  0.0872 

END  $ ' END  of  data' 

START 

PLOT  PROB, ' lo ' =0 , ' hi '=1 . , 'xtag ' ='  -Cum.OOH  (uM) ' 

END  $ ' END  of  file' 

'FIG  26  B' 

PROC  NDATA 
PREPAR  pconc,prob 

set  title= ' BBDR:  Simulated  Effects  of  TCE  on  PTyrPase  Activity' 

set  Nmax=50 

DATA 

PCONC  PROB 

0.500  0.9927 

1.000  0.9883 
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1.500 

0.9834 

2.000 

0.9781 

2.500 

0.9723 

3.000 

0.9659 

3.500 

0.9591 

4.000 

0.9517 

4.500 

0.9439 

5.000 

0.9356 

5.500 

0.9269 

6.000 

0.9177 

6.500 

0.9082 

7.000 

0.8982 

7.500 

0.8879 

8.000 

0.8772 

8.500 

0.8661 

9.000 

0.8548 

9.500 

0.8431 

10.000 

0.8312 

10.500 

0.8190 

11.000 

0.8066 

11.500 

0.7940 

12.000 

0.7813 

12.500 

0.7683 

13.000 

0.7552 

13.500 

0.7419 

14.000 

0.7286 

14.500 

0.7152 

15.000 

0.7016 

15.500 

0.6881 

16.000 

0.6745 

16.500 

0.6609 

17.000 

0.6472 

17.500 

0.6336 

18.000 

0.6200 

18.500 

0.6065 

19.000 

0.5930 

19.500 

0.5796 

20.000 

0.5662 

20.500 

0.5530 

21.000 

0.5398 

21.500 

0.5268 

22.000 

0.5139 

22.500 

0.5011 

23.000 

0.4885 

23.500 

0.4760 

24.000 

0.4636 

24.500 

0.4514 

25.000 

0.4394 

25.500 

0.4276 

26.000 

0.4159 

26.500 

0.4045 

27.000 

0.3932 

27.500 

0.3821 

28.000 

0.3712 

95 


28.500 

0.3605 

29.000 

0.3500 

29.500 

0.3397 

30.000 

0.3296 

30.500 

0.3197 

31.000 

0.3100 

31.500 

0.3005 

32.000 

0.2913 

32.500 

0.2822 

33.000 

0.2734 

33.500 

0.2647 

34.000 

0.2562 

34.500 

0.2480 

35.000 

0.2400 

35.500 

0.2321 

36.000 

0.2245 

36.500 

0.2170 

37.000 

0.2098 

37.500 

0.2027 

38.000 

0.1958 

38.500 

0.1891 

39.000 

0.1826 

39.500 

0.1763 

40.000 

0.1702 

40.500 

0.1642 

41.000  • 

0.1584 

41.500 

0.1528 

42.000 

0.1473 

42.500 

0.1420 

43.000 

0.1369 

43.500 

0.1319 

44.000 

0.1271 

44.500 

0.1224 

45.000 

0.1179 

45.500 

0.1135 

46.000 

0.1093 

46.500 

0.1052 

47.000 

0.1012 

47.500 

0.0974 

48.000 

0.0937 

48.500 

0.0901 

49.000 

0.0866 

49.500 

0.0833 

50.000 

0.0800 

END  $ ' End  of  data' 

START 

PLOT  PROB, ' lo ' =0 . ,  ' hi ' =1 . , ' xi ' =500 , ' xtag ' = ' -TCE  (umol/g  liver)' 
END  $'End  of  file' 
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